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THERMAL NEUTRON CAPTURE CROSS SECTION OF Ce'*! 
By L. P. Roy? AND L. YAFFE 


ABSTRACT 


The cross section of Ce? has been determined by the activation method and 
found to be 0.95+0.05 barns relative to a value of oc.5* = 36.3 barns. Disintegra- 
tion rates were determined by a 4x counter. The value obtained agrees with that 
of Katcoff et al. (1949) also obtained by the activation method but differs 
markedly from that of Pomerance (1952).* 


INTRODUCTION 


In the course of some related work (11) it became necessary to determine 
the capture cross section of Ce'*® as accurately as possible. Bothe (1), by the 
activation method, had obtained a value of 1.1 barns relative to a value of 
49 barns for the capture cross section o of Ho'®. Using the most recent value 
of oyo'% = 64 barns (7, 12), then ace = 1.4 barns according to Bothe. 
Katcoff et al. (5) obtained a value of 0.95+0.18 barns by the activation 
method using uranium metal as a flux monitor. Pomerance (8) used the pile 
oscillator technique to obtain a value of 1.8+0.2 barns for isotopically en- 
riched Ce!42. He used o,,y!*7 = 95 barns as a standard. Using the more recent 
value of 98.7 barns (2) Pomerance’s value becomes 1.9 barns. The spread in 
these values was too great to be used in our work and large enough to warrant 
a redetermination of o¢,e42 by the activation method using improved tech- 
n.ques. 

If naturally-occurring cerium were irradiated in a nuclear reactor the follow- 
ing first-order reactions would occur: 


(1) ssCe8* (0.193% of natural cerium) 
(n, 7) 


ssCe!37 





‘ Ex. 
> 57La!3? ————> ,,Ba!3? (stable) 
hr. > 400 yr. 


1Manuscript received March 28, 1956. 

Contribution from the Radiochemistry Laboratory, Department of Chemistry, McGill University, 
Montreal, with financial assistance from the National Research Council of Canada and the Atomic 
Energy Control Board. 

2Holder of a Shawinigan Chemical Scholarship (1952-53), a National Research Council Student- 
ship (1958-54), and a Province of Quebec Bursary (1954-55). Present address: Atomic Energy 
of Canada Limited, Chalk River, Ont. 

*Note ADDED IN PRroor.—Dr. Pomerance his been kind enough to inform us that his sample 
of enriched Ce? showed a detectable amount (<1500 p.p.m.) of cadmium. His calculations show 
that, if the sample contained exactly 1500 p.p.m. of cadmium, his value for the cross section of 
Ce'42 would be 1.040.2 barns. 
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(2) ssCe!8 (0.250% of natural cerium) 
(n, v) 


ssCce!39 ae s7La'®® (stable) 


(3) ssCe'*? (88.47% of natural cerium) 
(n, y) : 


ssCe!4! =" sgPr'4! (stable) 


(4) ssCe'#? (11.07% of natural cerium) 
(n, ¥) 


sacel43 





> g9Pri43 oe Nd!43 (stable) 
33 hr. 13.5 d. 

Owing to the multiplicity of radioactive cerium isotopes formed it is ob- 
viously very difficult to measure Ce!** directly. However it is possible to deter- 
mine the activity of the Ce!** produced by separating and measuring its 13.5 
day Pr'43 daughter. 

The differential equations involved are well known and it can be shown 
that the activity of Pr'** at any time ¢, after the irradiation is given by 


Arde 


x, Brle™ 1 gM) Cedee 2" , 
ae 


[1] Cy, = 


where Bry and Cr are the numbers of atoms of Ce!*? and Pr'43 respectively 
formed after a time of irradiation T and 


(2 By = MAS (ye?) 

and 

[3] Cr= mates (1—e™) —"eten (rng). 
2 1~ A2 


A is the number of atoms of Ce!*? present during the irradiation in a neutron 
flux of mv neutrons/cm.?/sec. and is assumed to be constant. oo, o; and Aj, A> 
are cross sections and disintegration constants respectively, the subscripts 
0, 1, 2 referring to Ce'*?, Ce!48, and Pr'** respectively. 


EXPERIMENTAL 


Samples of spectroscopically pure cerium oxide were sealed into small 
quartz ampoules and irradiated in a position in the NRX Chalk River heavy 
water reactor where the neutron spectrum was predominantly thermal. 
With each sample was enclosed a piece of fine weighed cobalt wire which was 
used as a flux monitor. 

After irradiation the cerium oxide was dissolved in boiling concentrated 
sulphuric acid to which a few drops of 30% hydrogen peroxide was added. 
Ten milligrams of praseodymium carrier was added and the rare earths 
precipitated by the cautious addition of ammonium hydroxide. The precipi- 





ROY AND YAFFE: CROSS SECTION ; 1025 


tate was dissolved in nitric acid and cerium separated by oxidation with solid 
potassium bromate and precipitation as ceric iodate. The supernatant, con- 
taining the praseodymium, was neutralized by cautious addition of con- 
centrated sodium hydroxide which precipitated praseodymium hydroxide. 
The precipitate was dissolved, cerium carrier was added, and the entire 
cycle of cerium removal and praseodymium precipitation was repeated six 
or seven times. The final precipitate of praseodymium hydroxide was dissolved 
in 1 ml. of 6 N hydrochloric acid and made up to standard volume. The 
absolute disintegration rate was determined by the use of a 47 counter as 
described by Pate and Yaffe (6). 

The cobalt wire monitor was dissolved in 1 ml. aqua regia, the solution 
made to standard volume, and its activity determined in a similar manner. 

The chemical yield of the praseodymium was determined by precipitating 
it from the master solution as praseodymium oxalate and drying it to 
Pre(C2O.)3 . 9H.O. 

RESULTS AND DISCUSSION 


The Pr!43 was identified by its half-life. A typical decay curve is shown in 
Fig. 1. Values obtained from this and other samples which were followed 


°° 





= 





| oe 


COUNTING RATE (C/M) 





3, 








Db 66 610 61S 620 25 30 35 40 «45 
TIME (DAYS) 


Fic. 1. Decay curve of Pr’. 


for about three half-lives ranged between 13.5 and 13.7 days. As can be seen 
from the previous discussion the only activity which could interfere from the 
cerium irradiation would be La!*? which would not be separated from the 
praseodymium by the chemical procedures employed. A simple calculation 
shows that, owing to its long half-life, the La'*’ contribution to the activity 
will be negligible. Freedom from other rare earth radioactive impurities is 
attested to by the fact that the half-life obtained is very close to the currently 
accepted value (3, 4, 9, 10). 

Table I shows the flux values obtained from the cobalt monitor. o¢,% was 
taken as 36.3 (13) after correction for the more recent value of o,4y*7 = 98.7 
barns (2). 
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TABLE I 
NEUTRON FLUX VALUES OBTAINED FROM COBALT MONITOR 











Sample No. Co*® atoms, Time, Co®, Calc. flux, 
x< 107-19 sec. X10-* — dis./sec.K10-5 n./cm.?/sec. K107 
1 1.17 3.70 6.1 9.24 
2 4 37 7.19 3.6 2.80 





In Table II are shown the values of the cross section of Ce!#? as calculated 
by means of equation [1]. 











TABLE II 
VALUES OF CROSS SECTION OF Ce! OBTAINED 
Sample No. A, h, AsCi,, oo, 

Ce'? atoms X 10718 sec. X1075 dis./sec. X 1075 barns 
la 8.24 1.214 8.44 0.945 
1b 8.24 1.214 8.46 0.948 
lc 8.24 1.214 8.44 0.945 
2a 9.31 2.760 7.98 0.965 
2b 9.31 2.760 7.94 0.961 
2c 9.31 2.760 7.98 0.965 





An average value of 0.95+0.05 barns was obtained, in excellent agreement 
with the value of Katcoff et al. (5). The error quoted is five times the standard 
deviation obtained and certainly takes into account systematic errors. It is 
difficult to understand the higher value obtained by Pomerance.* Two possi- 
bilities occur to us. First, the cerium sample used by him might have contained 
small amounts of high neutron-absorbing rare earth impurities. Secondly 
Ce!43 may form a short-lived isomer. The cross section for its formation would 
not be detected by our technique but would be measured by the pile oscillator 
method. 
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GRAFT COPOLYMERS OF STYRENE AND METHYL 
METHACRYLATE 


PART II: VISCOSITY BEHAVIOR}? 


By M. H. Jones 


ABSTRACT 


A study of the viscosity behavior of the graft copolymers described in Part I 
has been made with dilute solutions in benzene at 25°C. Although the slope 
constants of the Huggins equation increase with the frequency of branching when 
measurements are made in a capillary viscometer under ‘free fall’ conditions, this 
is shown to be attributable to the dependence of viscosity on shear gradient. 
At a constant shear rate, the Huggins k’ values approximate to those of linear 
polymers. It is suggested that the marked increase in viscosity observed with 
decreasing shear rate for the graft copolymers is due to molecular entanglement. 


INTRODUCTION 


Although the viscosity behavior of dilute polymer solutions has been 
studied actively for a number of years, most of the experimental evidence is 
concerned with polymers which are essentially linear in structure and little 
is known of the influence of branching. 

Through a modification of the Stokes equation to take into consideration 
interactions between solvent and polymer molecules, Huggins (7) has derived 
the expression 


Nep/C = [n]+k'[nP? C+’ [n}bC?+ ... 


relating the reduced viscosity to the concentration of the solution. Fo" 
sufficiently dilute solutions, only the first two terms of the series are significant. 
The interaction coefficient, k’, is considered to be a constant for a given solvent— 
polymer system and, in general, is independent of molecular weight. This is 
similar to the empirical equation of Schulz and Blaschke (12). Increases in 
the slope constant for high molecular weight fractions from polymerization 
carried to high conversions, where the rate of transfer to polymer resulting in 
the formation of branched molecules is more significant, have been attributed 
to branching (4, 5, 11). A similar explanation has been given to observed 
increases in k’ for soluble fractions of styrene — divinyl benzene copolymers 
(10, 17) where branching occurs through the unreacted double bonds of the 
diolefinic units in the polymer chain. Theoretical calculations (14) have 
supported this contention. However, the results of Thurmond and Zimm (16) 
show no detectable difference in the slope constant of the Huggins equation 
between fractions of a styrene-divinyl benzene copolymer and a poly- 
styrene prepared under identical conditions. Decreases in the exponent of the 
Houwink viscosity — molecular weight relationship for high conversion poly- 
butadienes (8) and soluble crosslinked polystyrenes (8, 16) have also been 
attributed to branching. 


1Manuscript received March 21, 1956. 

Contribution from the Department of Chemistry, Ontario Research Foundation, Toronto 5, 
Ontario. 

2Presented in part at the 6th Canadian High Polymer Forum held at St. Catharines, Ontario, 
April 14-15, 1956. 
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The series of well-characterized branched molecules described in Part I 
provide a means of studying some of these effects and observing their change 
with progressive increments of branching. 


EXPERIMENTAL 
Materials 
Samples of homopolymers and random copolymers of styrene and methyl 
methacrylate were prepared by bulk polymerization using a,a’-azo-bis- 
isobutyronitrile as initiator. The preparative conditions are recorded in 
Table I. Polymerizations were limited to small conversions at comparatively 


TABLE I 
LINEAR POLYMERS AND COPOLYMERS OF STYRENE AND METHYL METHACRYLATE 











Polymer Poly" Temp. [Initiator] Mole fraction % Poly" % PMMA 
CC.) (gm. 1.~') MMA in copolymer 

PS1 70 12.5 

PS 2 45 0.98 4.8 

PMMA 1 60 1.87 9.3 

PMMA 2 60 0.28 7.4 

Cop 1 60 0.30 0.811 fe 71 

Cop 2 60 0.30 0.682 6.0 60 

Cop 3 60 0.30 0.508 5.2 50 

Cop 4 60 0.093 0.682 5.2 58 





low temperatures; these conditions are such that branching may be con- 
sidered negligible. The compositions of the copolymers were determined by 
infrared analysis. The method of synthesis of the graft copolymers has been 
described in detail in Part I. 


Viscosity Measurements 

Viscosity measurements were made in benzene solution at 25.0°C. +0.005°C. 
using Ubbelohde and Ostwald type viscometers. A large bulb was incorporated 
in the Ubbelohde viscometer to permit dilution of the solution im situ. The 
Ostwald viscometer was similar to that used by Flory and his co-workers (6) 
and was designed to study the influence of shear rate on viscosity in the range 
encountered with the Ubbelohde viscometer. Kinetic energy corrections were 
calculated from the dimensions of the viscometers using the Wilberforce 
modification of the Poiseuille equation (1), i.e. 


ms, ar'hgpt _ ‘mpv 
a 8h Salt 


or 7= Apt—22. 


A value of 1.1 (1) was chosen for the constant m and the corrections were 
applied to the flow times. Hence the corrected flow time 


t* = ;-2 


At’ 








JONES: GRAFT COPOLYMERS. II . 1029 


For the Ubbelohde and Ostwald instruments, the maximum kinetic energy 
corrections for the solvent were 0.9% and 0.5% respectively. As these are small, 
the arbitrary choice of a value for the empirical constant, m, was considered 
justifiable. Flow times were determined with a stop watch graduated in 1/10 
sec. intervals and the reproducibility of successive readings was better than 
+0.1%. The dimensions of the two viscometers are given in Table II together 
with the corrected flow times for the solvent. 


TABLE II 
VISCOMETER CONSTANTS 











Viscometer v r l ey 
(cc.) (cm.) (cm.) (sec.) 
Ubbelohde 2.20 0.0230 13.3 115.1 
Ostwald: Upper bulb 2.37 0.0217 19.6 145 .9f 
Lower bulb 2.74 0.0217 19.6 222 .3T 





tFor solvent volume of 10.01 cc. 


Solution viscosities at different shear rates for a given polymer concentration 
were obtained by adding known volumes of the solution to the Ostwald 
viscometer. The maximum rates of shear, G, were derived from the expression 


G= rhgp/2lnn, 


where 7 is the viscosity of the solvent and 7, the relative viscosity of the 
solution. The hydrostatic head, h, of the liquid column in the viscometer was 
taken as the arithmetic mean of the maximum and minimum heights. Al- 
though this differs from the true mean as calculated by a formula of the 
Meissner type (1), the maximum difference for the present viscometer was 
less than 1.5%. 

Polymer solutions were prepared by shaking known weights of polymer 
with the solvent for 18 hr. to ensure complete solution. The solutions were 
filtered through sintered glass disks of medium porosity before transference 
to the viscometer. The usual precautions were taken to avoid contamination 
of the viscometers and solutions with dust particles. In general, reduced 
viscosities were determined at four or more concentrations in the range 
0.05-0.40 gm./100 cc. and the parameters of the Huggins equation obtained 
by the least squares method of calculation. 


RESULTS AND DISCUSSION 


The parameters of the Huggins equation for the series of graft copolymers 
determined with the Ubbelohde viscometer under ‘free fall’ conditions are 
given in Table III. Similar data on the linear polystyrenes, polymethyl 
methacrylates, and random copolymers are included for comparison. Usually 
duplicate determinations were made and the slope constants agreed within 
0.02 unit and the intrinsic viscosities within 1%. 

The k’ values for the linear polymers are constant at approximately 0.35. 
Even with small amounts of branching, i.e. one to two branches per molecule 
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TABLE III 
HUGGINS PARAMETERS FOR LINEAR AND GRAFT COPOLYMERS—VARIABLE SHEAR GRADIENT 











Polymer P, Fr. Average no. [n]es°c. k’ 
(backbone) (branches) branches per (dl. gm.~) 
molecule 

Al 2500 4500 1.3 3.26 0.40; 
A2 2500 1500 1.5 1.85 0.38; 
D1 3100 2300 2.6 3.10 0.545 
D2 3100 1700 3.4 3.10 0.587 
D3 3100 1600 3.0 3.17 0.61; 
D4 3100 1300 3.3 2.62 0.565 
El 3800 1500 1.8 2.99 0.43; 
E2 3800 1500 3.2 3.43 0.542 
E3 3800 1800 3.9 3.66 0.634 
E4 3800 1900 5.9 3.99 0.780 
PS 1 3.16 0.335 
FS3 1.73 0.358 
PMMA 1 1.84 0.35; 
PMMA 2 2.66 0.32, 
Cop 1 1.84 0.35: 
Cop 2 1.66 0.366 
Cop 3 1.57 0.345 
Cop 4 2.53 0.34: 





(Polymers Al, A2, and El) there appears to be a significant increase in k’ 
and, with more highly branched structures, the effect is pronounced. The 
progressive increase in slope constant with frequency of branching is most 
clearly demonstrated in the series E1-E4 where the variation in chain length 
is small compared to the increase in the number of branches. The reduced 
viscosity — concentration curves for these polymers are illustrated in Fig. 1. 
The broken lines indicate the slopes that would be obtained at the pertinent 
intrinsic viscosities for a k’ value of 0.35. Although these results were obtained 
with the same viscometer and consequently at a constant shear stress, the 
reduced viscosities at each concentration during dilution are measured at 
progressively increasing shear gradients since these are an inverse function 
of the solution viscosity. In general, the shear gradients for a viscosity plot 
varied in the range 500-1500 sec.—!. 

Studies of the influence of molecular weight and solvent on the shear de- 
pendence of viscosity for neutral polymers have been reported by several 
authors (3, 6, 9, 13). Their results indicate that shear effects may just become 
significant at molecular weights of the same order as those of the present 
series of polymers and an investigation of the influence of shear rate was 
made to determine the magnitude of the effect. Although the range of shear 
gradients covered was small, a considerable difference in behavior was ob- 
served between the linear and branched molecules. Solutions of the poly- 
styrene, polymethyl methacrylate, and the random copolymers, which were 
believed to have essentially no branching, behaved as Newtonian liquids 
(Figs. 2 and 3). The graft copolymers showed a marked dependence of viscosity 
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Fic. 1. Reduced viscosity versus concentration plots for graft copolymers E1-E4 with 
Ubbelohde viscometer. , 


on shear rate and this was clearly evident for samples having viscosities, and 
presumably molecular weights, similar to the. linear polymers. Reduced 
viscosity — shear plots for the graft copolymers with the greatest difference in 
intrinsic viscosity are shown in Figs. 4 and 5. The curvature of the plots 
precludes the usual practice of extrapolating to zero shear gradient to obtain 
reduced viscosities that are independent of shear effects, and values for each 
concentration were taken from the curves at an arbitrary shear rate of 500 
sec.—'. Huggins plots of these values for two of the graft copolymers (E3 and 
E4) are illustrated in Fig. 6 together with the curves obtained under ‘free 
fall’ conditions with the same viscometer. For the polymer E3 the curves have 
been displaced downwards 0.5 unit in the ordinate scale for clarity. The values 
of the Huggins parameters for these two polymers and for a linear polystyrene 
and polymethyl methacrylate are recorded in Table IV. The significant result 
is that the slope constants for the graft copolymers, when measured at a 
constant shear rate, decrease to approximately the value obtained for poly- 
mers having a linear structure. This, again, emphasizes the necessity of 
determining shear dependence in capillary viscometry of high polymer solu- 
tions and demonstrates the fact that non-Newtonian behavior may appear 
as a unique viscosity—concentration effect in a relationship of the Huggins 
type when results are obtained under ‘free fall’ conditions. It is obvious that 
the increases in slope constant for the graft copolymers obtained with the 
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Fic. 4. Reduced viscosity versus shear rate plots for various concentrations of polymer E4. 
Fic. 5. Reduced viscosity versus shear rate plots for various concentrations of polymer D4. 
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Fic. 6. Reduced viscosity versus concentration plots for graft copolymers E3 and E4 
with Ostwald viscometer. Ostwald upper bulb—“‘free fall’, -@-—; Ostwald lower bulb—‘free fall’, 
-O-; Ostwald, constant G = 500 sec.~!, -O-. 


Ubbelohde viscometer (Table II) are fortuitous. Although the data on the 
shear dependence of viscosity are limited, slope constants have also been 
determined at constant shear rates of 250, 750, and 1000 sec.—' where errors 
due to extrapolation should not be too serious. The values obtained for the 
polymer E3 are, respectively 0.37, 0.34, and 0.34, and for the polymer E4, 
0.39, 0.33, and 0.32. Although there appears to be a slight trend to increasing 
values of k’ with decrease in shear rate, the variation is probably within the 
experimental accuracy. 

It is to be expected that deviations from Newtonian behavior will be deter- 
mined by the ease of deformation and orientation of the molecules in solution 
and, consequently, will be related to the molecular weight and ‘goodness’ of 
the solvent since these factors govern the extension and softness of the expanded 
coil. However, it must be assumed that branched molecules, with their greater 
density of packing, will be less subject to deformation under flow than linear 
molecules but the observations show that the latter are independent of shear 
rate. Probably the explanation lies in the greater probability of molecular 
entanglement with the branched polymers because of the increased number of 
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TABLE IV 
INFLUENCE OF SHEAR RATE ON*HUGGINS PARAMETERS FOR LINEAR AND GRAFT COPOLYMERS 











Viscometer Shear rate [nlos°c. k’ 
(sec.—!) 
Polymer E3 
Ubbelohde 520-1380 3.66 0.63, 
Ostwald: Upper bulb 610-1750 3.84 0.533 
Lower bulb 470-1380 3.89 0.55; 
500 4.42 0.350 
Polymer E4 
Ubbelohde 510-1630 3 0.78 
Ostwald: Upper bulb 610-1620 4.21 0.665 
Lower bulb 460-1270 4.28 0.70; 
500 5.25 0.345 
Polymer PS 1 
Ubbelohde 760-1440 3.16 0.335 
Ostwald: Upper bulb 720-1780 3.13 0.333 
Lower bulb 570-1390 3.16 0.31; 
500 3.20 0.29; 
Polymer PMMA 2 
Ubbelohde 810-1550 2.66 0.32, 
Ostwald: Upper bulb 810-1620 2.71 0.28 
Lower bulb 640-1280 2.68 0.31: 
500 2.71 0.29, 





free chain ends. It must be remembered that for polymers in the present 
molecular weight range, the volume occupied by a single molecule is such that 
the concentrations employed for viscosity measurements are similar to the 
requirement for a system of close-packed spheres. From light scattering data 
(2), it has been estimated that the diameter of the sphere occupied by a 
polystyrene molecule of molecular weight 10° in benzene is approximately 
1000 A and this condition would be equivalent to a concentration of 0.24 
gm./100 cc. Thus, for most of the solutions used, the molecules are in close 
proximity providing a greater concentration of aggregates (15). A combination 
of factors, such as orientation and deformation of loosely connected coils, 
may then contribute to the shear dependence. 

For the Ubbelohde viscometer and the upper bulb of the Ostwald visco- 
meter, with a solution volume of 10.01 cc., the constant factor in the shear rate 
expression, i.e. rhgp/2ln, differs by only 0.5%. Thus, the slope constants 
determined at ‘free fall’ with these two systems should be practically identical. 
The fact that the values obtained for the graft copolymers with the Ostwald 
instrument are lower may be attributable to a greater degree of molecular 
disentanglement. The differences for the linear polymers are within the 
experimental accuracy. For the Ubbelohde determinations, dilutions were made 
in the viscometer and the solutions agitated intermittently for short periods 
(10-15 min.), whereas for the Ostwald measurements all solutions were shaken 
for 18 hr. before use. Prolonged agitation may be necessary to establish the 
aggregation equilibrium between single and double molecules for a given 
concentration. It was noted that the differences between reduced viscosities 
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for the two viscometers at similar concentrations and shear rates increased 
with decreasing concentration. 


In conclusion, it may be said that no detectable increase in the slope con- 


stant of the Huggins equation was found for a series of polymers with known 
amounts of branching when measurements were made at a constant shear rate. 


Al 


though this observation is contrary to theoretical predictions (14), it is 


recognized that the concentrations employed are such that polymer—polymer 
interactions will be of importance. Statistical treatments of the hydrodynamic 


be 


th 
in 


havior of long chain molecules are usually applied to the idealized condition 
at the particles are well separated in the solvent. The important difference 
behavior between branched and linear molecules appears to be their relative 


dependence of viscosity on the shear gradient. 
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THE VIBRATIONAL SPECTRA OF ISOBUTANE AND ISOBUTANE-d; ! 


By J. C. Evans? anp H. J. BERNSTEIN 


ABSTRACT 


The Raman spectra of liquid (CH3)s;CH and (CH3);CD were obtained photo- 
electrically and intensities and depolarization ratios measured. The infrared 
spectra of the corresponding gases were obtained from 3 to 35 uw. A vibrational 
assignment consistent with depolarization ratios and band contours has been 
made for the molecules on the basis of C3, symmetry and their spectra correlated. 


INTRODUCTION 


The vibrational spectrum of isobutane has been the subject of many 
investigations (7) and several ambiguities remain with regard to vibrational 
assignments (7). Sheppard and Simpson (7) suggest that further experimental 
work is required to decide upon the frequency assignment of the deformation 
modes and the asymmetrical skeletal stretching mode. In order to clarify the 
situation with regard to this molecule the Raman spectra of liquid HC(CHs); 
and DC(CH;); have been obtained with depolarization ratios, and the infrared 
spectra of the vapors have also been examined. The infrared spectrum of the 
heavy isobutane has been reported previously (2) but no attempt at assign- 
ment was made. 


EXPERIMENTAL 


Isobutane was obtained from a commercial cylinder and the deuterium- 
substituted compound prepared by Dr. L. C. Leitch. 

The Raman spectra of the liquids were obtained at —70°C. in a tube 
previously described (3) and depolarization ratios measured by comparing 
spectra obtained first with polaroid film wrapped with its axis parallel to the 
sample tube and then perpendicular to it. These ratios were corrected fcr the 
non-rectilinear irradiation set-up by the method of Rank and Kagarise (6). 
The mercury 4358 line was used to excite the Raman spectra, the Stokes lines 
of which are shown in Fig. 1. The photoelectric Raman spectrometer used (8) 
was operated under the following conditions: slits ~10 cm.—!, dynode voltage 
=500 v., with intermediate values of gain and response time. The current 
in the water-cooled mercury arcs was 17 amp. 

The infrared spectra were obtained with a Perkin Elmer Model 112 double 
pass infrared spectrometer with LiF, NaCl, and CsBr optics. The spectra 
shown in Figs. 2 and 3 were obtained in a 10 cm. cell at the various pressures 
and slit widths indicated. 


DISCUSSION 
Before attempting an assignment for this molecule it is worth while mustering 
all the information provided by theory. Symmetry considerations limit the 
discussion of the poiat group symmetry to either C; or C3,. Since steric hind- 
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rance probably rotates the CH; groups so that the planes of symmetry of the 
Cz, group are no longer present, one should actually consider isobutane to 
have the point group symmetry C; with the division of fundamentals into 
12A+12E with all modes active in both Raman effect and infrared absorption. 
For all practical purposes however we shall consider the molecule to have the 
higher symmetry of Cy to describe its vibrational spectrum, in accord with 
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Fic. 3. Infrared spectra of gaseous light and heavy isobutane in the region 1800-3100 cm.~!. 


previous investigators (7), for which the fundamentals fall into the grouping 
8A,+4A2+12E. In this case the 4A2 type fundamentals are forbidden in 
both Raman effect and infrared absorption. 

In order to understand the nature of the vibrational spectrum a little more 
clearly it simplifies matters to consider first the molecule in which the three 
H atoms of each CH; group are replaced by a single particle X. For this mole- 
cule then all bending vibrations involving X are external bending modes of 
the CH; group; the only internal mode for the CH; group in this molecule 
corresponds to the stretching vibration vcx. For the molecule HC(CX); of 
point group C3, there are 5A4:1+1A2+6E modes. It is readily seen that the 
distribution of the bending modes of the X group and, correspondingly, the 
CH; group is 14,:+1A2+2E. In Table I the modes of vibration of isobutane 
have been classified according to the approximate description of internal and 
external modes of the CH; group and also in terms of other groups in the 
molecules. Since the internal modes of the CH; group are well known (7) and 
give rise to bands at ~3000 cm.~', 1450 cm.—', and 1350 cm.—', the observed 
spectrum will be simplified to the extent that the 10 active internal modes of 
the CH; group will appear in the spectrum at the three regions indicated. 


The remaining bands of tie spectrum should belong to the three external 
C 


bending modes of the CH; group and the six modes of the H—-C—C group; 
Ne 

the E type torsional mode being rarely observed and the A: type not allowed. 

From considerations of this kind we readily see that the isobutane molecule 
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TABLE I 
Species No. CH; external CH; internal 
A, Vi vcH unsym. 
v2 ¥CH Sym. 
v3 VCH 
V4 dcu, unsym. 
V5 5cH, sym. 
v6 dccH,; 
v7 vcc 
vs dccc 
E V9 vcH unsym. 
V10 yvcu unsym. 
Vil vcH sym. 
vi2 dcu, unsym. 
Vi3 dcy, unsym. 
Via dcu, sym. 
Vis ocx 
Vi6 dccu, 
Vit dccH, 
Vi8 ¥CC 
Vig dcce 
20 Torsion 
Az v21 vcH unsym. 
22 dccH, 
V2 5cH, unsym. 
vo4 Torsion 





which has 36 degrees of vibrational freedom would be expected to give an 
observed spectrum of only 12 bands. Further in the Raman spectrum the A; 
type will be polarized and the E type depolarized and in the infrared spectrum 
of the gas the A, type will be parallel with a PR separation* of about 32 cm.~! 
and a pronounced Q branch (1), whereas the E type bands are perpendicular 
with less pronounced Q branch (1) and varying spacing due to Coriolis inter- 
action (4). 

Finally the Teller—Redlich product rule for isotopic substitution (4) may be 
used as a final check on the assignment since the theoretical product rule 
ratiost for the A, species is 1.403 and for the E type, 1.382. 

The observed Raman spectra and infrared spectra are tabulated and the 
bands assigned for isobutane and heavy isobutane in Tables II and III 
respectively, and in Table IV the spectra have been correlated. 

It is apparent that the depolarization ratios and infrared band contours 
unambiguously determine the species of the skeletal deformation modes. 
The question as to which frequency is to be associated with the degenerate 
skeletal stretching mode cannot be answered unambiguously since strong 
coupling in the deuterated molecule makes it impossible to use the isotopic 
shift as criterion. 

*The PR doublet separation was calculated from the formula of Gerhard and Dennison, Phys. 
Rev. 43: 197. 1933. The moments of inertia of isobutane were calculated assuming vcu = 1.09, 
vec = 1.54, with all angles tetrahedral and gave 

Ig = Ip = 107X10~" c.g.5. units, 
Ic = 183X10~-* c.g.s. units. 
For deuterium-substituted isobutane 


I, = Ip = 110X10~ €.g.5. units. 
{The moments of inertia required here were those given in footnote (*). 
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TABLE II 
RAMAN SPECTRA OF LIQUID HC(CH3)3 AND DC(CHs)s 

















Isobutane Heavy isobutane 
Ava Assignment 
Av Int.? p* Av Int.? p* 
2960 (12D) 2958 116 D V9, Vi0, Vir (e) 2956 116 D 
2936 (5p) 2933 94 0.27 v1, (a1) “{ 2938 74 0.57 
2927 80 0.43 
2908 (6p) 2907 103 0.18 ve (a1) 2903 108 0.25 
2891 (2) 2889 87 0.38 2vi2 or 2v13 (A4) 2882 42 0.25 
2869 (12p) 2869 167 0.13 v3 (a1) 2146 50 0.41 
(vg tvs A1) 2864f 87 0.23 
2778 (2p) 2783 38 60.58 2yi2 OF vi2-+13 (A) 2782 20 0.57 
2714 (5p) 2716 19 0.10 2vis (Ai) 2712 19 0.13 
~1468 16 ? vs (a) 1463 30 0.8 
1452 (7b)D 1450 35 D vie, vi3 (e) 1450 28 D 
V5 (a) 1380 6 D 
1322 ~1373 3 ? vis (e) 1358 2 . 
to (8b)D 
1355 1327 13 D vis (e) 944 18 D 
1184 8 ? v6 (a1) 1166 10 P 
1172 (4b)D 1169 15 D vi6 (e) 1229 25 D 
965 (4b)D 966 16 D vi7 (e) 1065 6 D 
795 (10)P 799 76 «600.24 vz (a1) . 796 70 0.17 
437 (2)P 433 3 0.35 vs (a1) 426 3 0.23 
371 (4b)D 367 8 D vig (e) 370 9 D 
1098 (3b?) Not observed v3 excited by 44047 
1326f 4 
? 1289 4 
? 891 2 
v, excited by \4047 1196 1 
v2 excited by \4047 1137. 1 
v3 excited by 44047 1096 1 
? 2168 ~10 0.53 
? 2127 ~12 0.49 





*Raman spectrum given in Kohlrausch, ‘‘Ramanspektren’’, Edwards Bros., Ann Arbor, Mich., 
194 5, p. 209. 

>Peak heights are given on the same scale for both molecules. 

‘These depolarization ratios have been corrected for convergence error. P = polarized, D = 
depolarized = 0.87. 

*Fermi resonance with 2v, (A). ; 

tvat+vs (A1) is a possible explanation for the intensity of 87 at 2864 cm. in the spectrum of 
DC(CH;)3: also incomplete deuteration. However since the largest change in intensity on deuterium 
substitution occurs here, 1.e. from 167 in the light compound to 87 in the heavy, this band must be 
identified with the CH stretching mode where deuterium substitution is involved. 

tDue to incomplete deuteration, t.e. CH(CH3)3 impurity. 


3000-1500 cm.— Region 

From the intensities in Table II it is readily seen that the band at 2869 
cm.—' in the light compound has become very much weaker in the deuterated 
compound compared to the intensities of the other bands in this region. This 
enables one to identify this frequency with the stretching of the tertiary CH 
bond. In the infrared spectra the situation is not so simple and indeed it 
would be very difficult to make this assignment on the basis of infrared data 
alone. The assignments of the gas frequencies in the infrared are complicated 
by the extensive overlapping and therefore cannot be considered final. In 
Table III however this attempt has been made in accordance with expected 
contours. All observed bands in this region have been assigned except the very 
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TABLE III 
INFRARED SPECTRA OF GASEOUS ISOBUTANE AND ISOBUTANE-d; 
HC(CHs)s DC(CHs)s 
Assignment 
v,cm.! Int. Type PR v,cm.! Int. Type PR 
separation separation 
2980 2980 \ 
-— w. (i) ™40 2 (A1) —_ w. (Pp ~™~40 
~~ 2940 sh ™ 2940 
Pn s. (L) vo, vio, vr (€) oo s. (L) 
2920 s. (i) mn (a1) 2920 s. @))) 
2904 2904 
m~2920 ~2920 
2904 s. (|) ™26 v2 (a1) 2904 s. (@)) ~26 
2894 2894 
— 2166 
2880 m. (|) vs (a1) 2152 s. t¢\)) 32 
2872 sh 2134 
™~ 2824 sh —_ 
2806 w. (I) Quiz or v2+v13 (A1) 2788 w. (I) 
2785 2768 
re w. (1) vss (E) pot w. (1) 
2645 — 
2629 w. (I) vate (A1) 2620 w.  ({)) 
m~ 2616 2604 
2570 v.w. ? 2584 v.w. 
? 2550 v.w. 
? 2535 V.w. 
po \ w. vet+vu (E) 
2480 sh w. ? 2464 v.w. 
? 2446 v.w. 
2422 
2402 \ v.w. (L) vst+unperturbed m7 (E)? 2410 V.w. 
2371 | 
2357 w. (i) 2vs (A1) 
~2310 | ~2310 | 
2286 v.w. (||) va tv7 (Ai) ~™~2298 v.w. (|) 
o~2265 
2170 v.W. vs+v7 (Ai) 
2125 v7+vu (EB) 
2046 
2Xunperturbed mz (Ai)? 2030 m 41) 28 
2018 
1967 v.w. ve+vr (A1) 
vetviz, vat+ns (E) 1885 v.w. 
1910 v.w. vmat+ys (Ai) 
1821 v.w. vwst+vs (E) ~1814 v.w. 
1760 v.w. vs+ne (E) 
1720 v.w. mate (A1) 1710 v.w. 
1698 v.w. vs+unperturbed vs (E) 
1644 v.w. mo +unperturbed ms (A1)? 1650 v.w. 
? 1620 
1493 ~1492.—s) 
1477 s. (Dd 32 ™ (a1) 1477 s. (I) 31 
1461 1461 
1406 1410 
1394 m (I) ~35 vs (a1) 1394 8. a1) ~37 
~™1371 ~™~1373 
a oe CH vu (e) on m. (1) 
a m. (1) ms (e) 7 s. (L) 
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TABLE III (concluded) 
INFRARED SPECTRA OF GASEOUS ISOBUTANE AND ISOBUTANE-d; 
































HC(CHs)s DC(CHs)s 
Assignment 
vy,cm.-1 Int. Type PR vy,cm.-! Int. Type PR 
separation separation 
1177 m. ve (a1) 
~1166 m. ne (e) a s. (1) 
ae } s. (1) nis (e) red m. (1) . 
? ~~ 907 i 
890 w. ? 
811 ~ 810 
797 s. (I) 30 v7 (a1) 793 s. (i) 31 : 
781 779 : 
~ 445 448 : 
426 m. (|) ~33 vs (a1) 432 m. (jl) ~35 
~ 412 ~™ 413 . 
TABLE IV 
VIBRATIONAL ASSIGNMENTS 
HC(CHs)3 DC(CHs)s Product rule ratios 
- 
Raman Infrared Raman Infrared Theor. Obs. | 
liquid gas liquid gas (Raman) ; 
A, v1 2933 2920 2927 2920 1.40 1.39 i 
v2 2907 2904 2903 2904 
V3 2869 2880 2146 2152 
% 1468 1477 1463 1477 
¥5 1373 1394 1380 1390 
V6 1184 1177 1166 
v7 799 797 796 793 
vs 433 426 426 432 
E v9 2958 2962 2956 2962 1.38 1.375 
V10 2958 2962 2956 2962 
vil 2958 2962 2956 2962 
viz 1450 ~1475 1450 ~1475 
v13 1450 ~1475 1450 ~1475 
vis ~1373 1371 1350 1367 
Vis 1327 1330 944 952 
vie 11692 1166 1229 1233 
v7 966° 1065 —_ 
vis 9172 918 813 812 
Vig 367 — 370 — 





*Assignment of K. S. Pitzer and J. E. Kilpatrick, Chem. Revs. 39: 435. 1946. 


weak features at 1620, 2446, 2464, 2535, 2550, and 2584 cm.~' in the infrared 
spectrum of DC(CHs)3, and the very weak band at 2570 in the infrared spec- 
trum of HC(CH;);3. In the Raman spectrum of DC(CHs); the bands at 2127 
and 2168 cm.— are difficult to assign on the basis of the observed fundamentals. 
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Region Below 1500 cm. 


In this region there is little difficulty in assigning the bands between 1480 
and 1300 cm.~! but comparison of the spectra of the light and heavy molecules 
between 1300 and 800 cm.~! shows very little resemblance to what one would 
expect for this particular deuterium substitution. In Fig. 1 for example the 
strong band at 1169 in the light compound seems to have been shifted to higher 
wave number on deuterium substitution whereas the band at 966 cm.~ al- 
though shifted in the expected direction to lower wave number has moved too 
much if the shift in v7 is a criterion. In Fig. 2 the anomalous isotopic shifts are 
even more pronounced in the bands at ~1180 cm. and 918 cm.!. There can 
be no question that there is strong interaction or coupling between these 
vibrational modes. 

One possible set of interactions to account for the observed spectrum is the 
following. If the bands at 1169 and 1320 cm.~ in the light compound are 
coupled, the unperturbed position of 1169 may be in the region to give the band 
at 1229 cm.~' in the deuterated compound. The unperturbed position of the CH 
deformation mode at 1327 cm.—' would fall perhaps to ~900 cm.~' upon 
deuterium substitution. This band would be strongly coupled with the un- 
perturbed band at 917 cm.~! (which should be relatively undisplaced in the 
deuterium compound), to give rise to the bands at 813 and 944 cm.—". Since all 
bands are strongly coupled in the light compound the unperturbed position of 
966 cm.—! may give rise to the band observed at 1065 cm. in the heavy com- 
pound. The band observed at 1166 cm.~ in the heavy compound is the polarized 
band arising from the polarized band at 1184 cm.— in the light compound; 
this band does not participate in the coupling because it is of different sym- 
metry species from the others. If this coupling system is sensible the ratios 
to be taken in calculating the product rule are 


O17 966 1169 1327 
813 “* 1065 “* 1229 944 











Rigorously the product of the ratios of the frequencies of the unperturbed band 
should be taken but since the interaction displaces bands in pairs by equal 
amounts in opposite directions and since these displacements are small com- 
pared to the band positions the ratios of the observed frequencies might be 
expected to be very nearly the same as those for the unperturbed values. 

Since the change in frequency of all modes other than v5 is expected to be 
negligible on deuterium substitution the observed product rule ratio will be 
given by the product of the ratios of the frequencies of the coupled bands. 
This ratio is calculated to be 1.37) which is consistent with the theoretical 
value of 1.38 given above for the E type modes. It is apparent that assignment 
of one of 917, 966, and 1169 cm.—' to the degenerate skeletal mode is not forth- 
coming from the present data; however if the suggested coupling scheme is 
indicative of the true state of affairs the mode at 1169 cm.— which has an 
unperturbed frequency in the twelve hundreds is not a likely candidate for this 
assignment. 
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The observed infrared spectrum in the region below 1500 cm.~! has been 
assigned except for the weak band at 890 cm.—! in CD(CH3)3 whereas in the 
Raman spectrum of DC(CH;); the corresponding band at 891 cm. and the 
bands at 1289 cm.—! remain unexplained. The assignment is essentially the same 
as that of Pitzer and Kilpatrick (5) differing in the following respects: 

(i) The A, wag of the CH; group vg is given here as the Raman band at 
1184 cm.—! in preference to the band at 1098 cm.—! which we did not observe 
in the Raman spectrum. The latter band is probably v3; excited by (4047. 

(ii) We have assigned the CH stretching mode frequencies which were all 
taken as 2950 cm. by Pitzer and Kilpatrick for the purpose of thermodynamic 
calculations. 

The present assignment will not alter materially the results of the thermo- 
dynamic calculations (5) since the uncertainty introduced in assigning the 
Az wag (v22) is certainly as large as the present change in vibrational 
frequencies. 
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ADSORPTION OF OXYGEN BY A SILVER CATALYST! 


By W. W. SMELTZER,? E. L. TOLLEFson,? 
AND A. CAMBRON4 


ABSTRACT 


An apparatus is described to measure volumetrically the amount of gas 
adsorbed by a catalyst at constant pressure and at exposure times as short as 
0.1 min. The volumes of oxygen adsorbed by a technical silver catalyst at 
pressures of 200, 400, and 700 mm. and in the temperature range of 180° to 
300° C. have been measured. Surface area determinations of 0.33 to 0.38 sq. 
meter per gm. of catalyst by physical adsorption and chemisorption methods 
show that each silver atom of the surface is associated with approximately one 
oxygen adatom at a monolayer coverage. The initial rate of adsorption of oxygen 
is directly proportional to its pressure. Observed transitions in the oxygen 
adsorption rates indicate that more than one type of oxygen complex occurs on 
the catalyst surface. The Elovich equation provides the best approximation of 
the adsorption rate data but is of limited applicability in interpretation of the 
mechanism of adsorption. Values of 22-29 kcal./mole and 17-25 kcal./mole have 
been calculated for the apparent activation energy and isosteric heat of adsorption 
respectively. 


INTRODUCTION 


Investigations by Twigg (13, 14) and Orzechowski and MacCormack (8) 
on the oxidation of ethylene to ethylene oxide and carbon dioxide over a silver 
catalyst have emphasized the role played by the chemisorption of oxygen in 
the reaction kinetics. There is, however, meager information about the oxygen 
adsorption characteristics of silver at elevated temperatures. As part of a 
program designed to elucidate some of the properties of a technical silver 
catalyst used in the oxidation of ethylene, the rates of adsorption and volumes 
of chemisorbed oxygen have been determined as a function of pressure and 
temperature in this investigation. The results reported are for the operating 
temperatures and oxygen partial pressures used in reaction studies with a 
silver catalyst prepared by the method of Cambron and McKim (5). 

Thus far, probably the most extensive investigation has been made by 
Benton and Drake (2) who examined the adsorption of oxygen in the tempera- 
ture range of — 183° C. to 200° C. by a silver catalyst prepared from precipi- 
tated silver oxide. Two types of adsorption were found, namely, a practically 
instantaneous physical adsorption at —183° C. which decreased rapidly with 
increasing temperature and an activated adsorption which proceeded slowly 
at 0° C. and increased rapidly at higher temperatures to give a saturation 
volume nearly independent of temperature over the range of 150° to 200° C. 
For this latter type of adsorption, the values of the isosteric heat of adsorption 
and energy of activation were 23-15 kcal./mole and 12.7 kcal./mole re- 
spectively. 

1Manuscript received April 4, 1956. 
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With the exception of this investigation, little appears in the literature on 
the adsorption of oxygen by silver at elevated temperatures. Taylor (11) 
measured the rate of oxygen uptake by silver in the temperature range of 0° 
to 184° C. and found that the energy of activation decreased with increasing 
adsorption from 17 to 13 kcal./mole. Twigg (13) found that the change in the 
electrical resistance AR of a silver plated glass wool catalyst increased in the 
presence of oxygen at a pressure p according to the relationship 


[1] AR = avV/p/(1+bvV >), 


where a and 6 are constants. Applying values for the solubility of oxygen in 
silver reported by Steacie and Johnson (10), Twigg was led to the conclusion 
that the resistance effect was caused primarily by oxygen adsorption. Attempts 
have been made to determine the number of oxygen atoms complexed with 
each surface silver atom. Buttner, Funk, and Udin (4) measured the surface 
tension of solid silver in helium-oxygen mixtures between 870° and 945° C. 
Using the data with the Gibbs adsorption isotherm, these authors concluded 
that 1.4 atoms of oxygen were associated with each silver atom at 932° C. 
Armbruster (1) measured the adsorption of oxygen by silver foil of preferred 
crystalline orientation in the temperature range of —195° to 20° C. The data 
were satisfactorily represented by the Langmuir isotherm. Using the calculated 
monolayer volume, this author concluded that one oxygen molecule was 
associated with four silver atoms. Since few data are available on the kinetics 
of oxygen chemisorption on silver above 150° C. and because the oxygen-silver 
atom ratio shows apparent variations with temperature, the adsorption of 
oxygen at exposure times as low as 0.1 min. has been measured in the tempera- 
ture range of 180° to 300° C. in this investigation. 


EXPERIMENTAL 
Catalyst Preparation 


Catalyst I was prepared from an alloy of silver containing 11.3% calcium 
and 0.11% nickel. A powdered sample in the mesh size range —35 +65 was 
activated by a steam exposure at 305° C. for seven hours to oxidize the calcium, 
a leach in boiling 20% acetic acid for four hours, and finally a thorough washing 
in distilled water. This activation procedure was repeated. Chemical analysis 
showed that the catalyst contained less than 0.05% calcium and 0.001% 
nickel. A charge of 49.0 gm. of this catalyst was placed in the adsorption cell. 

Catalyst II was prepared from an alloy of silver containing 10.0% calcium, 
0.14% nickel, and'0.14% copper. A powdered sample in the mesh range —20 
+35 was activated by the above procedure. Chemical analysis showed that 
the catalyst contained 0.08% nickel, 0.13% copper, and 0.08% calcium. A 
charge of 50.0 gm. of this catalyst was placed in the adsorption cell. 


Catalyst Reference State 


In order that the catalyst surface would be in the same reference state 
before each adsorption run, a standard evacuation temperature of 300° C. 
was chosen since a reasonable rate of desorption occurred without injury to 
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the catalyst by sintering. A standard pumping time of 20 hr. was adopted after 
it was established that the volume of oxygen adsorbed after 40 hr. pumping 
was the same, within experimental error, as that after 20 hr. for the same 
experimental conditions. 


Apparatus and Procedure 

The apparatus for the determination of the gas volumes adsorbed by the 
silver catalyst at various temperatures and oxygen pressures consisted essen- 
tially of an adsorption cell connected to a gas burette and a mercury mano- 
meter. With the exception of the pressure control device, the vacuum system 
was of conventional design and consisted of pressure gauges, an oxygen and 
argon purification train and storage, a mercury diffusion pump, and other 
standard accessories. 

The furnace, which was positioned around the adsorption cell by a slide, 
consisted of a 12 in. X 4 in. diameter aluminum cylinder of 1 in. wall thickness 
wound with a monel heating element and insulated with magnesia. The 
electrical current through the heater was regulated by a Bristol temperature 
controller which maintained the catalyst temperature within +1°C. The 
temperature of the adsorption cell was measured by two standardized chromel- 
alumel thermocouples fixed inside the furnace. 

Fig. 1 illustrates the adsorption section of the apparatus. The adsorption 
cell was connected to the gas burette by a short tube and a standard 3 mm. 
stopcock. To allow for end effects of the furnace on the apparent volume of 
the adsorption cell, the dead space calibrations were made at the experimental 
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temperatures with argon. The automatic pressure control device was included 
in this section of the apparatus to permit burette readings in rapid succession 
at constant pressure. For this purpose, manometer Me was equipped with 
electrical contacts P and Q and these were connected to the grid circuit of a 
tube driving a relay. The output of this single pole double throw relay was 
such that the action of the grid circuit could be reversed. If the output circuit 
was closed, the solenoid S was energized and raised the glass-enclosed iron slug; 
simultaneously the mercury level was lowered away from the sintered glass 
disk in the pressure control device. Thus, air flowed through the disk and raised 
the mercury in the gas burette when a pressure existed below the control 
pressure. 

Prior to the performance of an adsorption experiment, the standard de- 
sorption procedure was performed. The furnace temperature was then set to 
the desired value and stopcocks A, C, and D were closed. The mercury in 
burette B was lowered and the manometer M:2 was adjusted so that the column 
of mercury would contact point P when the system was brought to the run 
pressure. Stopcock E was then closed and oxygen was admitted through stop- 
cock ‘D from the storage bulb until the mercury in manometer M2 registered 
slightly less than the desired run pressure. At this stage, pressure was applied 
to the mercury surface in flask H by using the relay arrangement. When the 
stopcock G was opened, the mercury in the gas burette rose until it was 
automatically stopped when electrical contact was made between the lead P 
and the mercury column in manometer Mp2. After values of the temperature, 
pressure, and volume had been recorded, the mercury in the burette was 
raised by a volume equal to that of the dead space of the adsorption cell. An 
appropriate air pressure differential was then established across manometer M, 
and stopcock A was opened. The volumes of adsorbed oxygen at constant 
pressure were recorded by burette readings at 0.1 min. intervals or longer. 


Materials 


The silver—calcium alloys were prepared from high purity Royal Mint silver 
and high purity calcium supplied by Dominion Magnesium Limited. The 
minor alloying constituents were prepared from high purity electrolytic copper 
and Merck nickel shot of 98-99% purity. 

The argon and oxygen gases were obtained from commercial cylinders of 
the Dominion Oxygen Company with a specified purity of at least 99.7%. 
Prior to storage in the adsorption apparatus these gases were passed through a 
purification train consisting of ascarite to remove carbon dioxide, heated 
platinized asbestos in the case of oxygen to remove hydrogen, heated copper 
in the case of argon to remove oxygen, phosphorus pentoxide to remove water 
vapor, and a trap at liquid air temperature to remove residual condensible 
vapors. 


RESULTS 
The catalyst surface should maintain the same area and properties during 
the course of experimentation in order that an analysis can be made of the 


oxygen adsorption results at different temperatures and pressures. Fig. 2 shows 
the volumes of oxygen adsorbed by Catalyst I with a given set of experimental 
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of the catalyst at 300° C. for different time periods. 


Adsorption of oxygen by Catalyst I at 260° C. and 400 mm. pressure after evacuation 


conditions after the standard evacuation. The volume adsorbed at 400 mm. 
pressure and 260° C. after the catalyst was maintained at 300° C. for 100 days 
is approximately 81% of the initially adsorbed volume. Furthermore, the 
decrease in the adsorption capacity of the catalyst during the last 55 days is 
less than 5% of the volume adsorbed. The extensive results presented herein 


for Catalyst I were taken during this period. 


Illustrative plots of the volumes of oxygen adsorbed by Catalyst I versus 
time for each set of adsorption conditions studied in the temperature range of 
180° to 300° C. are given in Figs. 3 and 4. There is an initial rapid rate of 
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Fic. 4. Oxygen adsorption —time curves for Catalyst I at pressures of 200, 400, and 
700 mm. and temperatures of 220°, 260°, and 300° C. 


adsorption which increases with increasing temperature. This rate then de- 
creases to a low value and in the temperature range of 200° to 300° C. the 
volumes of adsorbed oxygen approach their equilibrium values after exposures 
of 25 min. For example, the adsorption in this temperature range after 15 min. 
was at least 98% of the total volume adsorbed after 25 min. The small increase 
in oxygen uptake during the last 10 min. may be due partially to the solubility 
of oxygen in silver. Isobars of the volumes adsorbed by the catalyst after 
25 min. exposure in oxygen and, also, these volumes corrected for the equi- 
librium solubility of oxygen in silver according to the values given by Steacie 
and Johnson (10) are illustrated in Fig. 5. Values are not given for 180° C. 
because adsorption was not complete in the experimental time. 

Results for the initial rapid adsorption of oxygen by Catalyst II after ex- 
posures of 0.1 min. are illustrated in Fig. 6. Although the rate of adsorption at 
temperatures above 220° C. was too fast to be measured accurately because 
of the initial pressure drop in the system, the curves indicate that this initial 
rate is proportional to the oxygen pressure. Reliance cannot be placed on a 
comparison of the slopes of the curves at different temperatures, as the adsorp- 
tion experiments were completed after the catalyst had been evacuated for 
only one week at 300° C. In a similar manner, the volumes of oxygen adsorbed 
by Catalyst I after 0.3 min. of exposure are, within the experimental precision, 
proportional to the oxygen pressure. 








1052 CANADIAN JOURNAL OF CHEMISTRY. VOL. 34, 1956 









OXYGEN ADSORBED -CC. 
> 
° 


30 — 


| 














200 220 


~ 
oS 
| | | ‘2 
240 260 280 300 
TEMPERATURE °C. 


Fic. 5. Isobars of oxygen adsorption by Catalyst I. 


@ A @ volumes adsorbed after 25 min. at 200, 400, and 700 mm. 


spectively with no solubility correction. 


O A O volumes adsorbed after 25 min. at 200, 400, and 700 mm. 


spectively with solubility correction. 





VOLUME ADSORBED- CC. 





, | 











° 200 


400 600 
PRESSURE - MM. 





pressure re- 


pressure re- 


Fic. 6. Volumes of oxygen adsorbed by Catalyst II after 0.1 min. and temperatures of 


220°, 260°, and 300° C. 





SMELTZER ET AL.: OXYGEN ADSORPTION , 1053 
DISCUSSION 
Surface Area and Silver-Oxygen Atom Ratio 


The isobars of Fig. 5 show that the volume of oxygen finally adsorbed is 
larger at lower temperatures at constant pressure. This indicates a process 
involving a state of equilibrium between the adsorbate and the gas phase. In 
an attempt to determine the catalyst surface area and some properties of the 
adsorbate, these results are analyzed according to the Langmuir adsorption 
isotherm. To substantiate this approach, results privately communicated by 
Drs. A. Orzechowski and K. E. MacCormack on oxygen and krypton adsorp- 
tion by a catalyst prepared in the same manner as used in this investigation 
are analyzed according to the Langmuir (7) and Brunauer, Emmett, and 
Teller (B.E.T.) (3) gas adsorption methods of obtaining surface areas re- 
spectively. 

One form of the Langmuir isotherm equation is 


[2] ‘ = Sate ’ 

E mOpP Vin 
where # is the pressure, V the volume of gas adsorbed, V,, the volume required 
to complete a monolayer, and 6 is a constant. Plots of 1/V vs. 1/p for the 
values recorded in Fig. 5 are given in Fig. 7. Linear curves satisfy the results 
for the temperatures of 200°, 220°, and 260° C. and the pressure range studied. 
This is valid for both the total uptakes of oxygen and these volumes corrected 
for solubility of oxygen in the catalyst. The monolayer volume calculated 
from the intercepts in the temperature range 200°-260° C. varies from 4.8 to 
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5.6 cc. and, if oxygen solubility is taken into account, from 4.3 to 4.9 cc. Thus, 
Vm values deviate approximately 15% over the temperature range investi- 
gated. If the determination at 220° C. is taken as a mean value, Vm is 5.3 cc. 
or 4.7 cc. depending upon whether an oxygen solubility correction ‘is made. 
These values (0.11 or 0.10 cc. per gm. of catalyst) correspond to areas of 0.41 
or 0.37 sq. meter per gm. if it is assumed that the cross-sectional area of an 
oxygen molecule is 14.1 sq. A. 

Orzechowski and MacCormack determined the adsorption isotherm with 
krypton as the adsorbate at 77.5° K. and, also, the volume of oxygen adsorbed 
at 220° C. by a 40 gm. sample of catalyst under the experimental conditions 
employed in this investigation. A surface area of 0.36 sq. meter per gm. of 
catalyst was obtained from these results when they were analyzed according 
to the B.E.T. method with the assumption that the cross-sectional area of the 
krypton molecule was 21.8 sq. A (9). A Langmuir plot of these results gave a 
monolayer volume of 4.0 or 3.5 cc. (0.10 and 0.09 cc. per gm. of catalyst) de- 
pending upon whether or not oxygen solubility was considered. These latter 
values, which correspond to surface areas of 0.38 and 0.33 sq. meter per gm., 
are in agreement with the values determined by the B.E.T. method and, also, 
with the values of 0.41 and 0.37 sq. meter per gm. determined in this investi- 
gation. 

This permits a simple approximation to be made as to the number of oxygen 
atoms complexed with each surface silver atom. ‘If the catalyst surface consists 
of cubic face-centered crystal faces, the effective area of a silver atom is 
8.32 sq. A. The B.E.T. determination of 0.36 sq. meter corresponds to 4.3 X 10" 
atoms of silver. This number is equivalent to either 0.04 or 0.08 cc. of oxygen 
depending upon whether one or two atoms of oxygen are adsorbed by two 
silver atoms. The monolayer volumes ranging from 0.09 to 0.10 and 0.10 to 
0.11 cc. of oxygen per gm. of catalyst determined from the results of Orze- 
chowski and MacCormack and this investigation, respectively, indicate that 
one oxygen atom is adsorbed by each silver atom under saturation conditions. 
Despite this agreement, the calculation is only a first approximation as only 
one type of crystalline face has been considered. 

The agreement of the values for the monolayer volume determined by a 
Langmuir and B.E.T. analysis does not substantiate the view that simple 
Langmuir considerations satisfy the characteristics of oxygen chemisorption 
on silver. This is so because, first, the calculated monolayer volume varies with 
temperature, second, the energy of adsorption is not constant, and, third, 
another analytical expression may be used to represent the experimental 
results. The Langmuir plots of Fig. 7, which do not yield a common intercept, 
illustrate the first limitation. An examination of the isosteric heat of adsorption, 
Q, may be used to illustrate the second limitation. Q may be evaluated from 
the equation 


_1(ainp 
” on3lan).. 


where R is the gas constant and p is the equilibrium pressure of a specified 
volume V of adsorbate at a temperature 7 on a surface of constant area A. 
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Although the accuracy is poor, the values calculated for the isosteric heat of 
adsorption (using the volumes of oxygen adsorbed shown in Fig. 5 without 
corrections for oxygen solubilities) vary from 25 to 17 kcal./mole for surface 
coverages of 68 to 82%. These values, which are in agreement with those 
reported by Benton and Drake (2), illustrate that the isosteric heat of adsorp- 
tion is not constant even at large surface coverages. 

The agreement of the experimental results, which may be fitted to the 
Langmuir isotherm equation, with the expression 


[4] V = av/p/(1+bv/), 


where V is the volume of adsorbed gas, p is the equilibrium pressure, and a 
and b are constants, illustrates the third limitation. This expression is similar 
to equation (1) found by Twigg (13) to satisfy the variation of the electrical 
resistance of silver in the presence of oxygen. To distinguish the applicability 
of these isotherm equations, measurements must be extended over a larger 
pressure range. We stress that the Langmuir isotherm equation has been 
applied arbitrarily over a limited pressure and temperature range to calculate 
the monolayer volume. 


Adsorption Kinetics 

In recent papers, Taylor and Thon (12) and Landsberg (6) have shown that 
the rates of chemisorption of several gases on solid surfaces may be represented 
by the Elovich equation, 
[5] dq/dt = ae~*%, 


where g is the amount of gas adsorbed, ¢ is the time, and a and a@ are constants. 
The integrated form of this equation is 


[6] q = (2.3/a) log(t+to) —(2.3/a) log to 


if tg = 1/aa. Thus, with a correctly chosen fo, the plot of g as a function of 
log(t+¢) should give a linear curve with slope 2.3/a. Also, a, which represents 
the initial rate of adsorption for g = 0, may be evaluated. If tp is too small to 
be assigned a value, a may be evaluated from 


[7] qe = (1/a) In(1+aat) 


if g: represents the amount of gas adsorbed after time ¢. Orzechowski and 
MacCormack (8) have fitted the data of Benton and Drake (2) for the chem- 
isorption of oxygen by silver to an equation of this type, one valid for the 
approximate range 0.3 < @ < 0.5 and another for the range 0.5 < @ < 0.7 
where @ is the fraction of the monolayer. In each range, a was directly pro- 
portional to the oxygen pressure. 

Although the adsorption versus time curves of Figs. 3 and 4 show that the 
Elovich logarithmic rate equation may not be expected to completely fit the 
adsorption results, it is advantageous to apply this equation in a discussion of 
the adsorption kinetics. Plots of the adsorbed volumes on Catalyst I versus 
the logarithm of time for the temperature range of 180° to 260° C. and an 
oxygen pressure of 200 mm. are shown in Fig. 8. There are, at the temperatures 
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Fic. 8. Adsorption data at 200 mm. pressure for Catalyst I plotted according to the 
Elovich equation. 


of 180°, 200°, and 220° C., abrupt changes in the slopes of the curves within 
the range of 2.2 to 3.0 cc. of adsorbed gas to regions of practically constant 
slopes. These latter slopes decrease at higher surface coverages in a manner 
which is more pronounced at higher temperatures. For the experimental con- 
ditions studied, the initial abrupt changes occur only at coverages corre- 
sponding to less than 3 cc. of adsorbed gas. It is assumed that a and a may be 
evaluated from the constant slopes and the volumes of gas adsorbed in these 
regions (equation [7]). These values and the surface coverages for which the 
curves are of approximately constant slope are recorded in Table I. Values are 
recorded only for the temperature range of 180° to 220° C. as the adsorption 
rate is too rapid at higher temperatures to determine the functional relation- 
ship with time. 

The calculated values for an initial rate of adsorption, a, increase with 
pressure and temperature. The dependence of the rate on pressure cannot be 
determined because the evaluations are not sufficiently accurate. However, 
the results presented in Fig. 6 indicate that this rate is directly proportional 
to the pressure. Moreover, it is to be noted that the calculated values for this 
initial rate are much larger than the experimentally determined values. In 
Fig. 9 plots are shown of the logarithms of the calculated and experimentally 
determined values for the initial rate of adsorption versus the reciprocal of the 
absolute temperature. As a measure of the experimentally determined rate, 
the reciprocals of the times required to adsorb 3 cc. of oxygen on Catalyst I 
at 200 mm. pressure of oxygen were used. The apparent energies of activation 
calculated from the slopes of the least squares lines are 27-29 kcal./mole and 
22 kcal./mole for the calculated and experimentally determined initial rates 
of adsorption respectively. This latter value of 22 kcal./mol2 is larger than the 
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values of 17 to 13 kcal./mole reported by Taylor (11) and Benton and 


Drake (2). 


The surface coverages recorded in Table I illustrate that the curves of 
constant slopes are confined to the approximate range 0.5 < @ < 0.9 if the 
experimental results are plotted according to the Elovich equation. In this 
analysis, fg was assumed to be equal to zero as positive values of this constant 
gave no better fit of the data to the equation. The above range corresponds to 


TABLE I 


VALUES OF CONSTANTS @ AND @ OF ELOVICH EQUATION AND SURFACE 
COVERAGE @ FOR WHICH EQUATION IS APPLICABLE 








6, 
fraction of 





Temperature, Pressure, a, a, 

a mm. 1/cc. cc./min. monolayer 

180 200 1.5 30 0.5-0.8 
400 1.4 “0 0.5-0.9 
700 13 55 0.6-0.9 

200 200 1.5 100 0.5-0.8 
400 1.4 150 0.5-0.8 
700 1.4 160 * -0.9 

220 200 1.7 320 0.6-0.7 
400 1.6 630 *-0.8 
700 1.5 670 * -0.9 





* Not observed. 


that of 0.56 < 6 < 0.7 found by Orzechowski and MacCormack (8) from an 
analysis of the data on oxygen adsorption on silver at lower temperatures. 
Also, the abrupt changes ia the slopes of the Elovich curves at a surface cover- 
age of approximately 0.5 are similar to a break in the curve for the adsorptioa 
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of hydrogen on zinc oxide found by Taylor and Thon (12). These authors 
suggest that abrupt changes in the slopes of the Elovich curves indicate a 
change from one kind of site to another at a certain stage of adsorption. 
This investigation supports the view that more than one kind of site is 
involved in the chemisorption of oxygen on silver. Figs. 3 and 4 show that the 
curves for adsorption at 180°, 200°, and 220° C. have transitions due to dif- 
ferences in the adsorption rates at 0.6, 0.2, and 0.2 min. respectively after the 
beginning of the adsorption runs. These abrupt changes in the adsorption-time 
curves and the curves for data plotted according to the Elovich equation occur 
at approximately half coverage of the catalyst surface. This suggests that 
initially the surface is rapidly covered by oxygen so that each oxygen atom is 
associated with two silver atoms. A relatively slow rearrangement of adatoms 
may then occur which permits further adsorption until one oxygen atom is 
associated with each silver atom. If this is correct, the different kinds of 
adsorption sites may be associated with differences in the bonding of oxygen 
with silver. For example, on the basis of surface geometry and bond angles, 


132A 
o—-0 
a aaa 
S/ 
w~ 
Ag Ag 288A Ag 





two species could exist on a catalyst surface at the interatomic cubic lattice 
spacings of 4.08 A and at the spacings for closest approach of silver atoms of 
2.88 A. Consequently, equilibria between the silver surface atoms and 
oxygen gas may be represented by the equations 


[8] 4A (eurtace) + O2(¢) = 2Ag20 (compiex)s 
[9] 2Ag20 (compiex) + Ox = 2Ag2O2(compiex)- 


The observed transitions in the adsorption rates then indicate approximately 
the termination of process [8] and the commencement of process [9]. In the 
temperature range of 180° to 220° C. process [8] is rapid whereas process [9] 
is relatively slow. Both of these processes are rapid at 260° and 300° C. as the 
transitions were not observed. This postulate of two surface oxygen complexes 
is only a crude approximation as no significance is given to the influence of 
variations in'the energy properties of different crystalline faces on the kinetics 
of adsorption, and cognizance is given to only the abrupt changes in the slopes 
of the Elovich curves. Although these latter abrupt changes occur at surface 
coverages of approximately 0.5, this does not disprove the proposal by 
Orzechowski and MacCormack (8) that different kinds of adsorption sites may 
be associated with different crystalline faces. 





SMELTZER ET AL.: OXYGEN ADSORPTION , 1059 
GENERAL DISCUSSION 


In an attempt to account for the observed phenomena, it has been postulated 
that the catalyst surface area may be determined by the Langmuir and B.E.T. 
gas adsorption methods and that at least two oxygen complexes are formed on 
the catalyst surface. These postulates lead to a qualitative interpretation of 
the adsorption kinetics and emphasize opposing views on the chemisorption 
of oxygen on a silver catalyst. Twigg (13, 14) and Orzechowski and Mac- 
Cormack (8), in their investigations on the oxidation of ethylene with a silver 
catalyst, have proposed that only single oxygen adatoms are involved in the 
adsorption and reaction kinetics and that proposed mechanisms are applicable 
to homogeneous and heterogeneous surfaces. Twigg assumes that the formation 
of ethylene oxide and carbon dioxide requires one and two favorably situated 
oxygen adatoms respectively. Selective poisoning is attributed to an alteration 
in the configuration of these adatoms. Orzechowski and MacCormack propose 
that both reactions involve single oxygen adatoms and that different energetic 
considerations account for the formation of ethylene oxide and carbon dioxide. 
According to their postulated reaction mechanism, selective poisoning is 
attributed to preferential poisoning of different crystalline planes or to the 
formation of impurity centers which influence the activation energies of the 
several reaction paths. On the other hand, transitions in the adsorption-time 
curves of this investigation at surface coverages of approximately 0.5 indicate 
that at least two oxygen complexes occur with one and two oxygen atoms 
equivalent to two silver atoms. If significance may be attached to the calcu- 
lated initial rate of adsorption (Table I) for the formation of: the latter com- 
plex, this rate is of a sufficient magnitude to influence the reaction kinetics of 
ethylene oxidation. Selective poisoning may then be attributed to the prefer- 
ential alteration in concentration of these different complexes. This concept of 
paired oxygen adatoms is different from that proposed by Twigg in that 
specific account is taken of the surface geometry and complete dissociation of 
adsorbed molecular oxygen need not occur. This study emphasizes that 
further experimental data are required for an understanding of the oxygen 
chemisorption phenomena on silver and that extreme caution must be used in 
representing the adsorption of oxygen by an analytical rate equation in reaction 
studies involving a silver catalyst. 
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THE REACTION OF DEUTERIUM ATOMS 
WITH ETHYLENE! 


By S. Tosy? AND H. I. ScHIFF 


ABSTRACT 


Deuterium was dissociated on a hot tungsten filament and the atom concen- 
tration measured by isothermal calorimetry. The recombination coefficient of 
deuterium atoms on a glass surface, coated with metaphosphoric acid, was found 
to be 3.81075, and similar to that found for hydrogen atoms. The reactions of 
H-atoms and D-atoms with ethylene were found to be very rapid. The effects 
on the yields of the products and on their isotopic composition of variations of 
reactant flow rate, atom concentration, pressure, and atom-detector position 
were studied. The major products were methanes, ethanes, and ethylenes, with 
minor amounts of propanes and butanes. The methanes were always highly 
deuterated while the ethanes were slightly deuterated. A mechanism is proposed 
to explain the observations based on a flow pattern in the reaction zone. The 
possibility of differences in the reaction rates of variously deuterated intermediates 
is also discussed. 


INTRODUCTION 


The homogeneous reaction between atomic hydrogen and ethylene occurs 
so rapidly that little is known of the detailed mechanism. Investigations using 
isotopic hydrogen have been largely qualitative (19, 11). This paper reports an 
attempt to obtain quantitative data on the products of the reaction between 
deuterium atoms and ethylene. 


EXPERIMENTAL 


The method used to produce and detect the deuterium atoms was similar 
to that used by LeRoy for hydrogen (2, 3, 7). About 2% of the hydrogen or 
deuterium in a fast flowing gas stream was dissociated into atoms on an elec- 
trically heated tungsten filament. The detector was a tightly wound platinum 
spiral, mounted on a rack which could be moved by a pinion attached to a 
standard tapered joint. Two teflon gaskets placed on the inner joint prevented 
the “‘Apiezon”’ grease from being squeezed out during rotation and kept the 
system vacuum-tight. This arrangement had the advantage of eliminating 
sliding contacts and permitted a constant lead resistance regardless of detector 
position. The wattage dissipated by the detector was determined from its 
resistance and the voltage drop across it. 

To conserve deuterium, the gas was circulated during an experiment by 
connecting the vent of the rotary pump through a refrigerated trap to a return 
line. Blank experiments showed that the exchange between the deuterium and 
the pump oil was less than 0.5% for six hours’ circulation. 

The reaction vessel was 22.5 cm. long and 25 mm. in diameter and was 
coated with metaphosphoric acid. The speed of gas circulation was 1400 
cm.’/sec., corresponding to a linear flow rate of 280 cm./sec. through the 
reaction chamber and a maximum reaction time of 0.08 sec. 

1Manuscript received A pril 3, 1956. 

Contribution from the Department of Chemistry, McGill University, Montreal; Quebec. 
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Deuterium was prepared from 99.7% deuterium oxide obtained from Atomic 
Energy of Canada Limited. The deuterium oxide was allowed to react with a 
carefully degassed mixture of zinc and freshly dehydrated calcium oxide. 
Deuterium was then evolved at 500°C. (9, 12) by the reaction: 


Ca(OD)2. + Zn — De + CaZnOs>. 


The yield of deuterium depended on the rate of heating of the furnace due to 


the competing reaction: 
Ca(OD)2 — CaO + D.O 


which occurs at temperatures above 350°C. Yields as high as 90% of the theo- 
retical value were obtained with rapid furnace heating. Mass spectral analysis 
showed the product to be 96% De, the remainder being almost entirely HD. 
Since the reaction under investigation also produced HD as one of the pro- 
ducts, additional quantities of D2 were periodically added to the system to 
prevent the deuterium content of the circulating gas from dropping below 96%. 

The ethylene, obtained from the Ohio Chemical Company, was ‘‘trap-to- 
trap’’ distilled several times. The center fraction was 99.9% pure, with ethane 
as the only detectable impurity. About four millimoles of ethylene-d, was 
kindly supplied by Dr. L. C. Leitch of the National Research Council, Ottawa. 

Deuterium atoms were allowed to react at room temperature with the 
ethylene in a conventional fast flow system. The effects of variation of reactant 
flow rate, atom concentration, deuterium pressure, and atom-detectcr position 
were investigated. 

The products were condensed in two large, liquid-nitrogen cooled traps, the 
second of which contained silica gel to retain the methanes. The total ‘“‘con- 
densible”’ fraction was first measured volumetrically. The sample was then 
divided in two and the mass spectrum of one portion determined. Unsaturates 
were removed from the other portion with acid mercuric sulphate solution 
(8) and the remaining saturated hydrocarbons analyzed mass spectrometri- 
cally. The analysis of the unsaturates was then obtained from the difference 
between the spectra of the treated and the untreated portions. 

Since standard samples of all the deuterated products were not available, 
approximate procedures were required to interpret the mass spectra. From a 
mixture known to contain only methane-d; and methane-d,, the probability 
factors a and 6 for the ease of removal of H- and D-atoms respectively from the 
parent molecule were measured for the mass spectrometer used. The values 
of these parameters for the other methanes were then estimated by linear 
extrapolation (20), and the mass spectra for all four deuteromethanes were 
calculated. 

The deuteroethanes were analyzed at reduced electron energies (25) using 
data kindly supplied by Dr. D. P. Stevenson of the Shell Development 
Company. Analyses were also performed at 70 volts electron energy, using 
published data (20) for ethane-d,, -d2, -d3, and -ds, and calculated spectra 
(24) for ethane-d, and -ds. Since highly deuterated ethanes were present in 
relatively small amounts, the errors involved in using the calculated spectra 
were not serious. The agreement between analyses at normal and reduced 
electron energies was satisfactory. 
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The deuteroethylenes were analyzed with the aid of published spectra (6). 
Small quantities of deuteropropanes and deuterobutanes were also found. 
Their total quantity was estimated from the sum of the ion currents in the 
“C,”’ and ‘‘C,”” mass number region. 


RESULTS 

The Recombination of Hydrogen and Deuterium Atoms 

Dingle and LeRoy found (7) that straight lines could be drawn through both 
first and second order plots for H-atom recombination, and so were unable to 
assess the relative importance of heterogeneous and homogeneous recombination 
reactions. We obtained similar results for both H- and D-atoms. Average values 
for the “‘pseudo”’ first order rate constant k’, and the “‘pseudo”’ second order 
rate constant k’’, for H- and D-atoms are given in Table I. The value for 


TABLE I 
PSEUDO RATE CONSTANTS FOR H-ATOM AND D-ATOM RECOMBINATION 








ky’ = 3.82+0.12 sec.-! ky’ 


kp’ = 2.77+0.08 sec.-! } *P 





- ky” = (2.52+0.06) X 10° cm.? sec.~! mole! ky” 
kp” = (1.66+0.04) X 10° cm.* sec.~! mole! 





ky’ can be compared with Dingle and LeRoy’s value (7) of 4.80 sec.—!. The 
ratio of either the first or the second order rate constants for hydrogen and 
deuterium should be less sensitive to wall conditions than the individual values, 
and clearly indicates that H-atoms recombine faster than D-atoms. There 
are no data with which these ratios can be compared, but Amdur (1) gives a 
value of 1.36 for kp’”’/ky’”’, the ratio of the third order rate constants at 30°C. 

The heterogeneous recombination of atoms can also be expressed in terms of 
a recombination coefficient y, which may be defined as the ratio of the number 
of atoms recombining on a surface to the total number striking that surface 
(21), and is, therefore, independent of the geometry of the apparatus. It can 
be calculated from the relation 


y = 2k’r/V 


where V is the mean atomic velocity and 7 the radius of the tube. Table II 
shows the values obtained in this and in other investigations. It will be seen 
that the values of y for H- and D-atoms are equal within experimental error. 
This suggests that any difference between the first order recombination rates 
is due only to the differences in atomic velocities. 


The Reaction of H- and D-Atoms with Ethylene 

Preliminary experiments with H-atoms showed that ethane and methane 
were the main products of the reaction, along with small amounts of propane 
and butane. The amounts of various products formed in the reaction with 
D-atoms are given in Table III. Fig. 1 shows the percentage of each product 
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TABLE II 


RECOMBINATION COEFFICIENT FOR GLASS SURFACE 
COATED WITH METAPHOSPHORIC ACID 





VOL. 34, 1956 





Author 











Atoms 7 X105 Temp., °C 
Poole (17) H 3.9 Room temp. 
Smith (22) H 2 15 
Dingle and LeRoy (7) H 4.5 5 
This work H 3.8 25 
This work D 3.9 25 
TABLE III 


VARIATION WITH REACTANT FLOW 














Partial Partial 
Total press. press. Percentage products 
Run Temp., press., atoms, ethylene, 
<<. mm.Hg mm.Hg mm.Hg Methanes Ethanes d-Ethylenes C3, C4 

D-21 24.0 1.36 0.0246 0.0182 97.9 2.1 0 0 
D-20 23.5 1.37 0.0191 0.025 98.7 1.3 0 0 
D-17 23.0 1.34 0.0222 0.0510 84.2 12.6 2.9 0.3 
D-27 25.8 1.36 0.0222 0.0599 73.0 — — — 
D-14 22.0 1.36 0.0199 0.132 70.5 18.5 8.7 2.2 
D-12 25.5 1.38 0.0211 0.235 59.2 22.5 14.1 4.2 
D-11 (25.5 1.41 0.0186 0.239 59.4 23.3 12.7 4.6 
D-24 25.0 1.37 0.0216 0.312 53.9 25.0 16.3 4.9 
D-18 23.5 1.36 0.0217 0.414 52.4 25.9 19.1 2.6 
D-26 25.5 1.36 0.0199 0.491 52.0 26.6 i 3.6 
D-25 26.5 1.36 0.0241 0.525 51.2 26.7 15.7 6.4 
D-15 24.0 1.34 0.0210 0.559 49.4 24.6 20.4 5.6 
D-16 23.5 lao 0.0221 0.625 51.6 24.3 19.6 4.4 





as a function of ethylene flow rate 


normalized for a D-atom flow rate of 1.0 


micromoles per second; the amount of each of the deuterated methanes, 
ethanes, and ethylenes is plotted against reactant flow rate in Figs. 2, 3, and 
4 respectively. 
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Fic. 4. Variation of ethylenes with ethylene flow rate. 
O, C:2H,/10; O, C:H;D; 4, C2H2D2; @, C2HD;; a, C2D,. 


The extent of deuteration of any product can be conveniently expressed in 
terms of the ‘‘deuterium number’”’, i.e. the average number of deuterium atoms 
per molecule of that product. The variation of deuterium number with ethylene 
flow rate is shown in Fig. 5. The deuterium numbers for the ethylenes do not 
include unexchanged reactant, and no values are shown at low flow rates, 
where the amount of unexchanged ethylene was too small for quantitative 
analysis. In one series of experiments the deuterium pressure was varied while 
the partial pressure of D-atoms was maintained at 0.028 mm. Hg by controlling 
the dissociator temperature; the partial pressure of ethylene was 0.017 mm. Hg. 
The results are plotted along with deuterium numbers in Fig. 6. The deuterium 
number was pressure dependent only in the case of ethane—with increasing 
pressure ethane-d,; became more abundant than ethane-d. 

Three experiments were performed in which the atom concentration was 
varied while the deuterium and ethylene pressures were kept constant. Table 
IV shows the results, and, for comparison, some published data for the re- 
action of ethylene with H-atoms produced by mercury photosensitization 
and by the discharge method. 

Since all atoms striking the detector are recombined, an attempt was made 
to alter the reaction time by placing the detector in the reaction chamber 
during an experiment. Fig. 7 shows the results of a number of experiments in 
which the detector position was varied from just below the exit tube of the 
reaction vessel (its normal position during other experiments) to just below 
the ethylene inlet jet. 
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Fic. 5. Variation of deuterium number with ethylene flow rate. 
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TABLE IV 
VARIATION WITH ATOM CONCENTRATION 














Percentage 








Ethane Butane 





Reference Atoms 
Moore and Taylor (16) H 
This work D 
This work D 
This work D 
H. I. Schiff H 


(unpublished work)? 


atoms 

Ratio— —— 
reactant Methane 

ca. 1073 2 

0.58 57.6 
1.32 fa.2 
2.28 76.0 
ca. 20 92.7 


14 84 
36.8 5.5 
24.4 3.3° 
22.5 1.5° 
7.3 0 





*Traces of acetylene were also found. 
>Includes some propane. 
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Fic. 7. Variation of products with detector position. 


These results indicate that the reaction zone did not extend more than 2 cm. 
from the inlet jet which corresponds to a reaction time of less than 0.007 sec. 


The reactiois: 


H + C.H, seas C2H; ’ 
D + CH, — C.H,D 


[1] 
[1a] 


must therefore be very fast and have collision yields greater than 10-4, in 
agreement with previously reported values (4, 5, 15) for reaction [1]. 
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DISCUSSION 


Perhaps the most unexpected result of this investigation was the production 
of C,H, and C2H;D as the’ most abundant ethanes in the reaction between 
D-atoms and ethylene. The amounts of these products were much too large to 
be accounted for by the hydrogen impurity in the deuterium. Other possible 
explanations for their occurrence were therefore analyzed. The following 
possibilities were considered: 

(1) The preferential dissociation of Hz on the dissociator filament due to its 
higher dissociation equilibrium constant relative to D2. However, calculation 
showed that this would increase the H-atom concentration by at most 0.1% of 
the total atom concentration. HD, the largest impurity in the deuterium, 
dissociates less readily at 2000°A. than does Dz. 

(2) The exchange of D-atoms with the phosphoric acid wall poison. If the 
entire increase in HD content observed when the deuterium was circulated 
for several hours in the absence of ethylene is attributed to this exchange, the 
maximum production of H-atoms would correspond to less than 1% of the 
total atom concentration. 

(3) Back diffusion of ethylene to the dissociator. When helium was substituted 
for the deuterium 100% recovery of the ethylene was obtained. 

(4) Loss of methane from the silica gel trap and its subsequent dissociation on 
the tungsten filament. In view of the excellent mass balances obtained this could 
not have occurred to any appreciable extent. Moreover, the dissociation of the 
methanes would not be a good source of H-atoms since they were found to be 
highly deuterated. 

(5) Heterogeneous reactions occurring on the atom detector. Experiments with 
the detector completely removed gave identical results. 

(6) Exchange of the ethanes with the mercuric sulphate solution during the 
removal of the olefins. The unsaturates were removed from several samples by 
reaction with purified bromine; no differences were found. 

(7) Errors in mass spectral analyses due to the lack of standard gases. Although 
these errors are difficult to assess quantitatively, comparison of analyses per- 
formed both at high and at low electron energies showed that the order of 
abundance of the deuteroethanes was correct. 

Since none of these possibilities is capable of explaining the production of 
light ethanes, it is reasonable to conclude that these were indeed products of 
the reaction between D-atoms and ethylene. This can be shown to be stoichio- 
metrically possible. The reactions forming the major products can be formally 
represented by the equations: 


4C.H, + aD — C,DaHi_2 + (a—2) H for methanes, 
C,H, + BD — C,DsH,4_s +8H for ethylenes, and 
C.H, + yD — C,DyHe_, + (y—2)H for ethanes, 


where a, 8, and y represent the deuterium numbers of the products. In a 
typical experiment a = 2.69, 6 = 1.50, and y = 1.35. The yields of the pro- 
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ducts were 0.377, 0.155, and 0.184 «M./sec. for the methanes, ethylenes, and 
ethanes respectively. Hence hydrogen atoms could have been generated by 
the reaction at a rate of (0.690.377) + (1.501.55) — (0.65.184) = 
0.37uM./sec. Considerably more H-atoms can therefore be produced in the 
deuteration of the products than are required for the formation of C2H. and 
C.H;D. This was found to be true in all cases. 

However, any mechanism for the reaction must be able to account for these 
light ethanes along with the other products. In the case of the reaction of H- 
atoms with ethylene the following elementary reactions may be considered: 


H + C:H, > C:H; [1] 
C:H; ~ C:H, +H [2] 

H + C:H; > CoH6* [3] 
C:H.* > C:H; +H [4] 
(M+) CoH6* — CoHe (+M) [5] 
C:H.* > CH, + H2 (6] 
C:H.* > 2CH; (7! 

H + CH; > CH,* [8] 
(M+) CH,* — CH, (+M) [9] 
2CH; — C2H.* [10] 


Corresponding reactions can occur with D-atoms. 


Formation of Ethylenes 
There are several possible explanations for the formation of the deutero- 
ethylenes. The primary step undoubtedly is 
D + C:H, > C2H.D. {1a] 


Since the deuteroethyl radical probably possesses excess energy it may readily 
dissociate again: 
C.H,D => C:H;D a H. [2a] 


The more highly deuterated ethylenes could be formed from repetition of these 
two reactions. This mechanism suggests that the yield of deuteroethylenes 
might be pressure dependent. This was not the case over the pressure range 
studied. 

Alternatively, the deuteroethyl radicals could undergo disproportionation 
(10). However, the very small butane yields indicated that disproportionation 
cannot be a major source of deuteroethylene. 

The third possibility is that reaction [1a] is followed by: 

D + C:H,D > C:H;D + HD. [6a] 


However, the repetition of [la] and [6a] does not provide an acceptable ex- 
planation for the formation of a// the deuteroethylenes since this sequence 
does not produce any H-atoms. In some of the experiments the amount of 
H-atoms which could be derived in the formation of the other products was 
insufficient to account for the amounts of C,H, and C2,H;D. This objection 
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could be overcome if the deuteroethyl radical first underwent exchange 
reactions of the type 
D + C.H,;D = C.H3;D.2 ~ H [4a] 


prior to forming ethylene by the counterpart of reaction [6a]. There is some 
evidence for [4a] from the study of the reaction of D-atoms produced by mercury 
photosensitization with ethane (23). Deuterobutanes were found which in- 
dicated exchange of ethyl radicals before recombination. 

Little is known about the likelihood of reaction [6] (or [6a]). It is exothermic, 
but would appear to be less likely than reaction [4], [5], or [7] by the ‘‘Prin- 
ciple of Least Motion’”’ (18). However, no unequivocal choice can be made 
between the first and third explanations for the formation of the deutero- 
ethylenes. 


Formation of Methanes and Ethanes 

For reasons already indicated, appreciable quantities of ethane could not be 
formed by disproportionation of ethyl radicals. Moreover, normal dispropor- 
tionation could not account for the formation of light ethanes. An unlikely 
‘double disproportionation”’ of the type 


C:H,D + C:H,D = CoH. + C.H2D, 


would be required. The remaining possible sources of ethane are reactions [5] 
and [10], both of which are preceded by reaction [3]. Methane formation, 
however, is also preceded by reaction [3]. The problem then arises of account- 
ing for the fact that the methanes were always found to be highly deuterated 
while the ethanes were only slightly deuterated. A possible explanation can 
be found by considering the two most likely fates for the methyl radicals in 
the system: 

H + CH; — CH, [11] 


The distribution of atoms in the reaction zone can be represented’ diagram- 
matically as: 

l l 

| Total atom concentration—high | 

| D/H ratio —high ¥ 

| wee eee eee eee ee eee ee Gas 

| stream 

| Total atom concentration—low 

| D/H ratio —low 





The total atom concentration decreases down the reaction tube as a result 
of both reaction and recombination. But the proportion of H-atoms in the gas 
stream will increase owing to the exchange of the D-atoms with ethylene and 
with free radicals. 

Reaction [11] will be favored over reaction [12] in the upper part of the re- 
action zone where the atom concentration is highest. But since this is where 
the D/H atom ratio is also highest, the methanes will be highly deuterated. 
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Under the conditions of low atom concentration at the lower part of the re- 
action zone reaction [12] will occur more frequently than reaction [11]. But 
the ratio of H- to D-atoms has increased so that the ethanes will be less 
deuterated than the methanes. 

This picture is in agreement with the observation (Table IV) that the yield 
of ethane decreases, and the yield of methane increases with increasing atom 
concentration. The reverse variation of yields with increasing ethylene con- 
centration is likewise consistent with this explanation. 

The large increase in methane formation with pressure. (Fig. 6) may be 
attributed to the effect of the deuterium molecule acting as the third body in 
reaction [9]. The third body requirement is probably more stringent in this 
reaction than in reaction [10]. The reason for the increase in deuterium 
number of the ethanes with pressure is not apparent. It conceivably could be 
due to excited ethyl! radicals formed in reaction [la] reacting with Dy». 

The experiments in which the atom detector was kept in the reaction cham- 
ber also suggested that the methanes are formed very near the ethylene inlet 
tube while the ethanes are formed lower in the reaction zone. Thus the meth- 
anes were unaffected by the presence of the detector while the deuteration of 
the ethanes was increased. 

Although this scheme is qualitatively in accord with the experimental 
observations, no account has yet been taken of the possibility that the reaction 
rates of atoms and radicals may depend on their isotopic composition. If the 
rates were independent of isotopic composition the deuterium distribution 
for any given product should be controlled by simple probability, and depend 
solely on the relative amounts of H- and D-atoms available for reaction. For 
example, the plot between the amount of each deuteroethane against the 
number of deuterium atoms in the molecule should give a simple probability 
curve, the maximum of which should correspond to the deuterium number. 
The heterogeneous reaction between ethylene and molecular deuterium over a 
heated nickel catalyst (14, 25) produced deuteroethanes which fitted a simple 
probability curve with a maximum at C.H,Ds». Keii, however, has pointed out 
(13) that this ag.eement may have been largely fortuitous. In the present 
investigation the less deuterated products were more abundant than would be 
expected on probability grounds. 

Further evidence for the dependence of reaction rate on isotopic constitution 
was obtained from experiments in which H-atoms were reacted with ethylene- 
d,. In the reaction between D-atoms and C.2H,, the most abundant methane 
produced was CD3;H. In the reaction between H-atoms and C,D, the most 
abundant methane was CH2D2, and not CDH3. Similarly, although CH, 
was the ethane formed in greatest amount in most of the D+C2H, experi- 
ments, C.D, was not the most abundant ethane formed under similar conditions 
in the reaction of H with C,D,. Owing to the large number of steps required to 
form some of the products, even a small difference in reaction rates of isotopi- 
cally different species could account for the differences in the isotopic distribu- 
tions of the products. 
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Because of the complexity of the reaction no decision could be reached as to 
what extent the deuterium distribution was due to the earlier explanation 


5 based on the flow pattern in the reaction zone, and to what extent it was due 
, to differences in reaction rates between variously deuterated intermediates. 
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THE REACTIONS OF ACTIVE NITROGEN 
WITH ALKYL CHLORIDES'! 


By B. Dunrorp,? H. G. V. Evans,? a4np C. A. WINKLER 


ABSTRACT 

The reactions of active nitrogen with methyl, ethyl, vinyl, propyl, and iso- 
propyl chlorides yielded hydrogen cyanide and hydrogen chloride as the main 
products. Small amounts of cyanogen and a polymer were formed from all the 
halides, and all except methy! chloride also yielded small amounts of C2 and C; 
hydrocarbons. The observed changes in amounts of products recovered with 
different reactant flow rates were characteristic of a fast reaction in which com- 
plete consumption of either reactant occurs when the other is present in excess. 
Mechanisms for the reactions are suggested on the basis that relatively long- 
lived complexes are formed in the initial attack of a nitrogen atom on the alkyl 
chloride. 


INTRODUCTION 

A number of reactions of active nitrogen, produced in a condensed electrical 
discharge, have been studied in this laboratory. Emphasis has been placed 
largely, though not exclusively, on the reactions with hydrocarbons (2, 3, 9, 
11, 13, 14, 17, 18, 19). Hydrogen cyanide is the main nitrogen-containing pro- 
duct of these reactions, which implies that the reactions with saturated hydro- 
carbons involve direct attack at a completely shielded carbon atom. They are 
much faster than would be expected from the principle of least motion enun- 
ciated by Rice and Teller (15). 

In an attempt to elucidate further the mechanism of the reactions, it 
appeared desirable to investigate the effect of a chlorine substituent on the 
reactivity. Accordingly, the reactions of active nitrogen with methyl, ethyl, 
vinyl, propyl, and isopropyl chlorides have been studied. 


EXPERIMENTAL 

The reactions were studied with conventional fast-flow techniques (11). 
Active nitrogen was produced by a condensed discharge and mixed with the 
halides in a reaction vessel that was sufficiently far from the discharge tube 
to prevent back diffusion into the discharge. The flow and activation of the 
nitrogen were kept constant in experiments with a given halide and the effects 
of change in flow rate of the halide determined. The reaction vessel could be 
surrounded by a furnace to study the effect of temperature. Condensable 
products were trapped with the usual arrangement of liquid air traps, and were 
analyzed by low temperature fractionation (12) and with a mass spectrometer.’ 
The hydrogen cyanide, cyanogen, and hydrogen chloride produced in the 
reactions were estimated by titration (16, 20). The hydrogen produced in the 
reaction with methyl chloride was determined mass spectrometrically on 
samples of the gas issuing from the vacuum pump. 

1 Manuscript received March 28, 1956. 

Contribution from the Physical Chemistry Laboratory, McGill University, Montreal, Quebec, 
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analyses. : 
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RESULTS 

All the reactions produced a peach-colored flame of fairly high luminosity, 
which was confined to the upper part of the reaction vessel at the higher flow 
rates. 

The main products of all the reactions were hydrogen cyanide and hydrogen 
chloride; small amounts of cyanogen were also formed. In all the reactions, 
significant amounts of a brown polymer were formed on the walls of the 
reaction vessel but not in the traps. This polymer was quite inert, and had 
properties somewhat resembling those of paracyanogen. Its composition varied 
appreciably with the different alkyl chlorides and to a lesser extent with 
different flow rates of a given alkyl chloride. It contained nearly equal amounts 
of carbon and nitrogen (on a mole basis) together with smaller amounts of 
chlorine and hydrogen. 

The formation of polymer on the walls of the reaction vessel made it difficult 
to obtain reproducible results at higher flow rates. (The concentration of active 
nitrogen is very sensitive to wall conditions, and the polymer itself was 
probably formed by a wall reaction.) 

Methyl Chloride 

The rates of production of hydrogen cyanide, hydrogen chloride, and 
hydrogen, and. the rate of consumption of methyl chloride, all varied linearly 
with methyl chloride flow rate up to a flow rate of about 3X10-* mole/sec., 
and were independent of temperature in this linear region (Fig. 1). In this 
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Fic. 1. Relations between methy! chloride flow rate and rates of methy! chloride consump- 
tion and product formation. 


region, only 89% of the methyl chloride was consumed, independently of flow 
rate or temperature. At higher flow rates the yields of products approached a 
limiting value which was the greater, the higher the temperature. 

The curves in Fig. 1 have only qualitative significance, since deposition of 
polymer changed the surface condition of the thermocouple wells at the’center 
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and top of the reaction vessel and the thermocouple readings were therefore 
only a rough measure of the reaction temperature. Quite marked changes in 
the readings were observed both during an experiment and in successive 
experiments. In addition, the reaction was sufficiently exothermic to create 
marked temperature gradients in the reaction vessel. The temperatures at the 
center of the reaction vessel have therefore been recorded, rather arbitrarily, 
as ‘low’ and “high” corresponding roughly to 100° and 300°C. respectively. 

The amount of cyanogen formed was quite erratic, but it was always less 
than 3% of the methy! chloride consumed. A careful search failed to reveal the 
presence of any cyanogen chloride, methane, or ethane in the products. 
Ethyl, Vinyl, and Propyl Chlorides 

In these reactions, small amounts of hydrocarbons were formed as secondary 
products. Vinyl chloride gave only acetylene, ethyl chloride gave a mixture of 
C, hydrocarbons, and the propyl chlorides gave mixtures of C2 and C; hydro- 
carbons. The composition of these products did not change significantly with 
flow rate in the limited range where samples sufficiently large for analysis 
could be recovered. Typical values are: for ethyl chloride, 94% ethylene, 
5% acetylene, 1% ethane; for n-propyl chloride, 34% propylene, 41% ethyl- 
ene, 17% acetylene, 6% propane, 2% ethane; for isopropyl chloride, 42% 
propylene, 33% ethylene, 18% acetylene, 6% propane, 1% ethane. 

Figs. 2 and 3 show the changes in amounts of products recovered with dif- 
ferent reactant flow rates. The different limiting yields with the various 
reactants may be attributed to changes in operating conditions, particularly 
the wall conditions in the reaction vessel. The shape of the curves for hydrogen 
cyanide production from the alkyl halides is quite similar to that of the 
analogous curves for the olefins (9, 17, 18), and is typical of a fast reaction 
in which there is complete consumption of either reactant when the other is 
present in excess. Although the break in the curve for hydrogen cyanide 
production is not as sharp for vinyl chloride as for the other chlorides, the 
general character of the reaction and comparison with the similar behavior 
of acetylene (19) leave no doubt that the reaction of active nitrogen with 
vinyl chloride is fast. 

The reactions of the higher halides were even more exothermic than the 
reaction with methyl chloride, and the difficulties of temperature measure- 
ment were aggravated to such an extent that it was not possible to obtain 
meaningful results. A few experiments with ethyl chloride in a heated reaction 
vessel indicated that the yields of products were independent of temperature 
within experimental error. 


DISCUSSION 


There is some evidence that active nitrogen contains more than one chemi- 
cally reactive species (7, 8, 14). However, considerations outlined in a further 
paper (4) suggest that one species, probably atomic nitrogen, is responsible 
for most of the reaction observed with the hydrocarbons and related com- 
pounds, and it is on this basis that the present results will be discussed. 
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The observed initial linear variation in yield of hydrogen cyanide with 
methyl chloride flow rate, and its temperature independence, appear to be 
possible only for a fast reaction. It can be estimated from the experimental 
conditions that this behavior must correspond to a collision yield greater than 
about 5X10-°. Similar results would probably be obtained for the higher 
halides, but the low flow rate region was not investigated as thoroughly for 
these reactions. 

The recovery of methyl chloride at low flow rates might be attributed mainly 
to decomposition of nitrogen atom-— methyl chloride complexes back to 
methyl chloride and nitrogen on the cold surface of the trap, although some 
channelling of methyl chloride in the flow through the reaction vessel might 
have occurred. Competition between a fast reaction of nitrogen atoms with 
methyl chloride and fast recombination of nitrogen atoms catalyzed by methyl 
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Fic. 3. Relations between rates of product formation and flow rates of propyl and isopropyl 
chlorides. 


chloride or one of the reaction products can explain the observation that the 
vields were independent of temperature at low flow rates and increased with 
temperature at high flow rates. Similar, but more pronounced effects have been 
found in the reaction of methyl cyanide with active nitrogen where the tem- 
perature measurements were more trustworthy (5). A mathematical treat- 
ment of the problem is in preparation. 

The features of the methyl chloride reaction outlined above suggest that a 
relatively long-lived methyl chloride — nitrogen atom complex was formed. 
The assumption that such complexes are formed in the reactions of the higher 
halides appears to be necessary to explain the formation of the secondary 
products. Further, the close similarity of the reactions of the higher halides 
to those of the unsaturated hydrocarbons suggests the possibility that hydrogen 
chloride is eliminated during formation of the complex, to yield a nitrogen 
atom complex with the corresponding unsaturated hydrocarbon. This is more 
fully discussed in the subsequent paper (4) on the basis of which it is possible 
to suggest the following mechanisms for the formation of the products obtained 
from the various active nitrogen — alkyl chloride reactions: 
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Methyl Chloride 
CH;Cl + N — N.CH;Cl —- HCN + HCl + H 
H +H -—H, 


Ethyl Chloride 
C:H;Cl + N > N.C.H,Cl ~ HCl + N.C:H, 
N.C:H, ~ HCN + CH; 
CH; + N—HCN + 2H 
N + N.CsH, > Ne + CoH, 
N + N.CsH,y > N2 + CoH4* (possibly CH2CH>2*) 
C2H,.* — He + C2H.* (possibly CH=CH*) 


Propyl Chloride 
C3;H;Cl + N — N.C;H;Cl — HCl + N.C;3He 
N.C;3Hs ~ HCN + C,H, +H 
N.C:He > Nz + CsHe 
N + N.C;Hs > N2 + C3He* (possibly CH;—CH—CH:"*) 
C3H«* — CH;—CH*: + CH 
CH;—CH?*: > CoH, 
CH;—CH*: — C:H:2 + H2 


Vinyl Chloride 
CH.CHCI + N — N.CH:CHCI — HCI + N.C:H2 
N.C:H2 + N > N2 + C:He2 
N + N.C2H2 — 2HCN 


It will be noted that the formation of the secondary hydrocarbon products 
is accompanied by recombination of two nitrogen atoms. This recombination 
appears to be the only source of sufficient energy to form these products and, 
for their formation in significant amounts, it is obviously necessary that a 
relatively long-lived complex should be formed. Thus, the higher halides 
might be expected to show the same change in yields with temperature at 
higher flow rates that was observed for methyl chloride. In spite of experi- 
mental difficulties, there was some indication that the effect was less marked 
for the higher halides. This might indicate a lower activation energy for these 
reactions. On the other hand, the spontaneous decomposition of complexes 
containing more than one carbon atom would make available carbon-contain- 
ing intermediates that should react very rapidly with nitrogen atoms (1) and 
the direct reaction should be favored over the recombination. 

The cyanogen produced in the alkyl chloride reactions was probably formed 
by the combination of two cyanogen radicals. The recombination of two nitro- 
gen atoms makes available some 225 kcal. of energy, and there is a quite high 
probability that a substantial part of this energy is transferred from the 
newly forming Nz molecule to the internal degrees of freedom of the rest of the 
complex (4). This can bring about dissociation into very reactive fragments. 
For methyl chloride, for example, reactions of the following type are all 
exothermic: 
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N + N.CH;Cl — Ne + CH.Cl + H 
N + N.CH;Cl — Nz + HCl + CH + H2 
N + N.CH;Cl — N2 + CH, +H + Cl. 


It seems rather likely that cyanogen radicals will be formed by reaction of such 
fragments with still another nitrogen atom, e.g. 

N + CH: —~CN + He 

N+CH —CN +H. 


These reactions are sufficiently exothermic to give the electronically excited 
cyanogen radicals that are responsible for emission of the reaction flames (4). 
For the higher halides it is possible that some cyanogen was formed directly 
by such reactions as 

N + N.C2H, — (CN)2 + 2H2. 


However, this type of reaction involves a great deal of rearrangement, and 
more likely are reactions such as 
N + CH: ~ CN + He 
2CN — (CN)>. 


The fact that cyanogen production passes through a maximum at low chloride 
flow rates seems to be good evidence that more than two atoms of nitrogen are 
involved in its formation. If only two nitrogen atoms were involved, cyanogen 
production from the ethyl chloride reaction should be similar to the ethylene 
production, at least at high flow rates where the amount of ethylene reacting 
with nitrogen atoms can be neglected. 

According to the mechanism suggested above, the production of unsaturated 
hydrocarbons should pass through a maximum as the chloride flow rate is 
increased. The chloride flow rate at which this maximum should occur will 
depend upon the relative rates of the reactions of the chlorides and the un- 
saturated hydrocarbons with nitrogen atoms. Apparently, the maximum yields 
were not attained in the region studied and this, together with the observation 
that little or no unsaturated hydrocarbons were recovered at low flow rates, 
gives some indication that nitrogen atoms react faster with the unsaturated 
hydrocarbons than with the alkyl chlorides. 

The presence of highly reactive intermediate species in the reactions of the 
halides with nitrogen atoms, and particularly the presence of free CN radicals, 
probably accounts for the formation of polymer by a wall reaction involving 
mainly CN radicals with smaller amounts of the other intermediates. This 
would mean that recombination of two nitrogen atoms was necessary for poly- 
mer formation, in agreement with the observation that the polymer was 
formed near the top of the reaction vessel, especially at higher flow rates. It 
would also account for the differences in composition of the polymer from the 
different reactions. 

The mechanisms outlined above involve production of hydrogen atoms, and 
the small amounts of saturated hydrocarbons formed in the reactions of the 
higher halides probably resulted mainly from hydrogenation of the correspond- 
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ing olefin. Failure to detect methane as a product indicates that the con- 
centration of atomic hydrogen was low in all the reactions. Hydrogen chloride 
(10) and methyl! cyanide (6) are efficient recombination catalysts for hydrogen 
atoms. It seems quite possible that alkyl halides would behave like methyl 
cyanide in this respect, and in any event there was probably enough hydrogen 
chloride formed in the reactions to keep the concentration of hydrogen atoms 
at a low value. 

The occurrence of secondary reactions involving more than one nitrogen 
atom makes it necessary to be cautious in drawing conclusions about the rate 
of the primary reaction, particularly since many of the products of these 
secondary reactions are known to react fast with nitrogen atoms. However, it 
seems likely (4) that the main energy barrier for the reactions is surmounted 
when the complex is formed. This conclusion would have to be modified if the 
* secondary reactions formed products capable of reacting with alkyl chlorides 
to set up a chain reaction. For the methyl chloride reaction, hydrogen and 
chlorine atoms are possible chain carriers. However, the absence of any detect- 
able amounts of methane, ethane, chlorine, cyanogen chloride, or ethylene 
chloride in the products makes it seem unlikely that there was an important 
contribution from chain reactions carried by hydrogen or chlorine atoms. It is 
possible to write a reaction scheme for a chain carried by CN radicals: 


N + CH;Cl — N.CH;Cl 
N + N.CH;Cl — Nz + CH:Cl + H 
N + CH:.Cl — CN + HCI +H 
CN + CH;Cl ~ HCN + CH.CI. 


The third reaction is exothermic and should be very fast (4), but it is doubtful 
that the last reaction would be fast enough to be effective as a chain carrier. 
Absence of a significant contribution from a chain reaction is also indicated by 
the reactions of the higher halides. In these, the recovery of unsaturated 
hydrocarbons with the same number of carbon atoms as the parent halide 
provides good evidence that the main path of the reaction was not hydrogen or 
chlorine atom abstraction followed by nitrogen atom attack on the resulting 
radical. 

It appears probable, from the foregoing discussion, that the reactions of 
nitrogen atoms with alkyl chlorides involve direct attack at a shielded carbon 
atom without requiring a high activation energy. 
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THE VIBRATIONAL SPECTRA OF ACETALDEHYDE 
AND ACETALDEHYDE-d; ! 


By J. C. Evans? AND H. J. BERNSTEIN 


ABSTRACT 


The Raman spectra of liquid CH;CHO and CH;CDO were obtained photo- 
electrically and depolarization ratios measured. Infrared spectra of the vapors 
and solids were also obtained. A vibrational assignment is given and the spectra 
correlated. The vibrational spectra of the acety] halides reported in the literature 
are also correlated. 


Several infrared and Raman investigations of acetaldehyde and acetalde- 
hyde-d, have been made but complete agreement as to the vibrational assign- 
ments has not been reached (3, 4, 6, 7, 8, 9, 12, 14). Pitzer and Weltner’s re- 
assignment (9) is consistent with the thermodynamic data but has not been 
definitely established. The reported depolarization ratios of some of the Raman 
bands (4) are inconsistent with this assignment and, recently, a different 
assignment for the low lying bands has been proposed (6). 

The most promising source of further relevant information seemed to be the 
vibrational spectrum of acetaldehyde-d, (CH;CDO) and an investigation of 
the Raman and infrared spectra of this molecule has been made. The Raman 
spectrum in the liquid phase and the infrared spectra in the vapor and solid 
phases have been examined. Similar studies were made on acetaldehyde. 

The vibrational assignment for acetaldehyde which is indicated by these 
results is very similar to that made by Pitzer and Weltner; the revisions are 
mainly in the higher wave number region and have negligible effect upon the 
calculated thermodynamic functions. 

The spectra of the related molecules, the acetyl halides, may now be corre- 
lated with some success. 


EXPERIMENTAL 


Acetaldehyde (Eastman Kodak) was refractionated. Acetaldehyde-d; was 
prepared and fractionally distilled by Dr. L. C. Leitch. 

A White grating spectrometer with photoelectric recording (16) was used to 
obtain the Raman spectra of the two compounds in the liquid phase. De- 
polarization ratios were measured (2) and corrected for convergence error (10). 
The infrared spectra were obtained with a Perkin Elmer Model 12C (double 
pass) spectrometer with LiF, CaF,, NaCl, and KBr optics. A low temperature 
cell of the type described by Wagner and Hornig (15) was used to measure 
the infrared spectra of the solids. 

The spectra are illustrated in Figs. 1 to 4 and the data are collected in 
Tables I and IT. 


1Manuscript received April 16, 1956. 

Contribution from the Division of Pure Chemistry, National Research Council, Ottawa, Canada. 
Issued as N.R.C. No. 3998. 

2National Research Council Postdoctorate Research Fellow 1953-55. 
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2762 
2749 


1422 


1390 


1349 


1120 
1103 


888 sh 


883 


TABLE I 





CH;CHO 
Raman 
liquid Poorr 
3001 0.87 
2964 ~0.76 
2917 0.10 
2843 0.22 
2751 0.31 
2731 0.36 
~1740 sh ? 

1714 0.33 
1426 ~0.8 

1391 ~0.65 
1342 ~0.8 

1109 0.34 

911 0.45 

885 p. or dp. (?) 

767 0.87 

512 0.40 


+} 21352 = 








21743 = 3486 (A’) 
1743 +1400 = 3143 (A’) 
a’; v(CHs3) 

a’ ; v(CHs) 


a : v(CHs3) 
a’ ; v(CH) 


2800 (.4’) 


2X 1400 A’) 
2704 (A’) 


and other combinations of 
1352, 1400, and 1430 cm."! 
1743 +(830) = 2573 (A’’) 
1352+1111 = 2463 (A’) 
919+1352 = 2271 (A’) 


9999 


eons 


2X1111 = (.4’) 


919+1111 = 2030 (A’) 


2X919 = 1838 (A’) 
919 + (830) 


a’; v(C=0) 


(CH3) 


a’ and a”; def. 
a’; C—H rocking 


a’; def. (CHs) 
a’; v(C—C) 
888+tors. (~230 cm.~) 


a’; CH; rocking 


a’’; CH: rocking 


A”; 509+tors. (?) 
a’; C—H wagging 


a’; sC—C=O 





In Tables I and II estimated positions are used for the coupled modes in determining the wave 
numbers of overtones and combination tones. These estimated values are bracketed. 
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TABLE II 
CH;CDO 
Infrared 
Raman Assignment 
Vapor Solid liquid Pear. 
3490 
3405 2X1743 = 3486 (A’) 
3462 
3173 1743+1431 = 3174 (A’ and/or A”) 
3096 1743+1353 = 3096 (A’) 
3028 3008 2998 ~0.86 a’; v(CHs) 
3000 
2970 2968 2965 ~0.76 a’ ; »(CHs;) 
2946 
2921 2917 0.13 a’ ; v(CH3;) 
2840 2859 21431 = 2862 (A’ and/or A”) 
2733 2764 1353 +1431 = 2784 (A’ and/or A”’) 
2704 2X 1353 = 2706 (A’) 
2600 1743+(890) = 2633 (A’) 
2528 1431+1109 = 2540 (4’, A”) 
2476 
2465 1353+1109 = 2462 (A’) 
2450 
2416 . 
1743+668 = 2411 (A”’) 
2400 
2222 2234 1431+802 = 2233 (A’ and/or A”) 
2200 2202 2X1109 = 2218 (A’) 
2175 2165 1346+833 = 2179 (A’”’) 
2148 (1040) +1109 = 2149 (A’) 
2117 sh 1116+2X513 = 2139 (A’) 
2120 ) Fermi doublet 
2103 2097 0.35 | 
2096 | 2X (1040) 
2071 | a’; (C—D) 
2071 2069 0.36 
2050 
2036 
1353+668 = 2021 (A’”’) 
2018 
1970 
1958 1950 1109+(890) = 1999 (A’) 
1945 
1714 sh ? (810) +(890) = 1700 (A”’) 
1757 
1743 1690 1702 0.21 a’; v(C=0) 
1730 
1442 
1417 1426 0.80 a’; a’; CH; def. 
1420 
1366 
1353 1346 1343 0.68 a’; CH; def. 
~1345 
1244 v.v.w. 2X668 = 1236 (A’) 
1192 v.v.w. ? 
1118 
1109 1116 1111 0.12 a’; C—D rocking 
1094 
Overlapped 1081 1080 0.50 a’; v(C—C) 
1043 v.w. 1040 833 +tors. (?) 
1025 


1023 2X513 = 1026 (A’) 
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TABLE II (Concluded) 














—— Raman Assignment 
Vapor Solid liquid aie. 
924 v.v.w. 668+tors. (?) 
913 
865 
849 ~880 858 0.46 a’; CH; rocking 
837 
816 
802 833 ~820 (overlapped) ? a’’; CH; rocking 
799 

770 ? 

668 672 674 ~0.86 a’; C—D wagging 
513 
500 513 505 0.36 a’; sC—C=O 
487 





DISCUSSION OF RESULTS 

Following previous authors, symmetry C, is assumed for both molecules. 
Each molecule should have ten a’ and five a” fundamental modes of vibration. 
All modes should be infrared and Raman active, the Raman bands arising 
from the a’ modes being polarized and those from the a’’ modes being de- 
polarized. Moments of inertia were calculated* and the structures expected 
for the principal types of infrared bands were determined according to Badger 
and Zumwalt (1), some extrapolation of their results being necessary in this 
case. The a’”’ modes should give rise to bands with strong central peaks (C type) 
while the a’ modes should give rise to hybrid bands, i.e. bands which consist 
of doublets with or without a weak central peak. The separation of the doublets 
was calculated, using the symmetric rotator formulae (5), to be ~27-28 cm. 
for both CH;CHO and CH;CDO. The symmetric top approximation is not a 
very good one in these cases but several doublets of approximately this sepa- 
ration were observed. 


The C—H Stretching Modes 


Comparison of the spectra of the two molecules in the 3000 cm.— region 
gives immediately those bands which arise from the aldehyde CH group and 
the overtone and combination tone levels which are in Fermi resonance with it. 
The infrared spectra of the vapors show overlapping bands which are not 
readily identified individually but the solid spectra and the Raman spectra of 
the liquids are more clearly resolved. The depolarization ratios in the Raman 
spectra distinguish between the a’ and a”’ asymmetric CH; stretching modes at 
2965 and 3001 cm. respectively. Previous assignments in this region are 
markedly different. 


*Dimensions used in calculating the moments of inertia were: C—H=C—D = 1.09 A; C—C 
= 1.50 A; C=0 = 1.22 A; angle CHO = 120°; CH; angles tetrahedral (13). 
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C—D Stretching Mode 


Strong Fermi resonance between the fundamental and an overtone or 
combination tone is evident from the appearance of two almost equally 
intense, similarly polarized bands in the Raman spectrum of liquid CH;CDO. 
The infrared spectrum of the solid shows similar features with further overtone 
or combination bands in the vicinity, but a reduction in the extent of the 
resonance is indicated by the infrared spectrum of the vapor. A band with 
doublet structure of separation ~21 cm. is the pronounced feature present. 
The identity of the level which resonates with the fundamental »(CD) will be 
discussed in conjunction with the assignment of the modes in the 1000 cm." 
region. 


Region below 1800 cm. 


The new data did not warrant any major change in Pitzer and Weltner’s 
assignment of the fundamentals which lie in this region, but they did provide 
support in some instances and confirmation in others. Fig. 1, giving the Raman 
spectra of the two liquids, shows that the in-plane aldehyde C—H rocking 
mode is at 1391 cm.~!, and in CH;CDO it is at 1111 cm.—. In the infrared 
spectra of the vapors, the bands overlap but it is apparent that the most pro- 
nounced difference between the two spectra in the CH; deformation region is 
at ~1410 cm.—!. The aldehyde group CH rocking mode is therefore identified. 

The out-of-plane wagging modes of the aldehyde C—H and C—D groups may 
confidently be assigned to the two depolarized Raman bands at 767 and 
674 cm.—! respectively. Both bands are very weak in the infrared but both 
have the expected C type structure. It has been suggested recently that this 
C—H mode in CH;CHO coincides with the skeletal bending mode near 
500 cm.~, part of the evidence for this being the reported depolarization ratio 
of 6/7 for the Raman band (4). However, this band, which was hardly affected 
by deuterium substitution at the aldehyde C—H group, was observed here to be 
polarized in both molecules. 

The remaining bands, in the 800-1200 cm. region, may be readily corre- 
lated in the two molecules but for an understanding of the differences between 
the spectra it is necessary to postulate coupling of the various vibrational 
modes. 

The C type bands at 867 cm. in CH3;CHO and at 802 cm.“ in CH;CDO 
arise from the one a” mode which is expected in this region, a CH; rocking 
mode. The corresponding Raman bands should be depolarized but in both 
cases overlapping with adjacent polarized bands prevented this from being 
confirmed. The unexpectedly large difference in position of this band in 
CH;CHO and CH;CDO (~65 cm.—') may be explained in the following 
manner. In CH;CHO coupling of the two a’’ modes, C—H wagging and CH; 
rocking, which are both close to 800 cm.~', forces the two energy levels apart, 
one increasing in wave number and the other decreasing. In CH;CDO the 
coupling is removed and the apparent shift in the CH; rocking mode is appreci- 
ably more than the shift expected for an uncoupled mode of this kind. Two 
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modes remain to be assigned in this region, the C—C stretching and the CH; 
rocking modes, both being a’ modes. The situation is further complicated in 
CH;CDO by the presence of the a’ mode, C—D rocking in this region. In 
CH;CHO there are two bands, the higher of which at 1109 cm. in the Raman 
spectrum has been assigned by Pitzer and Weltner to the »(C—C) mode. This 
high value was favored because of the measured, short C—C bond (13). The 
other band, at 911 cm.~' (Raman), is then the CH; rocking mode. 

At first sight, comparison of the spectra of CH;CHO and CH;CDO would 
indicate a reversal of this assignment since the band at 911 cm.~! in CH;CHO 
is at 858 cm.~—' in CH;CDO and the band at 1109 cm. in CH;CHO is at 
1080 cm.—! in CH;CDO, that is, we have the surprising result that deuterium 
substitution has a greater effect upon the CH; rocking mode than upon the 
C—C stretching mode. However, if because of this apparent anomaly, the 
assignments were reversed, a further important feature of the spectrum of 
CH;CDO, the Fermi doublet at ~2080 cm.—', would still remain unexplained. 
Since it is not possible to find a suitable overtone or combination of ‘the ob- 
served bands in CH;CDO to account for this doublet, it must be assumed that 
the overtone of a perturbed fundamental is involved. A plausible suggestion 
is that the 1109 cm. band in CH;CHO arises from »(C—C) and that on deu- 
terium substitution it would fall to ~1050 cm.—!. However, interaction with 
the a’ CH; rocking mode gives rise to two levels at 1080 cm.—' and 858 cm.~', 
and thus the unexpectedly large change in the CH; rocking mode is explained. 
An equally plausible explanation is that the main interaction is between the 
C—D rocking mode and the CH; rocking mode, the band at ~2080 cm. being 
then the first overtone of the C—D rocking mode, which, if uncoupled, would 
lie at ~1050 cm.—. It is noteworthy that in CH;CHO, »v(C—H) is in strong 
Fermi resonance with two overtones or combination bands, one of which is 
very probably the overtone of the C—H rocking mode. However, the assign- 
ments of overtone and combination bands in Table II are based on the first 
explanation, although no serious difficulty would be encountered in using the 
alternative explanation. 

Product rule ratios were calculated* assuming negligible interaction be- 
tween the torsional modes and the other a’’ modes. Theoretical and observed 











TABLE III 
Species Theoretical ratio Observed ratio 
a’ 1.91 1.89 
a”’ 1.20 L322" 





*The observed ratio 1s, of course, dependent upon the un- 
certainties in estimating the position of the unresolved band 
825 cm.~ 


*Moments of inertia used were: 
CH;CHO: I4 = 16.6; Ip.= 78.8; Ic 
CH;CDO: I4 = 20.8; Ig = 


= 90.1X10- gm. cm.? 
79.0; Ic = 94.4X%10-° gm. cm? 
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ratios are collected in Table III. The contribution made by the torsional modes 
has been factorized out of the theoretical value for the a’’ modes (9). 

It may be seen from Tables I and II that the interpretation of several weak 
features required combination with the unobserved torsional mode. Further, 
in the infrared spectrum of solid CH;CDO very weak features at 770 and 
2117 cm.—', and in the vapor at 1192 cm.—', remain unexplained. 


Solid Spectra 

Most bands are at different positions in the solid and vapor spectra, some 
being at higher wave numbers in the solid and others at lower wave numbers. 
These changes provide an indication of the mechanism of intermolecular 
interaction between acetaldehyde molecules. The most significant are the 
changes undergone by the bands arising from the stretching modes of the 
carbonyl group and the aldehyde C—H and C—D groups. v(C—H) and »(C—D) 
are at higher wave numbers in the solid spectra indicating that hydrogen 
bonding between the C—H or C—D and the C—O groups is not involved in the 
interaction. It appears that the interaction is similar to that occurring in 
formaldehyde (11) and that similar arguments may be invoked to explain the 
changes in »(C—H), »(C—D), and the other vibrational modes with phase 
change, complicated, however, by the fact that interaction between vibrational 
modes is present in CH;CDO. 

There is also a difficulty in assigning the strong band at 1040 cm.“ in solid 
CH;CDO. Its considerable intensity suggests that it might be a fundamental. 
However this would lead into a very much more complicated interpretation 
of the Raman bands and the vapor spectrum than that already given. 


TABLE IV 
O 


ASSIGNMENT OF RAMAN SPECTRA OF CH;—C—X 








X =H D 





F Cl Br I 
{vCHs; 2964 2965 2994 p 2991 (3001) 
|yCHs 2917 2917 2940 2938 p 2929 2924 
| 2843) 2097 | 
|yC—X ait 815 594 p 555 p 540 
| 2731/} 2069 
vC=O 1714 1702 1840 1802 p 1812 p 1793 
| CH; 1426 1426 1428 1422 0.9 1418 0.85 1418 
a’ 46C—X 1391 1111 434 or 590 §=352por438p 304por338p 261 or 285 
|5CH3; 1342 1343 1376 1361 p 1352 0.96 1350 
|6CH3; 911 858 (1004) 957 p (938) d? (904) 
lyCC 1109 1080 1106 1100 p 1084 p 1067 
{6CCO 512 505 590 or 434. «438 por352p 338 or 304 285 or 261 
(vCH; 3001 2998 3038 3018 d 3015 
a” )56CHs; 1426 1426 1428 1422 0.9 1418 0.85 1418 
5CH; 885 ~825 (1032) p? (988) (943) 
(6CX 767 674 (208) (260) p? 





Brackets indicate tentative assignments. 


Bands definitely polarized or depolarized are marked p and d respectively. 
Numerical values of depolarization ratios are given for some bands where ambiguities might 


arise, and also question marks. 
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7 J? 
Vibrational Assignment for Molecules CH3C 


xX 


The vibrational assignment for CH;CHO and CH;CDO is of great assistance 
in understanding the spectra of the acid halides for which X is F, Cl, Br, and I. 
In the case of acetyl chloride the Raman spectrum with depolarization ratios 
was obtained with the photoelectric instrument previously mentioned. For 
the fluoride (12), bromide (12), and iodide (7), however, the data were taken 
from the literature and were in all cases Raman data only. Our results for 
acetyl chloride are in good agreement with those of Reference 12 with the 
exceptions that the band at 260 cm.~ is not reported in Reference 12 and that 
the band at 957 cm.—! was found here to be polarized. In Table 1V the spectra 
of these molecules in the liquid phase are assigned, where possible, according 
to the schematic representation given in column 1. 
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HYDROLYSIS AND HYDRAZINOLYSIS OF ESTERS 
OF N,N-DIMETHYLDITHIOCARBAMIC ACID 
A METHOD FOR THE PREPARATION OF MERCAPTANS'! 


By MARSHALL KULKA 


ABSTRACT 


A method has been developed for the preparation of mercaptans from alkyl 
and aralkyl chlorides. This consists of the condensation of sodium N,N-dimethyl- 
dithiocarbamate with halides followed by alkaline hydrolysis or hydrazinoly sis 
of the resulting N,N-dimethyldithiocarbamates (IV). In the hydrazinolysis of IV, 
an insoluble solid by-product was formed which was identified as 3-hydrazino- 
4-amino-5-mercapto-4,1,2-triazole. Although this method was found to have 
wide scope a few failures were encountered. Nitro-derivatives of IV were sensitive 
to alkali or hydrazine and yielded only tars on attempted degradation. B-Phen- 
oxyethyl N,N-dimethyldithiocarbamates (IX) were found to possess a labile 
alkyl-oxygen bond so that in the presence of alkali or hydrazine only the corre- 
sponding phenols were formed. The labile nature of the alkyl-oxygen bond 
ceased abruptly when the alkyl chain of IX was increased so that 3-p-chloro- 
phenoxypropy! and 4-p-chlorophenoxybutyl mercaptans could be prepared in 
high yields. 


In connection with our search for new insecticides a variety of mercaptans 
were required as iittermediates. Of the most useful methods for the preparation 
of mercaptans even the thiourea method (21) is limited in scope. While 
thiourea reacts readily with alkyl and aralkyl bromides to form the iso- 
thiuronium salts, the reaction with many chlorides is extremely slow or does 
not proceed at all. Furthermore no nitrobenzyl mercaptans can be prepared 
by this method. The purpose of this investigation was to find a method for the 
preparation of mercaptans which would supplement the thiourea method. 

v. Braun (5) observed that some aralkyl dithiocarbamates (I), which could 
be easily obtained from the corresponding aralkyl chlorides, when heated alone 
or with alkali gave the corresponding mercaptans (II) in high yields. One 
objection to the use of dithiocarbamates (I) as a source of mercaptans is the 


RS:,CNH,2———> RSH + HSCN 
I II 


fact that they (I) are not stable enough to be purified by distillation. The 
esters of N,N-dimethyldithiocarbamic acid on the other hand may be distilled 
without decomposition and therefore they were adopted for degradation 
studies. 

It was found that sodium N,N-dimethyldithiocarbamate (III), which is 
readily obtained from carbon disulphide, dimethylamine, and sodium hy- 
droxide, reacts easily with chlorides including p-nitrobenzyl chloride, aryl- 
oxyalkyl chlorides, and p-nitrochlorobenzene, compounds which do not react 
with thiourea under normal conditions. The resulting N,N-dimethyldithio- 
carbamates (IV) which were obtained in high yields (Table I) were not only 


1Manuscript received April 4, 1956. 


0 Contribution from the Dominion Rubber Company Limited Research Laboratories, Guelph, 
ntarto. 
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thermally stable but also exhibited considerable resistance to hydrolysis. 
Prolonged heating with aqueous acids or in acetic hydrochloric acid solution 
did not produce any appreciable change while heating with aqueous alkali 
caused only partial hydrolysis. However when IV was heated under reflux for 
five hours or longer with five moles of aqueous-alcoholic alkali, the mer- 
captide (V) was formed in high yield together with sodium sulphide, sodium 
carbonate, water, and dimethylamine. A variety of mercaptans were success- 
fully prepared (Table II) from the corresponding chlorides by this two stage 


RCI + NaS2CN(CHs)2 — RS2CN(CHs)2 + NaCl 
III IV 





IV + 5NaOH —— RSNa + NaS + Na2CO3 + 2H20 + NH(CHs)2 
V 


method and a few failures were encountered. p-Nitrobenzyl N,N-dimethy]l- 
dithiocarbamate (IV, R = p-NO,C,H,CH,) like p-nitrobenzyl isothiuronium 
chloride is sensitive to alkali and attempted hydrolysis yielded only tar. Also 
p-xylylene bis-N,N-dimethyldithiocarbamate, p-(CH3)2NCS.CH2CsH,CH2S,- 
CN(CHs3)2, remained unchanged after prolonged heating with excess alkali. 
Another method of degradation of IV, namely hydrazinolysis, was investi- 
gated. When benzyl N,N-dimethyldithiocarbamate (IV, R = CsH;CH2) was 
heated under reflux with excess hydrazine hydrate in ethanol, an 82% yield 
of benzyl mercaptan (II, R = CsH;CH:) was obtained. In addition a highly 
insoluble white compound precipitated from the reaction mixture and proved 
to be 3-hydrazino-4-amino-5-mercapto-4,1,2-triazole (VIII). This (VIII) must 


RS2CN (CHs)2 + H2NNH2—— RSH + HzNNHCSN(CHs)s 
IV II VI 


VI + HaNNH:2 ——> H2NNHCSNHNH: + NH(CHs)2 
VII 





N 
wal “ 
2vVlI—>_—s_—isf| 1 + H:NNH; + HS 
H.NN CNHNH; 
VIII 


be formed from the thiosemicarbazide (VI) which with hydrazine could 
produce thiocarbohydrazide (VII) and dimethylamine. Two moles of VII or 
alternatively one mole of VI and one mole of VII could undergo cyclization to 
form VIII. That this is the mechanism for the formation of VIII is indicated 
by the fact that authentic samples of diethyl xanthate or thiocarbohydra- 
zide (19, 20) when heated with hydrazine in ethanol formed VIII. 

It is noteworthy that where alkaline hydrolysis failed, hydrazinolysis suc- 
ceeded in the case of p-xylylene bis-N, N-dimethyldithiocarbamate. p-Xylylene 
dimercaptan was obtained in 77% yield after prolonged treatment with hydra- 
zine. However, attempted hydrazinolysis of p-nitrobenzyl and p-nitrophenyl 
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N,N-dimethyldithiocarbamates like attempted alkaline hydrolysis yielded 
only resinous materials. 

In view of the known anomalous reactions of some phenoxyalkyl halides it 
was necessary to investigate the conversion of these halides to the phenoxyalkyl 
mercaptans. Gordon (10) and Lyle, Dewitt, and Pattison (16) observed that 
2-phenoxyethyl halides did not produce the Grignard reagents in the presence 
of magnesium and ether. Instead the alkyl-oxygen bond was severed and 
phenol was formed. On the other hand phenoxyalky! halides with alkyl chain 
greater than three carbon atoms behaved normally in respect to Grignard 
formation (16). Also, abnormal results were obtained in the reaction of 2- 
phenoxyethyl bromide with sodium cyanide but not in the case of 3-phenoxy- 
propyl bromide which yielded y-phenoxybutyronitrile (10). Crane and 


- 
Rydon (6) have found that compounds of the type CsH;0CH2CH2SRz, but 
not of type CsH;0CH:CH2SR, are split by dilute alkali to give initially 


+ 
R,SCH=CH; and sodium phenoxide. 

In this investigation all phenoxyalkyl halides reacted normally with III to 
give the corresponding phenoxyalkyl dimethyldithiocarbamates (Table 1). 
However when the latter were subjected to alkaline hydrolysis or hydra- 
zinolysis it was soon found that 2-phenoxyethyl N,N-dimethyldithiocarba- 
mates (IX) gave no mercaptans. Instead the alkyl-oxygen bond of IX was 
severed and the corresponding phenol (X) was obtained in high yield. The 


¢_NOCH:CHS:CN(CH) +H:NNH:—> @ SOH + (CH:CH:S)x + NH(CH:): 
IX x 


+ 


same results were obtained when §-p-t-butylphenoxyethyl bromide was treated 
with thiourea and the resulting isothiuronium salt heated with aqueous alkali. 
The sulphur-alkyl bond in the compound obtained by substituting sulphur 
for the oxygen atom in IX is also labile in nature. Thus the reaction of. 6-p- 
chlorophenylmercaptoethyl N,N-dimethyldithiocarbamate with hydrazine in 
ethanol gave p-chlorothiophenol in 50% yield. 

It was expected that steric hindrance might play a part in retarding the 
cleavage of the alkyl-oxygen bond in IX. Therefore 8-(2,4-dichlorophenoxy) 
ethyl and 6-(2-t-butyl-4-methylphenoxy)ethyl N,N-dimethyldithiocarbamates 
were prepared and then subjected to hydrazinolysis. The products of the re- 
actions were 2,4-dichlorophenol and 2-t-butyl-4-methylphenol respectively 
showing that the ortho bulky groups do not strengthen the alkyl-oxygen bond. 

The weak nature of the alkyl-oxygen bond of IX ends abruptly when the 
length of the alkyl chain is increased. Thus 3-p-chlorophenoxypropyl and 4-p- 
chlorophenoxybutyl mercaptans were prepared in 85 and 87% yields by the 
hydrazinolysis of the corresponding phenoxyalkyl N,N-dimethyldithiocarba- 
mates. There is also a difference in the light stability of IX and its higher 
homologues. The phenoxypropyl and phenoxybutyl N,N-dimethyldithiocarba- 
mates remain unchanged on standing while IX turns pink or red very quickly 
on exposure to light. 
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It has been already reported (14) that when o0-(2-chloroethoxy)benzyl 
chlorides (XI) are treated with thiourea and the resulting o-(2-chloroethoxy) 
benzylisothiuronium chlorides subjected to alkaline cleavage, no mer- 


=f) OCHCHC! a OCH.CH:S:CN(CH;)s 
’ | 
LU 


ate CHCl CH2S,CN (CHs): 


captans are formed because simultaneous cyclization occurs to give the seven- 
membered ring compounds the 2,3-dihydrobenzo(f)-1,4-oxathiepins. With 
sodium dimethyldithiocarbamate (III), 2-(2-chloroethoxy)-4-methylbenzyl 
chloride (XI, X = 4-CH3) vielded the bis-dithiocarbamate (XII) but alkaline 
cleavage of this did not yield a dimercaptan nor 7-methyl-2,3-dihydrobenzo(f)- 
1,4-oxathiepin. Instead a complex mixture was formed whose components 
were pot identified. 


XII 


EXPERIMENTAL 


Sodium N,N-Dimethyldithiocarbamate 

To stirred 40% dimethylamine (340 ml.) was added dropwise first carbon 
disulphide (185 ml.) and then sodium hydroxide (120 gm.) in water (300 ml.), 
the reaction temperature being kept at 10—15° by cooling. The reaction mixture 
was stirred at room temperature for two hours and then concentrated in vacuo 
to about half the original volume and allowed to cool. The precipitate was 
filtered, washed with a little cold acetone, and air-dried. The yield from three 
successive crops was 361 gm. This product, which is a hydrate melting at 
106—108°, was pure enough to be used in the preparation of the N, N-dimethyl- 
dithiocarbamic esters. 

Preparation of the N,N-Dimethyldithiocarbamates (Table I) 

To a solution of sodium N,N-dimethyldithiocarbamate (0.3 mole) in acetone 
(400 ml.) was added the halide (0.25 mole), and the reaction mixture was 
heated under reflux for periods depending upon the reactivity of the halide. 
With alkyl bromides and benzyl chlorides the reaction was exothermic and 
required little heating. Three hours’ refluxing time was required in the case 
of the phenoxyalky] chlorides, and in the extreme case of p-nitrochlorobenzene 
several days were required. The precipitated sodium halide was filtered, the 
acetone was distilled off from the filtrate, and the residue was either distilled 
or crystallized from methanol. 


Hydrolysis of N,N-Dimethyldithiocarbamates to Mercaptans (Table II) 

This is exemplified by the hydrolysis of benzyl N,N-dimethyldithiocarba- 
mate to benzyl mercaptan. To a solution of sodium hydroxide (15 gm.) in 
water (30 ml.) was added benzyl N,N-dimethyldithiocarbamate (7) (15 gm.) 
and ethanol (150 ml.), and the reaction mixture was heated under reflux for 
five hours. Dimethylamine was given off through the condenser and could be 
caught in a dry ice trap, and sodium carbonate precipitated from the reaction 
mixture. The ethanol was removed and the residue acidified. Carbon dioxide 
and hydrogen sulphide were liberated. The precipitated mercaptan was ex- 
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tracted with benzene, the extract washed with water, and the solvent removed. 
The residue distilled at 77-79° at 12 mm. pressure yielding a colorless odiferous 
liquid (7.2 gm. or 86%) of benzyl mercaptan. 
Hydrazinolysis of N,N-Dimethyldithiocarbamates to Mercaptans (Table II) 
This is exemplified by the hydrazinolysis of 3-p-chlorophenoxypropyl 
N,N-dimethyldithiocarbamate to 3-p-chlorophenoxypropyl mercaptan. A solu- 
tion of 3-p-chlorophenoxypropyl N,N-dimethyldithiocarbamate (12 gm.), 85% 
hydrazine hydrate (20 ml.), and ethanol (80 ml.) was heated under reflux for 
30 hr. Dimethylamine was given off through the condenser and a white solid 
precipitated from the reaction mixture. The white solid (1.1 gm.) was filtered, 
washed with methanol, water, methanol, and hot benzene, and dried. It melted 
at 233-235° with decomposition and did not depress the melting point of 
3-hydrazino-4-amino-5-mercapto-4,1,2-triazole (20). The filtrate was taken 
to dryness 7m vacuo. The residue was treated with water and extracted with 
benzene. The extract was washed with water and the solvent removed. The 
residue distilled at 153-155° at 10 mm. pressure as a colorless liquid with odor 
not too objectionable. The yield was 7.2 gm. or 85%, n3* = 1,5620. Anal. calc. 
for CghvOCis: C, 53.43: Hi, 5.43. Found: C, 53.37, 53.37; 1, 5.41, 5.14. 


TABLE II 
THE PREPARATION OF MERCAPTANS, RSH, FROM THE CORRESPONDING 
N,N-DIMETHYLDITHIOCARBAMATES, RSeCN(CH3)2, BY HYDROLYSIS 
(METHOD A) AND BY HYDRAZINOLYSIS (METHOD B) 











Method of ne 

Mercaptan, R = prep. B.p. " % Yield 

B 73 (10 mm.) 82 
p-CICsHsCH;z (15) A 105 (10 mm.) 95 
2,4-DiCICsH;sCHe (4) A 123-125 (10 mm.) 1.6000 88 
3,4-DiCICsH3sCHs (15, 4) A 132-134 (11 mm.) 1.6015 93 
2,4-Di(CH3)CsH3;CHe2) (13) Ae 100-102 (10 mm.) 1.5610 70 
C.HsCH2CHz (5, 12) A 95- 98 (12 mm.) 66 
CH;(CH2),CHe2 (2) A* 155-160 61 
CH;(CHe) uCH2 A* 182-188 (10 mm.) 65 
p-CICs5H,O(CH 2); Be 153-155 (10 mm.) ‘1.5620 85 
p-CIC,H,O(CHe), Be 170-172 (10 mm.) 1.5550 87 
p-SHCH2C.HsCH2SH (1) Bt M.p. 44-45 re f 





*The reaction time required was 24 hr. instead of 5. 

{The reaction time required was three days instead of 30 hr. This compound was 
also prepared in 86% yield by the thiourea method: 
ANALYSES 

*Calc. for CosHi2S: C, 71.05; H, 7.89. Found: C, 71.01; H, 7.78. 

’Cale. for CsHnOCIS: C, 53.338; H, 5.43. Found: C, 53.37; H, 5.41. 

‘Calc. for CyoHyOCIS: c 55.43; Bi. 6.00. Found: Cc. 55.41; A. 5.85. 


Preparation of Sulphides 

This is exemplified by the preparation of 3-p-chlorophenoxypropyl p-chloro- 
benzyl sulphide. 

To a solution of potassium hydroxide (0.7 gm.) in methanol (40 ml.) was 
added p-chlorobenzy! chloride (2.0 gm.) and 3-p-chlorophenoxypropyl mer- 
captan (2.3 gm.), and the resulting solution was heated under reflux ‘for 
10 min. The precipitated potassium chloride was filtered, the methanol re- 
moved from the filtrate, and the residue distilled. The colorless liquid (3.2 gm.) 











KULKA: HYDROLYSIS AND HYDRAZINOLYSIS 





1099 


which boiled at 257° at 10 mm. pressure solidified on cooling. The cake was 
pulverized, washed with methanol, and dried. The white product melted at 
37-38° alone or in admixture with the compound prepared from p-chlorobenzy1 
mercaptan and 3-p-chlorophenoxypropy! bromide. Anal. calc. for CisHigOC1.S: 
C, 58.72; H, 4.89. Found: C, 58.72, 58.78; H, 4.64, 4.95. 

The above method was also used for the preparation of the following 
sulphides: 

Bis-p-chlorobenzyl sulphide from p-chlorobenzyl mercaptan and p-chloro- 
benzyl chloride in 80% yield. It melted at 42-43°; literature (17) 41°. 

Benzyl 3-nitro-4-methoxybenzyl sulphide from benzyl mercaptan and 3-nitro- 
4-methoxybenzyl chloride in 78% yield, m.p. 55-56°. Anal. calc. for 
C,sHisNO;S: C, 62.30; H, 5.19. Found: C, 62.30; H, 5.20. 

p-Chlorobenzyl 8-p-t-butylphenoxyethyl sulphide from p-chlorobenzyl mer- 
captan and 6-p-t-butylphenoxyethy! bromide in 90% yield, b.p. 250° (10 mm.). 
Anal. calc. for CigH230CIS: C, 68.17; H, 6.88. Found: C, 68.07; H, 6.86. 

p-Chlorobenzyl n-dodecyl sulphide from n-dodecyl mercaptan and p-chloro- 
benzyl chloride in 80% yield, b.p. 235-236° (10 mm.), 22° = 1.5155. Anal. 


D 


calc. for CigH3:ClS: C, 69.82; H, 9.50. Found: C, 69.89; H, 9.35. 


Hydrazinolysis of B-2,4-Dichlorophenoxyethyl N,N-Dimethyldithiocarbamate to 
2,4-Dichlorophenol 

A reaction mixture of -2,4-dichlorophenoxyethyl N,N-dimethyldithio- 
carbamate (20 gm.), 85% hydrazine hydrate (20 ml.), and ethanol (100 ml.) 
was heated under reflux for 30 hr. The white precipitate was filtered from the 
cooled reaction mixture and the filtrate taken to dryness in vacuo. The residue 
was treated with water and extracted with benzene. The benzene extract was 
washed with water and the solvent removed. The residue distilled at 89-93° 
at 10 mm. pressure and the distillate (9 gm. or 86%) which solidified melted 
at 38—40° alone or in admixture with 2,4-dichlorophenol. The white precipitate 
was extracted with boiling benzene leaving a benzene-insoluble product 
(1.5 gm.) which melted at 233-235° with decomposition and did not depress 
the me!ting point of 3-hydrazino-4-amino-5-mercapto-4,1,2-triazole (20). The 
benzene extract was concentrated and cooled. A white amorphous malodorous 
substance (0.4 gm., m.p. 80—90°) precipitated which gave analyses agreeing 
with those of ethylene sulphide polymer. 

In the same way §6-p-chlorophenoxyethyl, B-p-t-butylphenoxyethyl, B-p-tolyl- 
oxyethyl, and B-(2-t-butyl-4-methylphenoxy)ethyl N,N-dimethyldithiocarbamates 
hydrazinolyzed to p-chlorophenol, p-t-butylphenol, p-cresol, and 2-t-butyl-4- 
methylphenol respectively. 


Alkaline Cleavage of 8-p-t-Butylphenoxyethyl Bromide and N,N-Dimethyldithio- 
carbamate 

A solution of 8-p-t-butylphenoxyethyl bromide (55 gm.), thiourea (20 gm.), 
and ethanol (200 ml.) was heated under reflux for 24 hr. To this was added 
sodium hydroxide (20 gm.) in water (100 ml.) and the resulting reaction 
mixture was heated under reflux for 5 hr. The cooled reaction mixture was 
diluted with water, acidified, and extracted with benzene. The benzene extract 
was washed with water, the solvent removed, and the residue distilled. There 
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was obtained an alkali-soluble distillate (22 gm.), b.p. 120-125° (11 mm.), 
which when crystallized from benzene melted at 98-100° either alone or in 
admixture with p-f-butylphenol. This phenol was also obtained when £-p-t- 
butylphenoxyethyl N,N-dimethyldithiocarbamate was heated with alcoholic 
alkali as above. 
3-Hydrazino-4-amino-5-merca pto-4,1 ,2-triazole 

A reaction mixture of thiocarbohydrazide (2 gm.) (11), ethanol (10 ml.), and 
85% hydrazine hydrate (5 ml.) was heated under reflux for 24 hr. The white 
precipitate was filtered, washed with methanol, water, methanol, and hot 
benzene in that order. The white powder (0.85 gm.) melted at 233-235° with 
decomposition. Stollé and Bowles (20) give a melting—decomposition point 
of 228°. A similar prolonged treatment of diethyl xanthate also gave 3-hydra- 
zino-4-amino-5-mercapto-4, 1,2-triazole. 
4-p-Chlorophenoxybutyl Chloride 

A reaction mixture of p-chlorophenol (65 gm.), 1,4-dichlorobutane (100 ml.), 
a solution of potassium hydroxide (30 gm.) in water (100 ml.), and benzene 
(200 ml.) was heated under reflux for three days. The organic layer was 
separated, washed with aqueous alkali and with water, and the solvent and 
excess 1,4-dichlorobutane were removed in vacuo. The residue distilled at 
159-162° (10 mm.) as a colorless liquid, 22? = 1.5340, yield 45 gm. Anal. calc. 
for CipH»OCI.: C, 54.80; H, 5.48. Found: C, 54.35, 54.71; H, 5.26, 5.39. 

The residue from the distillation when crystallized from benzene—methanol 
yielded 21 gm. of white prisms of 1,4-bis-p-chlorophenoxybutane, m.p. 107-108°. 
Anal. calc. for CjgH:g02Cle: C, 61.64; H, 5.14. Found: C, 61.74; H, 5.12. 
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THE APPLICATION OF INTERPOLATION THEORY TO DATA 
ON THE ADSORPTION OF NITROGEN, OXYGEN, 
AND ARGON ON RUTILE! 


By J. M. Honic anp L. H. REYERSON? 


ABSTRACT 


The interpolation theory has been applied to data on the adsorption of nitrogen, 
oxygen, and argon on rutile at 77° K. It has been demonstrated that the oxygen 
and argon data cannot be interpreted on the basis of the Langmuir model, 
modified by the introduction of surface heterogeneity effects. The interpolation 
theory verifies, however, that the distribution of adsorption energies among 
surface sites, obtained by inversion of the nitrogen isotherm, is physically 
significant, and that nitrogen is suitable for BET surface area calculations. 


INTRODUCTION 


The inadequacy of the Langmuir model in the interpretation of gas adsorp- 
tion phenomena is well known. Nevertheless, a variant of this model, the BET 
multilayer adsorption theory, has found widespread use in the estimation of 
surface areas, though it is recognized that the results cannot be completely 
reliable. Attempts have been made to improve on the Langmuir model by 
taking account of the energy heterogeneity of the surface (2, 3, 10, 11, 12, 13, 
17, 20, 21, 22, 23). While this is a step in the right direction, the results are 
still unsatisfactory because lateral interactions are left out of consideration. 
On the other hand, a reasonably rigorous treatment of both effects (12) or 
even a rather specialized approach along these lines (1, 8, 9). leads to com- 
plicated equations. In practice one is therefore often forced to return to the 
use of the BET equation in the analysis of gas adsorption isotherms. 

In order to have confidence in the results of such a simple approach, one must 
determine how well the data conform to the restrictive assumptions of the 
Langmuir or BET model. Thermodynamic considerations (see for example (4, 
5, 16)) have been quite helpful for this purpose. Recently a new approach has 
been cited in the literature (14, 15) in which it was shown that all data con- 
forming to the Langmuir model, the BET model, or to a theory dealing with 
surface heterogeneity effects must lie within well-defined limits on a plot of 
surface coverage versus pressure. Details for the construction of these limits 
by use of the so-called interpolation theory (18) have been provided. Such an 
analysis should aid in deciding whether the simple adsorption theories are 
applicable to a given adsorbate—adsorbent system. For, if a given set of datum 
points is found to fall outside the interpolation limits, then the Langmuir model 
and the various modifications which include surface heterogeneity effects are 
hopelessly inadequate. In such cases it is obviously not justifiable to use the 

1Manuscript received April 12, 1956. 

Contribution from the Department of Chemistry, Purdue University, Lafayette, Indiana, and 
the Department of Chemistry, University of Minnesota, Minneapolis 14. Presented in part 
before the Division of Colloid Chemistry at the 128th Meeting of the American Chemical Soctety, 
at Minneapolis, September 11-16, 1955. 

2One of us (LHR) gratefully acknowledges funds received from the graduate school in partial 


support of this work. Both authors are indebted to Mr. Chi Chang for performing the calculations 
on which Figs. 1, 2, and 8 are based. 
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isotherm data to compute a distribution of adsorption energies among surface 
sites, and even the BET method of determining at what point on the isotherms 
monolayer coverage is complete is open to criticism. 

The purpose of the present article is to provide examples showing how the 
interpolation procedure operates. The method will be utilized to settle a 
problem raised by Honig and Reyerson (13) in connection with a study on the 
adsorption of nitrogen, oxygen, and argon on rutile. 


THE INTERPOLATION THEORY 


Any adsorbent—adsorbate system which conforms to the modified Langmuir 
theory is characterized by an adsorption isotherm of the form (11): 


n _ PA dF(e) 
= 0D = fae Te 


where » and 7, are the number of adsorbate molecules held on the surface at 
a pressure p and at monolayer coverage respectively, @ is the fraction of the 
surface covered at pressure p and temperature T, dF(e) = f(e)de is the fraction 
of sites associated with adsorption energies in the range ¢ to e+de, €) and em 
are respectively the minimum and maximum values of ¢ encountered on the 
surface, and @ is a quantity defined by: 





(2] B(e, T) = (2umkT/h?)®?kT j,(T)/js(e, T) 


where m is the mass of the adsorbate particles, while 7, and 7, represent the 
partition functions of the particles in the gas and in the adsorbed phases 
respectively. Note that the latter quantity depends on e (7, 11, 12). 

Eq. [1] represents a one-to-one correspondence between @ and F. Rosen- 
bloom (18) has developed a theory in which he proves that one can select 
two distribution functions F,(e) and Fy(e) such that the corresponding equa- 
tions 6.(p, T) and @(p, T) form the limits, mentioned in the introduction, 
within which the data must fall if eq. [1] is to apply. For details of the con- 
struction of the extremal isotherms the reader is referred to the original 
publication (14). 

In what follows it suffices to note that the bounding isotherms 6, and 6 are 
adjusted to a given set of experimental data by the arbitrary selection of a set 
of reference points through which the functions 6, and 6 are required to pass. 
One then observes whether the remaining experimental points lie within the 
region which is bounded by the functions 6, and 6. This gives an indication as 
to the applicability of the Langmuir model or of models combining the Lang- 
muir postulates with the concept of surface heterogeneity. For, it is evident 
from the remarks that if any points fall outside this region, this is a clear 
indication that eq. [1] cannot be used in the interpretation of the data. 

The calculations cited below were based on the use of the following relations: 


“ .(d, T) = (62.—0:)p+[p2(1—82) — pi(1—@1)] , 





7) — ——__0102(Pi— 2) 
4] %(P, T) = Pil2(82—81) /p+ (O2p1—Oipe) © 








ar 
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These equations require the selection of two reference points (6, p:) and 
(2, P2) from among the data. It was shown in (14) that [3] and [4] result upon 
simplification of the expression for the “‘third-order’’ limiting isotherms. Being 
approximate, they bound a wider region in the @ vs. p diagram than the exact 
functions, and thus less precise information is obtained. On the other hand, no 
unknown parameters are involved, and the calculations are relatively simple. 

One should note that the quantity 7, cancels from both sides of [4] when 6 
is replaced by ”/nm. On the other hand, the value of m,, must be specified for 
use in [3]. It has been verified, however, that the dependence of 0, on nm» is not 
marked. It has also been established that none of the conclusions reached below 
would be invalidated by the selection of 7, values which differ radically from 
the set adopted below. 


APPLICATION TO NITROGEN, OXYGEN, ARGON ADSORPTION ISOTHERMS 

The interpolation theory will now be applied to the adsorption on nitrogen, 
oxygen, and argon observed at 77° K. (13). The isotherms are drawn as dotted 
lines through the experimental points in Figs. 1, 2, and 3. Two sets of extremal 
isotherms are entered as heavy curves on each figure. The reference data used 
in the computation of each set are listed in the captions of the diagrams. The 
areas bounded by each set of extremals are shaded and are labelled A and B 
respectively. The low pressure region of Fig. 1 has been expanded in the inset 
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Fic. 1. Limiting isotherms for the adsorption of nitrogen on rutile at K. Reference 
points for Area A: m = 0.125 mmoles/gm.; ~: = 0.13 cm. Hg; m2 = 0.149 mmoles/gm.; 
p: = 0.78 cm. Hg. Reference points for Area B: m = 0.137 mmoles/gm.; p: = 0.29 cm. Hg; 
n. = 0.177 mmoles/gm.; po = 4.12 cm. Hg; mm = 0.179 mmoles/gm. Open circles are experi- 
mental points, connected by dotted line; shaded regions are areas bounded by extremal iso- 
therms. 
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Fic. 2. Limiting isotherms for the adsorption of oxygen on rutile at 77° K. Reference 
points for Area A: m = 0.0346 mmoles/gm.; p; = 0.040 cm. Hg; mz = 0.0901 mmoles/gm.; 
p2 = 0.34 cm. Hg. Reference points for Area B: m = 0.0346 mmoles/gm.; f: = 0.040 cm. Hg; 
nz = 0.142 mmoles/gm.; p2 = 1.13 cm. Hg; am = 0.171 mmoles/gm. Open circles are experi- 
mental points, connected by dotted line; shaded regions are areas bounded by extremal iso- 
therms. 


to Fig. 1. The values of m, were taken from Table II of (13). Attention is 
directed to the interlacing structure of 6, and 6, which alternately form lower 
and upper bounds in the intervals marked off by the reference points p; and p>. 

Inspection of the graphs shows that the low pressure points fall within the 
bounds prescribed by the third-order theory. However, in every case there 
exists a point labelled C beyond which all experimental points fall outside the 
area A. These findings can be qualitatively related to the results of the BET 
theory which predicts the onset of multilayer formation beyond the point 
labelled M on each isotherm.* The formation of higher layers clearly violates 
the Langmuir postulate which excludes the case of multilayer adsorption. It 
is therefore necessary for the high pressure points to fall above the upper limit 
prescribed for the first layer by the interpolation theory. 

Further examination indicates, however, that the isotherms cross the the- 
oretical limits at a point C below the point M where monolayer coverage is 
essentially complete. Thus, for argon the point C occurs at = 0.10 mmoles 
(6 = 0.62), and for oxygen, at » = 0.11 mmoles (6 = 0.64). On the other 
hand, the nitrogen isotherm remains within the limiting bounds up to m = 0.16 
mmoles (6 = 0.90). 


3The presence of multilayer condensation is clearly established by unpublished work of the 
authors in the partial pressure range 0.1 < p/po< 0.9. All three adsorbates yield type II isotherms 
in the BET classification. 
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Fic. 3. Limiting isotherms for the adsorption of argon on rutile at 77° K. Reference points 
for Area A: m = 0.0264 mmoles/gm.; ~: = 0.08 cm. Hg; mz = 0.0808 mmoles/gm.; p: = 
0.69 cm. Hg. Reference points for Area B: m, = 0.0264 mmoles/gm.; pi = 0.08 cm. Hg; 
nz = 0.142 mmoles/gm.; 2 = 2.32 cm. Hg; 2m = 0.164 mmoles/gm. Open circles are experi- 
nalts points, connected by dotted line; shaded regions are areas bounded by extremal iso- 
therms. 


These findings can be checked by the selection of other reference points. 
It has been verified that when these are chosen from the data below the 
point C, the isotherm is always found to cross the third-order limiting curve 
near point C. However, if one of the reference points is selected from the data 
falling above point C, the limiting curves bound an area B which does not 
coincide with the experimental points except near (6, p;) and (62, p2). 


DISCUSSION OF RESULTS 


The results indicate that all adsorption data remain within their extremal 
bounds at low surface coverage, provided that the reference points are selected, 
in conformity with the simplifying assumptions of the interpolation theory, 
from the region below point C. 

However, while the nitrogen isotherm remains within its prescribed bounds 
until monolayer coverage is almost complete, the oxygen and argon isotherms 
deviate from their prescribed limits much sooner. Evidently the system nitro- 
gen: rutile conforms more closely to the modified Langmuir model than the 
other two systems. This behavior can probably be explained on the basis of a 
statement by Drain (6). Nitrogen possesses a quadrupole moment of apprecia- 
ble magnitude, which contributes to the physical binding of nitrogen molecules 
to the surface. In oxygen or argon this contribution is absent, and the binding 
energy to the surface is appreciably lower. As a result, factors such as lateral 
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interactions and surface mobility, which are ignored in the Langmuir model, 
play a more significant role for oxygen or argon than for nitrogen. 

Undoubtedly the experitnental conditions are such that even for nitrogen 
the modified Langmuir model is a poor approximation to reality. Nevertheless, 
the interpolation theory at least does not invalidate the use of the inversion 
technique in this case. 

These conclusions settle a problem raised by Honig and Reyerson (13). 
The authors showed that if the inversion procedure of Sips (19) is applied to 
the data plotted in Figs. 1-3, one obtains three different distribution functions. 
In particular, the rutile surface seemed to present an abundance of high energy 
sites toward nitrogen, but not toward oxygen or argon. Presumably, if the 
concept of energy heterogeneity has physical reality, the distribution functions 
obtained with the three gases on the same solid should have been similar. The 
interpolation theory now indicates that the restrictive assumptions basic to 
the inversion procedure do not hold for the oxygen or argon isotherms. Hence, 
in the light of information now available, the energy distributions calculated 
in (13) from the latter two isotherms have no physical significance. 

Perhaps one of the most interesting results is that the interpolation theory 
confirms the suitability of using nitrogen gas for surface area measurements of 
rutile according to the BET procedure. The selection of oxygen or argon as 
the adsorbate at 77° K. by Brunauer, Emmett, and Teller would obviously 
have been much less satisfactory for surface area studies. The interpolation 
method also indicates that the magnitude for n» obtained by the BET method 
must be approximately correct. For, it may be assumed that the failure of the 
nitrogen datum points to lie within the area A of Fig. 1 is due primarily to the 
onset of multilayer formation. Hence, the fact that the points C and M agree 
within 10% lends credence to the BET results. 
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THE ENTROPIES OF IONS IN AQUEOUS SOLUTION 
I. DEPENDENCE ON CHARGE AND RADIUS! 


By K. J. LAIDLER 


ABSTRACT 


Empirical equations recently proposed for the entropies of ions in aqueous 
solution involve a linear dependence on the first power of the valency z and, in 
some cases, an inverse dependence on the square of a modified radius. Such 
relationships have no fundamental validity since the entropies are related 
arbitrarily to a zero value for the proton. It is shown that the entropy changes 
for reactions of various ionic types indicate that the entropy must depend on ?, 
as proposed by Born on the basis of simple electrostatic theory. When the 
entropies of monatomic ions are related to the absolute scale, in which the proton 
has a value of —5.5 e.u., they are found to vary linearly with 2?/ru, where ru is 
the univalent radius defined by Pauling. The slope is found to be close to that 
predicted by the Born equation, and the non-electrostatic part of the entropy is 
shown to correspond to the movement of the ion in a free volume of 0.73 cu. A. 


INTRODUCTION 
Recently several workers (12, 9, 3, 4, 6) have proposed empirical equations 
which accurately relate the conventional partial molal entropies of ions in 
aqueous solution to parameters such as the valency, z, the effective radius, r,, 
and, in the case of complex ions, the number of ligands. The prototype of such 
equations is that of Powell and Latimer (12) who show that the entropies of 
monatomic ions are accurately represented by the formula 


1] §° = (3/2) R In(at.wt.) +37 —270 2/r2 


where the effective radius 7, is the crystal radius plus 2.0 A in the case of 
cations and 1.0 A in the case of anions. The discussion that follows will relate 
explicitly to the Powell-Latimer scheme but the arguments are applicable, 
mutatis mutandis, to all of the others. 


DEPENDENCE OF ENTROPY ON CHARGE 


Although it is undoubtedly the case that the equations permit accurate 
estimates to be made of the entropies of ions of various kinds, they unfortu- 
nately do not correctly represent the manner in which entropies depend upon 
charge and radius. That this is the case becomes apparent from a consideration 
of the following points: 

(1) The entropies employed in the empirical equations are arbitrarily 
related to the value of zero for the proton. Gurney (8) has recently obtained 
convincing evidence that the absolute entropy of the proton is actually 
—5.5 e.u., and if this is the case the entropies of ions of various types must 
all be modified by the subtraction of 5.5z in the case of cations and the addition 
of 5.52 in the case of anions. With these changes Powell and Latimer’s straight 
line plot splits into several parallel straight lines, and a different dependence 
of S° on z is clearly called for. 


1Manuscript received April 24, 1956. 
Contribution from the Department of Chemistry, University of Ottawa, Ottawa, Canada. 
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(2) Born (2) has proposed that the free energies and entropies of ions are 
the values given by simple electrostatic theory for a conducting sphere of 
charge ze and radius r. The resulting expressions, and the numerical values 
corresponding to ions in aqueous solution at 25° C., are 


[2] F°., = 2’e?/2r D = 2.05 2?/r kcal. per mole, 
[3] S°o1 = (2%e?/2r D) (6 In D/é T) = 9.42 2?/r cal. per degree. 


(In using the above numerical coefficients, the value of r must be in A.) On 
the basis of Eq. [3] one would expect a linear dependence of S° on 2?/r, rather 
than upon 2/r,.2 as found by Powell and Latimer. This is not, of course, a 
compelling argument, since Born’s assumption that ions behave like conducting 
spheres may not be justified. However it will now be shown that there is strong 
experimental evidence in favor of a dependence on 2?. 

(3) Since the variation of entropy with charge is much more important 
than thé variation with size (and this is especially true of the absolute entro- 
pies) one can, to a rough approximation, ignore size effects and see how entropy 
changes in complete reactions depend on the charges on the ions. Thus if one 
considers a reaction 


A?74+B7B = C?c+D?7p 


the entropy change, if the individual entropies depend on 2’, will be related to 
the quantity 
20° + 2p? —24?— Zp" 


which may be written as }-2?. If on the other hand the entropies depend on z 
the entropy change will be related to 


Zct+2Zp—Za— 2p 


which will be written as }>z. There are a number of reactions for which }>2? 
and >} -z differ in sign, and these should provide a means of differentiating 
between the two possibilities; apart from this, one can see whether trends in 
AS can be correlated with }-2? or with }vz. 

The results for a number of reactions are shown in Table I. The AS values 
were mostly calculated using the S° values listed by Latimer; these 8° values 
are all conventional ones, based on S.+° = 0.0, but the values of AS are 
absolute and do not depend on this convention. Only a selection of possible 
examples can be given, but those chosen are typical and show the general 
trends. It is clear that there is a close correlation between AS and }>2?, but 
none at all between AS and }-z. Fig. 1 shows a plot of AS against }>2?. If 
attempts were made to predict AS on the basis of }-z, the wrong sign would be 
obtained in many cases (e.g. reactions 11 to 14). One could, however, make 
rough estimates of AS by multiplying }>2? by about 10. 

This close correlation between AS and }-2’, covering a wide range of reaction 
types, and the lack of correlation with }-z, strongly indicates that entropy is 
a function of 2? and not of z. If the relationships proposed by Powell and 
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TABLE I 
ENTROPY CHANGES FOR REACTIONS OF VARIOUS IONIC TYPES 











AS 

Reaction >: 2 (e.u.) Ref. 

1. CSH;OH + OH- — C,H;O- + H:O 0 0 —- 82 10 

2 C.H;NH2 + H;0+* — C.H;NH;* + H.O 0 0 1.4 10 

3. 2H.0 — H;0* + OH- 2 2 — 18.7 9, 10 

4. NH,OH — NH,* + OH- 2 2 — 18.9 9 

5. HeCO; + H2O — H;0*+ + HCO;- 2 2 — 22.9 f) 

6. HSO,- + HO — H;0* + SO,2- 2 4 — 26.3 y 

7. HCO;- + HO — H;0* + CO,?- 2 4 — 35.0 9 

8. HPO,~ + HO — H;0+ + PO,'- 2 6 — 43.0 9 

9. UOH** — U**+ + OH- 2 8 — 52.0 1 

10. ClO,- + 2H* + 2Fe?+ — CIO;- + H.O + 2Fe** 0 8 — 73.5 9 
11. MnO, + 5Fe?* + 8H*+ — Mn?* + 5Fe**+ + 4H,0 —2 20 —213.6 9 
12. CroO7- + 3TI* + 14H* — 2Cr*+ + 3TI**+ + 7H2O0 —4 24 — 248.6 9 
13. Cr2O0;2- + 6Fe?+ + 14H* — 2Cr**+ + 6Fe*+ +7H:O —4 30 — 339.2 +) 
14. 2MnO,- + 5TI* + 16H*+ — 2Mn?* + 5TI** + 8H.O0 —4 30 —359.2 9 











-100 


4s 


-200 


-300 














Fic. 1. Plot of entropy change against >: for reactions of various ionic types. The 
numbers relate to the reactions listed in Table I. 


Latimer were accepted as of fundamental significance it would be necessary 
to regard the correlations shown in Fig. 1 as purely coincidental, arising from 
fortuitous variations in ion radii. Such would appear to be extremely unlikely 
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A NEW EMPIRICAL EQUATION 

In view of these conclusions it is clearly of interest to examine the absolute 
entropies of aqueous ions with a view to confirming their dependence on 2? 
and to finding out how they vary with ion size. 

Column 2 of Table II shows values for monatomic cations,* calculated 
from Latimer’s values in the manner indicated above. In the next column the 
value (3/2)R In(at.wt.) has been subtracted from the absolute values. This 
procedure was found by Powell and Latimer to improve the correlation, and 


TABLE II 
ABSOLUTE ENTROPIES AND UNIVALENT RADII 











Sake Univalent radius 

Ion (H+ = —5.5) 5°ars— (3/2) R In(at.wt.) 

Cst 26.3 11.8 1.69 
Rb* 24.2 9.0 1.48 
K+ 19.0 8.1 1.33 
Nat 8.9 — 0.4 0.95 
Agt 12.2 «a Ba 1.26 
Est —- 2.1 —- 7.9 0.60 
Ra?+ > — 14.1 BY ag 
Ba?+ —- 8 — 22.7 1.53 
Hg?+ — 16.4 — 32.2 1.25 
Cd?+ — 25.6 — 32.8 1.14 
Sr2+ — 20.4 — 33.7 1.32 
Ca?+ — 24.2 — 35.2 1.18 
Mg?* — 39.2 — 48.7 0.82 
Zn?t — 36.5 — 48.9 0.88 
TIs+ — 58.5 — 74.3 1.15 
In3* — 78.5 — 92.6 1.04 
Als+ — 91.4 —101.2 0.72 
Gast — 99.5 —112.1 0.81 
Ut — 100 — 116.3 1S” 
Pu‘t — 109 — 125.3 13° 





*Estimated by extrapolation. 


was found to do the same in the present study; its theoretical justification is 
considered below. Column 4 shows the “univalent’’ radii ry given by 
Pauling (11); these were found empirically to give better correlation than the 
ordinary crystal radii. This is understandable in view of the fact that the uni- 
valent radii represent the true sizes of the ions, whereas the crystal radii are 
reduced to various degrees owing to electrostatic attractions. 

In order to confirm the dependence of S°p, on 2? one may plot S°sps— 
(3/2) R In(at.wt.) against 1/r, and find by interpolation the entropies corre- 
sponding to some arbitrary radius. Taking 1.3 A as the radius, the entropies 
obtained are as follows: z = 1, 7.5 e.u.; 2 = 2, —27.0 e.u.; 2 = 3, —68 e.u.; 
and z = 4, —126 e.u. These values are shown plotted against both z and 2° 
in Fig. 2. The linear dependence on 2? is clearly established. 

Since Eq. [3] predicts a dependence of entropy on 2?/r, it is natural to test 
such a dependence first. Fig. 3 shows a plot of S°,,.—(3/2) R In(at.wt.) 
against 2?/r,. The linearity shown in this plot is as good as could be expected; 


*Anions are not considered in this paper, since the range of values is not sufficiently great to 
permit a test of the dependence on z and r. 
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corresponding to a univalent radius of 1.3 A. This plot establishes the fact that entropy varies 
linearly with 2? and not with z. 
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univalent radius. 


some improvement could be realized by adding a small constant value to all 
of the radii, but is hardly worth while in view of uncertainties in the univalent 
radii. It may be noted that small errors in the radii are highly magnified when 
z is large; an increase of 0.1 in the value for plutonium, for example, displaces 
the point by 0.9 to the left, and puts it to the left of the curve. 

The only point that is seriously off the line is that for aluminum. This may 
well be due to an error in the radius, and it is of interest that Goldschmidt (7) 
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gave a crystal radius (0.83 A) that is considerably greater than Pauling’s value 
(0.5 A). The use of a univalent radius calculated from Goldschmidt’s value 
would displace the point to the left of the line. Apart from this case 
Goldschmidt’s and Pauling’s values are mainly in good agreement; for lithium, 
however, Goldschmidt’s value is again larger than Pauling’s, and the change 
would displace the point away from the line. 

Attempts to find linear relationships using the crystal radii 7, rather than 
the univalent radii were not quite as successful. The addition of constant 
values to the crystal radii gave no advantage, and the use of functions like 
2/r.?, 2/(ryuta)’, 2°/r2, and z*/(r,+a)? showed no improvement over 2?/ry. 

It is therefore concluded that the function 2?/r, is the most satisfactory 
simple function. The equation of the straight line, obtained by the method of 
least squares (neglecting aluminum), is 


[4] S°ans = (3/2) R In(at.wt.)+10.2—11.6 22/ry. 

The mean deviation in entropy is 4.8 e.u., and the mean deviation in 1/r, is 
0.16; at a radius of 1.2 A this corresponds to a deviation in radius of about 
0.12 A, which is probably close to the actual uncertainty. 


DISCUSSION 


Several points of interest arise from this empirical equation. In the first 
place, the coefficient of the 2?/r, term in Eq. [4] is not far from that given by 
Eq. [3]. It thus appears that as far as entropies are concerned the simple 
electrostatic picture is much nearer to the truth than had hitherto been sus- 
pected. 

Other points relate to the significance of the term (3/2) R In(at.wt.) and to 
the intercept on the entropy axis. The contribution 


[5] S°n.e.s. = (3/2) R In(at.wt.)+10.2 


appearing in Eq. [4] is clearly the non-electrostatic part of the entropy. A 
very simple model (5) for a monatomic particle in the liquid phase regards it 
as having translational motion in a ‘‘free volume’”’ v; of the order of a molecular 
volume; the partition function corresponding to this type of motion is 


[6] f = {(2amkT)?/h3 }ve-F0/BT 


where Ep, is the difference in energy between the particle in the liquid phase 
and in the gas phase. If the ions are regarded as independent localized systems 
this partition function leads to the result that the entropy per mole is 


[7] S° = R[3/24+1n{ (2ermkT)*"v,/h?}] 
This is equivalent to 


[8] §°—(3/2) R In(at.wt.) = 3/2+(3/2) RIn{(2ekT)/2N}+R In x. 


Insertion into this equation of the numerical values (with T = 298.2° K.) 
gives rise to 
[9] §°—(3/2) R In(at.wt.) = 2.91+R In » 
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where 7; is the free volume in cu. A. 
In comparing this equation with Eq. [5] it is necessary to take into considera- 
tion the fact that Eq. [5] contains a ‘‘cratic’’ term (8) for the mixing of the 


ions with water; this term amounts to R In(1000/55.5) = 7.9 e.u. In order for 
Eqs. [5] and [9] to be equivalent it is therefore necessary that 


[10] 10.2—7.9 = 2.9+-R In 2. 


This leads to a value for 2; of 0.73 cu. A, which is entirely reasonable. The 
volume should, of course, be smaller than the volume of the actual cavities 
occupied by the ions, owing to the volumes of the ions themselves. 

The general conclusions to which the present work leads are therefore as 
follows. The absolute entropy of a monatomic ion in solution can be regarded 
as made up of a non-electrostatic and an electrostatic term. The first of these 
corresponds to the free motion of the ion in a free volume of about 0.73 cu. A, 
and the electrostatic term is represented approx /inately by Eq. [3], which is 
based on classical electrostatic theory. The view that entropies depend on z 
seems to be entirely without foundation. 

Later papers will deal more explicitly with the entropies of polyatomic ions. 
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THE ADSORPTION OF NITROGEN BY BROMINATED 
GRAPHITE AND CARBON BLACK! 


By H. L. McDERMOT AND J. C. ARNELL 


ABSTRACT 


The adsorption of nitrogen has been measured at 78° K. for three artificial 
graphites and one carbon black, which have been brominated to saturation and 
then debrominated by heat-treatment under vacuum or in the presence of an 
inert gas over a range of temperatures. During the bromine removal, the graphites 
were heated to a maximum temperature of 1000° C. and the carbon black to a 
maximum of 1200° C. When the graphite of smallest area was brominated it did 
not return to its original area even after heat-treatment at 1000° C., whereas the 
graphite of largest area exceeded its original area when it was brominated and 
heated at 1000° C. The surface area of the brominated carbon black was restored 
to the original value when it was evacuated at 400° C. After evacuation at 600° C., 
the surface area was greater than that of the original carbon black, but fell to a 

,Yalue less than the original area when the sample was evacuated at 1200° C. 
Chese observations are believed to be in accord with the known structures of the 
graphites and the carbon black. 


INTRODUCTION 

It has been known for a long time that graphite sorbs copious quantities of 
bromine when simply exposed to the vapor at room temperature. In 1933, 
Frenzel (2) found that-graphite took up 83.5% of its own weight of bromine 
from an atmosphere saturated with bromine vapor. This observation was 
confirmed by Riidorff (7) who undertook a systematic investigation of the 
sorption of bromine by a variety of natural graphites. He found that they all 
sorbed very nearly 83% of their own weight irrespective of the size of the 
graphite crystals. From this observation together with the large bromine 
sorption he concluded that bromine was not simply adsorbed on the surface 
of the crystals. Microscopic observation of single crystals of graphite revealed 
an expansion along the c-axis of the crystal accompanying bromination. X-ray 
diffraction photographs of the brominated graphite showed only 001 lines 
whose position indicated an identity period of 10.4 A compared to an identity 
period of 3.35 A for the untreated graphite. All of these experimental observa- 
tions were explained by postulating that the bromine was intercalated be- 
tween the sheets of hexagonally arranged carbon atoms which comprise the 
graphite crystal. Since the amount of bromine sorbed was independent of the 
size of the graphite crystals, it follows that virtually all the bromine was held 
in this fashion. 

Following the work of Riidorff, little work was done on the subject of 
brominated graphites. In particular, with the exception of an important 
investigation by Hennig (4), no work was done on the properties of the product 
remaining after debromination of the graphite. This product, which is a gra- 
phite containing a small amount of firmly held bromine, has been called a 
‘residue compound”’ by Hennig (4). He has shown that the composition of 

1 Manuscript received April 23, 1956. 
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these compounds is a function of the composition of the parent lamellar com- 
pounds and of the temperature to which the parent compound is heated during 
decomposition. Hennig measured the X-ray diffraction patterns of some of 
the residue compounds and found only a very slight increase in the c spacing 
of the graphite compared to that of the unbrominated graphite. He concluded 
from this that very little of the bromine of the residue compound was inter- 
calated in the graphite. He suggests that it was trapped by crystal imper- 
fections. 

The work to be presented logically follows that done by Hennig. It is a 
systematic investigation of some of the adsorptive properties of residue com- 
pounds prepared from three artificial graphites and one carbon black. In 
particular, the surface area and hysteresis effects have been studied by meas- 
uring the adsorption of nitrogen at 78° K. 


EXPERIMENTAL AND RESULTS 

Three graphites and one carbon black were chosen for study. The graphites 
used were artificial graphites manufactured by Acheson Colloids Ltd. and 
have been described previously by McDermot and Arnell (5). For the purpose 
of this work, they have been designated as G-1, G-2, and G-3. The carbon 
black, designated as Spheron 9, was provided through the courtesy of the 
Godfrey L. Cabot Co. Some of its properties have been given by McDermot 
and Lawton (6). All the samples were brominated by exposure to saturated 
bromine vapor in the absence of air at room temperature. The surface areas 
of the graphites and of the carbon black prior to bromination and the approxi- 
mate amounts of bromine sorbed by them are presented in Table I. 











TABLE I 
Surface area, Weight of bromine sorbed, 
meters?/gm. gm./gm. 
G-1 45 0.6 
G-2 96 — 
G-3 18 0.6 
Spheron 9 112 0.4 





It will be noted that there is no correlation between the weights of bromine 
sorbed and the surface areas of the graphites. Moreover the amounts of bromine 
sorbed are considerably less than those reported by Riidorff. 

After bromination the residue compounds were prepared by evacuating 
the samples at a series of temperatures beginning at room temperature and 
continuing up to about 1000° C. After each period of heating the samples were 
cooled in liquid nitrogen and the nitrogen adsorption measured at 78° K. In 
some of the early experiments, equipment was not available for heating the 
samples to high temperatures im vacuo and therefore they were heated in a 
stream of inert gas. For completeness the adsorption of nitrogen was measured 
for one graphite while it was still saturated with bromine. 
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The results for G-3 are presented first because most of the work was done 
with this graphite. The surface areas corresponding to different heat treat- 
ments are listed in Table II for three series of experiments carried out on this 
graphite. 











TABLE II 

Br retained, Surface area, 

Run Treatment gm./gm. meters?/gm. 
1 Saturated with Bre 0.6 5.2 
Evacuation at 100°C. — 10.5 
Evacuation at 300°C. 0.07 9.8 
2 Evacuation at 25°C. 0.068 8.6 
Evacuation at 100°C. 0.064 9.9 
Evacuation at 300°C. 0.055 10.1 
Evacuation at 600° C. 0 9.6 
Evacuation at 1000° C. 0 11.8 

3 Saturated with Brz to only 0.24 gm./gm. 

Evacuation at 25° C. 0.03 7.9 
Untreated graphite _ 18.4 





Some of the results from Run 1 are plotted in Fig. 1 along with the isotherms 
of the untreated graphite. The isotherm measured after the 300° C. evacu- 
ation is omitted because it was almost identical to that for the 100° C. evacu- 
ation. Fig. 2 depicts the isotherms measured after evacuation of the sample at 
room temperature and after heating it to 1000° C. 

In a previous paper dealing with these graphites (5), McDermot and Arnell 
demonstrated the existence of two types of hysteresis in the low temperature 
isotherms. The first type of hysteresis, which was believed to be due to the 
porous nature of the graphites, manifested itself as a broad hysteresis loop 
extending from saturation down to a relative pressure of approximately 0.45. 
This type of hysteresis has been called “‘pore’’ hysteresis. The second type of 
hysteresis, which was termed “swelling hysteresis’, was attributed to inter- 
crystalline swelling and was seen as an inability of the desorption curve to 
rejoin the initial adsorption curve below a relative pressure of 0.45, where 
“pore” hysteresis was manifestly impossible. 

Both of these hysteresis effects are evident in the isotherms for nitrogen 
adsorbed on G-3. In the case of the graphite saturated with bromine, the 
isotherm was found to be reversible over the whole range of relative pressures 
as might be expected if it is assumed that the pore space and the space re- 
sponsible for swelling are completely filled with bromine at saturation. Both 
types of hysteresis are absent. As bromine is removed from the saturated 
graphite the pores empty and the hysteresis loop, characteristic of the pore 
structure, reappears in the samples which have been evacuated at temperatures 
as low as 25° C. This is shown by the isotherm in Fig. 2 (triangles) measured 
after evacuation at room temperature. Swelling type hysteresis only reappears 
in the later stages of debromination where evacuation at high temperatures 
was used. In the case of G-3, swelling hysteresis was observed only after the 
sample had been heated to 1000° C. as is shown by the squares in Fig. 2. 
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An experiment was conducted to ascertain the effect of partial bromination 
of a sample of G-3. The sample was brominated to a concentration of 0.24 gm. 
bromine/gm. graphite and then evacuated at room temperature to a bromine 
content of 0.033 gm./gm. The isotherm measured after this treatment is com- 
pared to one measured after room temperature evacuation of a fully bro- 
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minated sample in Fig. 3. After evacuation at room temperature the bromine 
content of the saturated sample fell from 0.6 gm./gm. to 0.068 gm./gm. It is 
clear that the surface areas of the two samples are identical. It will be noted 
that the desorption curve obtained after evacuation of the partially bro- 
minated graphite is slightly higher than that measured after evacuation of the 
fully brominated sample. This is believed to be due to the fact that during 
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nitrogen adsorption a somewhat higher degree of saturation was reached in 
the former case than in the latter, rather than because of any difference in pore 
structure between the two samples. 

The effect of bromination on G-l1 is illustrated by the isotherms shown in 
Fig. 4. The isotherms were measured with a sample of G-1 which had been 
saturated with bromine, evacuated successively at room temperature and 
200° C., and then heat-treated at 500° C. and 1000° C. For the sake of clarity 
the desorption points at the higher pressures are omitted because this part of 
the isotherm is similar in shape for all of the graphites and is illustrated by the 
isotherm following the 200° C. evacuation. The lower desorption points are 
included because they again point out that the so-called swelling hysteresis, 
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which is absent in the room temperature isotherm, becomes apparent in the 
500° C. isotherm. It should be noted that unlike G-3, the surface area does 
return to almost its original value after the 1000° C. heating. The surface area 
corresponding to the various treatments along with the quantity of bromine 
retained by the graphite is given in Table III. 











TABLE III 
Treatment 

temperature, Br: retained, Surface area, 
yp gm./gm. meters?/gm. 

25 0.03 27 

200 0 32 

500 0 37 

1000 0 44 


Untreated graphite 46 
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Isotherms were measured with a sample of G-2 which had been saturated 
with bromine, evacuated at room temperature, and then heat-treated at 
500° C. and 1000° C. These isotherms along with that for the untreated gra- 
phite are presented in Fig. 5. For clarity many of the desorption points have 
been omitted. The most remarkable feature of these isotherms is the increase 
in surface area observed after the 1000° C. heat-treatment over that of the un- 
treated graphite. This is also brought out in Table iV which gives the surface 
areas corresponding to different heat-treatments. 
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TABLE IV 
Treatment 
temperature, Br2 retained, Surface area, 
a el gm./gm. meters?/gm. 
25 0.09 70 
500 0 92 j 
1000 0 121 
Untreated graphite 96 





The increase in area over that of the original graphite which was observed 
after the 1000° treatment was quite different from the behavior of the graphites 
with smaller surface areas. Accordingly an experiment was done to ascertain 
that the increase was indeed characteristic of the brominated graphite and not 
of the graphite. A sample of G-2 was heated to 1000° C. for the same length of 
time as was used for the brominated compound. The isotherm was measured 
for this sample and is shown in Fig. 6 along with that for the untreated G-2. 
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It can be seen that there is little difference between the two isotherms. 
Isotherms were measured for a sample of Spheron 9 which was brominated, 

evacuated successively at room temperature, 200°, and 400°, and then heat- 

treated at 600° and 1200° C. Representative isotherms are shown in Fig. 7 
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along with the isotherm of the untreated carbon black. All of the isotherms 
exhibited swelling hysteresis but this is only shown for one of the isotherms to 
avoid a confusion of points. It will be noted that this is the first carbon whose 
isotherms after room temperature evacuation have shown this type of hys- 
teresis. The surface areas corresponding to various temperatures along with 
the amounts of bromine retained by the Spheron are given in Table V. 

The table shows that the surface area increased to more than its original 
value at 600° C. and then with further heating to 1200° C. dropped to below 
its original value. 
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TABLE V 

Temperature, Bre retained, Surface area, 
we De gm./gm. meters?/gm. 

25 0.345 89 

200 0.091 98 

400 0.026 102 

600 0 131 

1200 0 94 

Untreated — 104 





DISCUSSION 


In the initial stages of debromination all the graphites behaved the same 
way. The adsorptive behavior of a fully brominated graphite is illustrated by 
the isotherm shown in Fig. 1 for a sample of G-3 saturated with bromine. It 
is clear,from the lack of hysteresis in this case that the bromine has filled the 
pores as well as the interstitial spaces which are believed to be responsible 
for swelling hysteresis. When the sample is evacuated at room temperature 
most of the bromine is removed, the surface area increases considerably, and 
part of the pore volume is emptied as is shown by restoration of the hysteresis 
loop. However, in all the graphites, the lower parts of the isotherms were still 
reversible after evacuation at room temperature showing that bromine con- 
tinued to block the adsorption sites that were responsible for swelling hyster- 
esis. McDermot and Arnell (5) have suggested that this so-called swelling 
hysteresis is due to a process of intercrystalline swelling. It is reasonable to 
suppose that the area between the crystallites responsible for this swelling was 
occupied by bromine and was thus denied to the nitrogen molecules. Hence 
the absence of swelling hysteresis and the cause for the loss of some of the 
surface area of the sample. 

In the case of Spheron 9 the surface area measured by the low temperature 
nitrogen adsorption method is not the area of the crystallites, but of agglo- 
merates of many crystallites. Moreover the crystallites of the carbon black 
are not regular crystals such as those of graphite but are more correctly parallel 
layer groups in which the orientation is random in an axis perpendicular to the 
layers of carbon atoms (1). Spheron 9 is not porous but it does display swelling 
type hysteresis (6). It is interesting to note in this connection that although 
the carbon black loses a considerable fraction of its original surface area on 
bromination and subsequent evacuation at room temperature, swelling 
hysteresis is still present in contrast to the behavior of the graphites. This 
point is shown very clearly by the isotherm outlined by the triangles in Fig. 7. 
Such behavior is remarkable because after evacuation at room temperature 
the carbon black retains more bromine than any of the graphites. McDermot 
and Lawton (6) have proposed that intercrystalline swelling is also responsible 
for the hysteresis effects that they have observed to accompany the adsorption 
of nitrogen by Spheron 9. Since the crystallites of the carbon black are com- 
pacted into agglomerates it is possible that spaces responsible for swelling 
hysteresis are of such a size that they may be penetrated by nitrogen molecules, 
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but not by the much larger bromine molecules. Hence, swelling hysteresis 
appears for the sample after evacuation at room temperature despite the fact 
that a relatively large amount of bromine is retained by the carbon black. 

It is of particular interest to examine in some detail the behavior of the 
graphites after evacuation at 1000° C. The graphite of lowest surface area G-3 
loses area when it is brominated and never regains it entirely even when heated 
to 1000° C. By contrast G-2, which has the highest surface area, shows an 
increase in area with respect to its original area when it is brominated and 
heated to 1000° C. G-1, whose surface area lies between the other two, just 
reaches its surface area after heating to 1000° C. The behavior of the graphites 
would thus appear to be related to the surface area of the starting material. 
Wilson and McDermot (8) have taken X-ray diffraction photographs of two 
of the Acheson graphites, of which one was G-3, and showed that the line 
broadening is consistent with their nitrogen surface areas. Hence the surface 
area can be taken to be a measure of the crystallite dimensions, in which case 
G-3 has the largest crystallites and G-2 the smallest. Viewed in these terms it is 
clear that the loss in surface area on bromination and subsequent removal of 
bromine at high temperatures such as was observed for G-3 can only be ex- 
plained by the filling of cracks or crevices by bromine or by crystallite growth. 
The first explanation would be expected to apply at low debromination 
temperatures because the bromine would be held by strong adsorptive forces 
which prevail in such fissures. However, it is reasonable to expect that a 
temperature of 1000° C. would suffice to drive off physically adsorbed bromine. 
Hence the loss of surface area after the high temperature treatment is likely 
due to crystallite growth resulting from the bromine and the high temperature 
treatments. In the case of G-2, the area of the brominated graphite which has 
been treated at 1000°C. is greater than that of the untreated graphite. By 
the same reasoning it must be assumed that the new surface area has been 
produced by a decrease in the average crystallite size due to an activation 
process resulting from the bromine and the heat-treatment. It has been proved 
experimentally that crystallite growth does not occur when G-2 is heated to 
1000° C. (see Fig. 6); therefore, the activation process is not opposed by a 
tendency toward crystallite growth. Hence the difference in behavior between 
the graphites must be due to different susceptibilities towards crystallite 
growth and towards activation by bromine. Furthermore it would appear that 
these susceptibilities are related to the crystallite dimensions because it is 
only in this respect that the graphites differ appreciably from one to another. 

It was observed that heat-treatment of brominated Spheron 9 at 600° C. 
produced a carbon black of much higher surface area than the original Spheron. 
This may have been caused by bromine activation at 600° C. which opened up 
new surface area for nitrogen adsorption. Further heating at 1200° C. di- 
minished the area to a value slightly less than that of the parent carbon black. 
In this instance, it is likely that the heat-treatment caused some crystal 
growth with a consequent loss of surface area. 

No attempt has been made in the course of this work to investigate the fully 
brominated compound in any detail. All of our attention has been focused on 
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the residue compounds left after the evacuation of a brominated carbon. 
However it is noteworthy that contrary to the findings of Riidorff (7), carbons 
of different crystallite dimensions sorbed different quantities of bromine. 
Riidorff (7) and Frenzel (2) found that all their graphites sorbed 0.80 to 
0.83 gm. of bromine/gm. of carbon. Goldsmith (3) reports a saturation value 
of 0.70 gm./gm. The values given in this paper range from 0.4 to 0.6 gm./gm. 
and do not bear any simple relationship to the area of the carbon. Thus it 
would appear that the quantity of bromine sorbed depends on the specific 
carbons and cannot be explained by any single process such as intercalation of 
bromine between the carbon planes. For example, in the graphites studied here, 
a part of the bromine was held in the pore structure and part was adsorbed on 
the surface. Judging by the ease with which most of the bromine could be 
removed, a large part was probably held by simple condensation between the 
particles. This idea is borne out by the observation that the powdered graphite, 
which was a free-flowing powder, became semisolid when saturated with 
bromine. This would be expected if the particles were stuck together by con- 
densed liquid. 

One experiment was done which showed that within the limits investigated 
the degree of bromination did not effect the surface area of the residue com- 
pound. A sample of G-3 which was only partially brominated yielded the same 
surface area after room temperature evacuation as a sample which had been 
saturated with bromine and then evacuated at room temperature, notwith- 
standing the fact that the quantities of bromine held by the two residue com- 
pounds were different. This is in accord with the notion that most of the 
bromine is held rather superficially and the only part that effects the surface 
area is that which is adsorbed or condensed between the crystallites. 
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INTENSITY IN THE RAMAN EFFECT 


Vv. THE EFFECT OF INTERMOLECULAR INTERACTION ON 
THE RAMAN SPECTRUM OF CARBON DISULPHIDE! 


By J. C. Evans? AND H. J. BERNSTEIN 


ABSTRACT 


The Raman spectrum of CS:2 has been obtained photoelectrically for the liquid 
and solutions of varying concentrations in cyclopentane, and depolarization 
ratios and intensities measured. The intensity of the forbidden bands v2 and »; 
relative to that of the Raman active band » decreases rapidly with dilution. 
The bands are present because of intermolecular interaction which probably 
occurs between the C atom of one molecule of CS: and only one of the S atoms 
of the other. 


The Raman spectrum of liquid CO, has been investigated recently (9) and 
the doubly degenerate bending mode v2 and the antisymmetric stretching 
mode »; which are forbidden by selection rules have been observed as weak 
bands along with the allowed totally symmetrical stretching mode »;. Although 
depolarization ratios of the forbidden bands were not measured it was argued 
from the band contours (9) in the liquid that » is depolarized and 3 polarized. 
These data have been interpreted on the basis of strong intermolecular inter- 
action in liquid CO, which distorts the molecule along its length so that the 
center of symmetry is essentially destroyed (9). 

One might expect a closely parallel situation with regard to spectral details 
and the effects of intermolecular interaction for carbon disulphide. In this 
paper we have investigated the Raman spectrum of CS,, liquid and in solutions 
of cyclopentane, in order to establish—(i) that the presence of the forbidden 
bands in the Raman spectrum is indeed due to the environment and (ii) the 
effect of the interaction on the shape of the molecule. 

The Raman spectrum of CS, (5), depolarization ratios (2), and relative 
intensities (6) have been obtained some years ago and the discussion of its 
spectrum and the effect of Fermi interaction given by Placzek (7) on the basis 
of these older data. Newer data have been compiled by Herzberg (4) and the 
positions of the bands given there are confirmed by our investigation to within 
+1 cm.—'. For purposes of calculating resonance interaction constants the 
more accurately measured data quoted by Herzberg have been used. The 
depolarization ratios of the forbidden bands » and v3; have been measured 
confirming Bhagavantam’s result (2) for v2 and »; (polarized). For purposes of 
comparing frequencies the infrared spectrum of the liquid was also obtained. 


EXPERIMENTAL 


The Raman spectra were photoelectrically recorded on a White grating 
Raman spectrometer (10) using 44358 as the exciting line. Depolarization 
ratios (pop.) Were measured by wrapping polaroid film around the Raman tube, 
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first with its axis parallel to the tube and then at right angles, and these ratios 
corrected for convergence (8) to give peor. Integrated intensities were measured 
relative to the 458 cm.— band of CCl, obtained under identical instrumental 
and optical conditions and were corrected for the spectral sensitivity of the 
photomultiplier o4, and for convergence error according to [1] 


las Tops(1+Peorr) , / {eaten \ 
(] Terue =| 1+ Pods ie 1+ pons ”" 458cem. —1 


where 100 is taken as the true integrated intensity of the 458 cm.—! band 
in CCl. 

When bands overlapped, no attempt to resolve them by curve analysis was 
made, the measured intensity being taken as the area under the contour. This 
is indicated by bracketing bands in the tables. The Raman spectra of three 
solutions of CS, in cyclopentane were also measured in which the ratio of the 
volume of CS, to the total volume was 0.65, 0.49, and 0.36. The infrared 
spectrtm was obtained with a Perkin Elmer 12C double pass spectrometer 
using CsBr, KBr, and NaCl optics. The results for CS, liquid are given in 
Table I and the Raman spectrum shown in Fig. 1. 


]x100 














TABLE I 
Liguip CS. 
Raman 
Infrared 
Assignment? Av? Av Pobs Peorr I true? v 
v2 390 0.75 0.63° 9.7; 393 Ss. 
vy’ 648.3 648) 0.41 0.2; 488 
v1 656.5 655 654 w. 
2v2 796.0 797 | 0.40 0.2; 94.2 789 v.w. 
Qv2' 804.9 804 / 809 v.w. 
v3 1515 0.79 0.63 6.94 Not 
measured 





4y,’ and 2v2' are satellite lines corresponding to the transitions (11'0), (03’0) — (01’0). 
’Data from Herzberg, Ref. 4. 

¢Measured by Bhagavantam to be 0.8 (Indian J. Phys. 7:83. 19382). 

4The 458 cm. band of CCl, ts 100 


Rough decomposition of the overlapped doublets at ~650 and 800 cm.~! 
indicates that the band at 648 cm.~! contributes about one-third of the total 
intensity and 804 cm.—! contributes about two-fifths of the intensity in its 
doublet. The ratio of the areas under the doublets is then 


Ts00/Ies0 = 94.2/488 = 0.193 


and a rough estimate of J2,,/J,, is 0.174. These ratios are to be compared with 
0.32 obtained photographically by Pienkowski (6). 

The intensities of the doublets at 650, 800 cm.—! and the forbidden band at 
390 cm.—! have been measured as functions of concentration and the results 
shown in Figs. 2, 3, and 4, respectively. It should be mentioned that depolar- 
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ization ratios of the doublets and frequencies did not change within experi- 
mental error as a function of concentration, nor did the contours of the doublets 
at 650 and 800 cm.—!. The » band position changed slightly as a function of 
concentration (from 390 cm.— in CS, to 387 cm. in dilute solution) whereas 
the depolarization ratio was too uncertain to indicate whether or not changes 
in it took place. The depolarization ratio under most favorable conditions was 
very difficult to measure for the 390 cm.—' band because of overlapping (at 
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wide slits) with either the exciting line or the band at 650 cm.~! or both. 
The uncertainty in the quoted depolarization ratio may be as much as 0.15 
for v. whereas it is about 0.03 for the other bands. 


DISCUSSION 


From Placzek’s formulae (7) for the separation of the main resonating com- 
ponents » and 2m, viz. X; = +/(6+4w*), and that for the satellites (»,’ and 
2v.'), viz. X2 = +/(&+8w*), we find, using the more accurate data quoted by 
Herzberg, that 6 = 120.0 cm.—' and w = 35.6 cm.~! for liquid CS». Further, 
according to the theory for resonance (5), 


Lavy = (Mas! VX +5 " 
- VX14+6 ayt+ V2 1—6 Qo 
where a; and a are the polarizability derivatives for the transitions 


(10°0) — (00°0) and (02°0) — (00°0) 





respectively. If a2/a; <1, which is usually true for the ratio of overtone to 
fundamental, 
Toy./Iy, = (X1—8)/(X1+6) = 0.0755. 


The calculated value for this intensity ratio is about half the observed value. 
Using the approximate value for J:,,/J,, obtained experimentally (0.174) and 
solving for az/a: gives —0.128 for this ratio. Since the uncertainty in the 
experimental intensity ratio is large, it is not even possible to say that the sign 
of the calculated ratio a2/a; is correct. The negative sign would imply of course 
that either a2 or a, was negative. Since the sum of the unperturbed levels 
°+2y° is the same as 1.+2v. = 1452.5 cm.~', and 2»,°—»° = 120.0 cm.—, 
we find that 27° = 786.2 cm.—' and »;° = 666.2 cm.—! for liquid CS». 

From Figs. 2, 3, and 4 we see that the bands at 650 and 800 cm.—' do not 
deviate as markedly from the linear concentration dependence as does the 
band at 390 cm.—!. The deviation from linearity can arise not only from the 
effects of intermolecular interaction but also because of the dependence of 
the intensity on the refractive index and the reflection loss (1). If account is 
taken of these factors and the corrections indicated in Ref. 1 made, the result 
that the most marked intensity fall-off with concentration is experienced by 
the 390 cm.~! band is not altered. This would be expected if the presence of 
this band is due to intermolecular interaction. It is interesting to compare 
also the relative intensities of these bands as a function of concentration (see 
Table IT). 

It is evident from Table II that the relative intensity of the allowed doublets 
at 650 and 800 cm.—' remains approximately constant whereas the relative 
intensity of the 390 cm.—! band approaches zero with dilution. This is direct 
evidence that this band arises primarily from the action of the intermolecular 
force field. Buchheim (3) has shown that the ratio of the intensity of this for- 
bidden band to that of the allowed band » should be proportional to the 
square of the field strength. If we assume a coulombic field between interact- 
ing CS molecules, and that the distance between interacting molecules is 
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TABLE II 
Relative intensity* 
Solution, 
vol. fraction 390 cm.~! 650 cm.~} 800 cm.—! 
of CS. doublet doublet 
1.0 2.00 100 19.3 
0.65 1.1; 100 18.1 
0.49 0.69 100 18.6 
0.36 0.50 100 17.0 





*The fourth power scattering term and the temperature dependent 
exponential term have not been included in the intensity. 
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Fic. 5. Relative intensity of the 390 cm.~! band of CS: asa function of (concentration) /3. 


proportional to the (molecular-volume)'/* then according to Buchheim the 
relative intensity of the 390 cm.—! band should be proportional to (volume- 
fraction)‘/*. This plot is shown in Fig. 5. 

The reason the forbidden bands », and »3; appear in the Raman spectrum is 
to be found in the effects of the interaction field. In CO, the field distorted the 
molecule in the direction of the OCO axis giving essentially C,,, symmetry (9) 
with its allowed Raman spectrum of two polarized (», v3) and one depolarized 
line (v2). In the case of CS. however both »; and » are found to be polarized. 
This suggests that the distortion of the CS: molecule is such as to leave the 
molecule with a plane of symmetry only (C;), in which case all three funda- 
mentals are allowed in the Raman effect and should be polarized. This would 
imply not only distortion along the SCS axis but also at right angles to it.This 
would correspond to the interaction between CS, molecules taking place be- 
tween a C atom of one molecule and only one of the S atoms of the other so 
that a real difference in field is experienced by the S atoms. Further, this 
difference is not only along the CS bond but also at right angles to it. On the 
basis of these spectroscopic data for liquid CS, and CO: one might expect the 
crystal structure of CSz to be different from that of COQy. 
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THE DIFFERENTIAL EQUATIONS OF CONSECUTIVE REACTIONS' 


By A. E. R. WESTMAN AND D. B. DELuRY 


ABSTRACT 


The differential equations of consecutive first-order reactions are formulated 
in a general manner. The solutions of these equations are shown to fall into a 
pattern, according to which the solutions to all such systems can be written down 
without resort to formal integration. The pattern is also convenient for making 
numerical calculations. 


A SPECIAL CASE 


Suppose that a molecule of tetraphosphate breaks down under hydrolysis 
to form one molecule of triphosphate and one of orthophosphate, that tri- 
phosphate breaks down into one molecule of pyrophosphate and one of ortho- 
phosphate, and that one molecule of pyrophosphate breaks down into two of 
orthophosphate. Let the symbols ./ 4 @ & refer to these phosphates and 
let a, 8, y, 6 stand for the numbers of phosphorous atoms in them (a@ = 4, 
8 = 3,7 = 2,6 = 1). Leta, b,c,d (a+b+c+d = 1) represent the proportions 
of the phosphorus atoms, at time ¢, which are present as ‘Y% -4 @ FZ and let 
do, bo, co, do represent their values when ¢ = 0. Then, if the reactions are first 
order and independent of each other and if their rate constants A, B, C* are 
known, the proportions a, }, c, d at any time ¢ are defined by the following 
system of differential equations: 


(1) —da/dt = Aa, 
[I] (2) —db/dt = Bb—(8/a)Aa, 
(3) —dc/dt = Cc—(y/8)Bb, 


with initial conditions a = ao, b = bo, ¢ = Co, d = dy when t = 0. The value 
of d can be specified either by a differential equation or more simply by the 
relation 
d = 1—a—b—c. 
With a view to generalization later on, the symbol (AZ) will be introduced 
here to stand for the fraction of the phosphorus lost to “that goes directly 
into 2. In terms of these symbols, the reaction described above is specified by 


(AB) = B/a, (SC) =0, (YD) = 1—(B/a); 
(BE) = y/B, (BL) = 1—(y/8); 
(GP) =1. 


The ratios B/a and 7/8 in Equations [I] will be replaced by (7:4) and (2@). 
The explicit solution of a set of equations like [I] offers no mathematical 
difficulty, but the resulting formulae are e!aborate and cumbersome and offer 
too many opportunities for mistakes in calculation. The purposes of this note 
1Manuscript received April 27, 1956. 
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are to put forward an orderly notation for writing the solutions and to point 
out that, once the underlying pattern of the solutions is perceived, the solutions 
may be written down directly without resort to formal integration of the 
differential equations. 

Direct integration of Equation (1) yields the solution 


a = ae*', 
which will be rewritten in the form 
@ = —a,(l—e-4*) +a. 


Again, integration of Equation (2) produces the formula 


6 = 4 


which will be ee in the form 





J G 
St, 








Bo +A) cage Lb nett 


b = ad ( tb B)s we —(G B) ste) |-sa-e? +65 


Solutions for c, d, ... can be expressed in the same manner, but the formulae 
soon become extremely elaborate. It is suggested here that a form of matrix 
notation simplifies the statements of these formulae and provides a pattern 
for writing directly the solutions to any such system of equations. One con- 
venient scheme is the following. Write in a column the arguments 1—e~4‘, 
1—e~*', etc., and finally, 1. Write in a row the initial amounts do, do, ... , as 
shown in Table I. Then the cell of the table corresponding to an argument and 
an initial amount receives the appropriate multiplier. 

The full pattern of the solution can be seen in the formula for c. The column 
of multipliers corresponding to ap all have the coefficient (:7F)(BP) which 
reflects the route by which ay has contributed to c. Likewise, they all contain 
the product of the rate constants AB, the two rates which are associated with 
this path. The denominators are selections from the differences A —B, A—C, 


TABLE I 











a 1—e4!? -1 





r A 
1—e~4t +(AB)5—, 





bine (AB) Zs ” 
; 1 
a AB 
let! — (AB) BOG ao 
AB C Dp B 
c 1-e8  —  +( AB) BO) GR Bo HBO) B= 
BE AB JAG B 
1-e tt (AB) BO) GH Bo (BO) B=! 
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B-—C, all taken according to the order of the sequence. The denominator in 
the multiplier of 1 —e~4' is the product of those differences in which A appears, 
that of 1—e~*' is the product of those in which B appears, and so on. Lastly, 
the signs are alternately plus and minus, starting with a minus sign at the foot 
of the column. Similar observations apply to the d) column. The final column 
always consists simply of the —1 and 1. 

Application of these rules provides the solutions to any such system of 
equations provided the conditions (-/@) = 0, etc. are met. A useful check 
when the multipliers are being calculated numerically is provided by the fact 
that the multipliers must sum to zero in each cell. This is easily demonstrated 
algebraically and corresponds to the fact that, as ¢ goes to infinity, the pro- 
portions of <4 ZB © etc. go to zero. (This must be modified if a rate constant is 
zero. See Appendix I1.) 

If, as an alternative to the assumption which led to System I, it is assumed 
that tetraphosphate hydrolyzes directly into pyrophosphate, which in turn 
breaks down into orthophosphate, System I must be replaced by another set 
of differential equations whose solutions are not included among those of 
System I. The solutions to the new system can, of course, be written out, but 
it seems worth while to lay out a general system of differential equations which 
includes as special cases those already mentioned and many others and to 
develop a rule for writing the solutions to the general system. 


THE GENERAL SYSTEM 


A reaction in which each of a sequence of substances -/ 4@ FD .. . breaks 
down to form some or all of the substances which follow it in the sequence, 
but not those which precede it, will be termed a consecutive reaction. In terms 
of the symbols and assumptions introduced earlier, the differential equations 
for this reaction are 


(1) —da/dt = Aa, 
[11] (2) —db/dt = Bbh—(.%A)Aa, 
(3) -—dce/dt = Cc—(%@)Aa—(BC)B», 


and so on. 

The pattern for writing the solutions to System I] is easily built up from the 
rules developed for System I. Equations (1) and (2) are the same in both 
systems, but Equation (3) in System II includes a term which provides for 
the path ./— ©. The contribution from this path is provided for in the solu- 
tion by the rule given earlier for a path like /— &, with F# replaced by 
and B by C. The ao column in the formula for ¢ is shown in Table II. 

Once the solutions to the general system are written out, they may be re- 
duced to any special case by giving the appropriate values to the symbols 
(4B), etc. and, when necessary, to the rate constants. For example, the 
situation tetra- — pyro- — ortho-, tri- — pyro-+ortho-, mentioned above, 
is specified by the values 

(YB) = 0, (AG) = 1, (AD) = 0; 
(BC) = 4/8, (BD) =1-(y/8); 


(6D) =1., 
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TABLE II 
ao 
_ AB afi sil 
1—e-At —(%B)( BO) a= aaa t(Y%O)g—e 
ws AB 
c l-e 8 +(AB)( BO) GR BIO 
“t oJ G BG\—AB__ SC) 
1 —¢-¢ -(%B)( BC) ae (Bao) ~ (YO) A 





Again, if it were assumed that, of the portion of /which reacts, a fraction p 
follows the path ./— & and the fraction 1— p follows the path /— @ the 
system is specified by 


(YB) = (8/a)p, (%C)=(1—-p), (AD) = p[l—(8/a)] ; 
(BC) = y/B, (BD) = 1-y/B; 
(CFD) =1. 


APPENDIX I 
The complete formula for d is given here to display the full pattern of the 
solutions to the general system. d is determined by the differential equation 


—dd/dt = Dd—(%D)Aa—(BD)Bb—-(CTD) Cz, 


where a, 0, c are given by Equations (1), (2), (3) of System II. 


APPENDIX II 


The general formula for d is here specialized to describe the case in which F 
is the final substance in the sequence. That is, does not break down (D = 0) 
and all the other substances break down, by all possible routes, into Z and 
nothing else. It follows that 


(CDP) =1, (B@)4+(BD) =1 (AB)4+(AC)4+(AD) = 1. 


The solution corresponding to these conditions is obtained from the general 
formula given in Appendix I simply by substituting D = 0. The form of the 
solution is altered somewhat by this substitution, because the argument 
1 —e—?! vanishes and therefore the corresponding row of multipliers is dropped. 
The resulting formula for d has the property that the sum of the coefficients 
in each cell is unity. This corresponds to the fact that, when ¢t becomes infinite, 


d = aptbho+co+do = 1. 
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THE INFRARED AND RAMAN SPECTRA OF ALLENE-1,1-d, ! 


By J. C. Evans,? J. K. WILMsHuRsT,? 
AND H. J. BERNSTEIN 


The infrared spectrum of allene-1,1-d, vapor has been observed from 2 to 
35 uw (Fig. 1) using a Perkin Elmer Model 12C double pass spectrometer 
equipped with LiF, NaCl, and CsBr optics and the Raman spectrum obtained 
for the liquid at — 50° C. (Fig. 2) with a White grating Raman spectrometer (8) 
with photoelectric recording. The depolarization ratios of the Raman lines were 
obtained by the method of Edsall and Wilson (3). 


1.0 





0.8} 
0.6} ~40°C -40°C 
04+ -40°C 


0.2) 








1.2 [2000 21002200 2300 2400 2500 2600 2700. 2800 ~=2900 ~— 3000 3100 3200 3300) 





Lob 
08+ “aac 
0.6}- 
0.4 


0.2 








Pa ere §=6RSIORY. = Sa (700° x: 


OPTICAL DENSITY 


0.8 


o4L 


0.2 








600 300 1000 i ns Fe 














0.2 Pan” 
i be 1 én 1 


360 400 500 — 700 





FREQUENCY cm-' 
Fic. 1. The infrared spectrum of HxXC=C==CD: vapor in a 10 cm. cell at the vapor pressure 
of the liquid for the temperature shown; (a) 6 cm. cell, liquid at —50°C. 


1{ssued as N.R.C. No. 4011. 
2National Research Laboratories Postdoctorate Fellow 1953-19565. 
3National Research Laboratories Postdoctorate Fellow 1955- . 
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Fic. 2. The Raman spectrum of liquid HxC—=C—=CDz (slits ~7 cm.~!). The intensity of 
the band at ~ 3000 cm.~! is taken as 100. dp = depolarized, p = polarized. 


TABLE I 
FREQUENCY ASSIGNMENT FOR ALLENE-1,1-d2 








Approximate 








Va Cay description CH=CH CH.—C=CD: CD —C—eDp,. 
a a, v, C—H stretching 2996 2995 R 2195 
ay a, v2 CHe bending 1440 1169 R 1228 
ay a, v3 C=C=C stretching 1076 924 R 874 
be a, vs C—H stretching 3005 2226 IR 2230 
be a, ve C=C=C stretching 1957 1941 IR 1921 
be a, vw CHe bending 1398 1416 IR 1034 
by do vs Torsion 865 (750) 615 
e bi,b2 ~=vs C—H stretching 3085 3105 IR, 2309 R 2330 
e bib. va CHerocking 1015 1058 R, 842 IR 830 
e€ bib2 ~vio CHe rocking 842 676 R, 842IR 667 
e bi,b2 v1, C==C=C bending 354 317 R, 317 R 306 
R = Raman band. 


IR “ infrared band. 
(750) = calculated from the Teller Redlich product rule (5). 


An assignment of all the fundamentals for allene-1,1-d, has been made on 
the basis of C.» symmetry and this assignment is given in Table I together 
with that of Lord and Venkateswarlu (6) for allene and allene-d,. 

The type a; vibrations in allene-1,1-d2, corresponding to the a; and be type 
vibrations in allene, give rise to polarized Raman bands and parallel type 
infrared bands and can be assigned immediately. The type 5; and » vibrations, 
corresponding to the type e doubly degenerate vibrations of allene, should give 
rise to depolarized Raman bands and perpendicular type infrared bands and 
are assigned by comparison with allene and allene-d,. A partial resolution of 
the two perpendicular bands at 842 and 1058 cm.—' was obtained. The az 
torsional frequency was not observed but was computed from the Teller 
Redlich product rule (5). All the remaining observed frequencies can be easily 
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assigned, in terms of these fundamentals, to overtone or combination bands 
(Table II). The Raman band at 1223 cm.— is probably »; (2995 cm.—) excited 
by the 4047 A mercury line. 

The Teller Redlich product rule (5), Decius-Wilson sum rule (2), and 
Bernstein—Pullin sum rule (1) are consistent with the present assignment 
(Table ITI). 

On attempting a direct comparison of the a; and b» frequencies in allene with 
the a; frequencies in allene-1,1-d2, it is found that an interesting change has 


TABLE II 
OVERTONE AND COMBINATION BANDS FOR ALLENE-1,1-d2 











Obs. frequency Assignment Calc. 
498* IR v9(b2) —v11(b2) 525 
626 IR Vil 634 

1158 R v9(b2) +r11(b1) 1159 
~1285 IR vs(b2) —[1008] 1301 
1334 IR 2Xvio 1352 
1433 R v10(b1) +4(a2) 1426 

~1566 IR v4(d2) +v9(b2) 1592 
1686 IR 2 Xvo(b2) 1684 
1935 IR v6(@1) to11 — vit 1941 
1977 R ? 

2173 IR 2 [1086] 2172 
2615 IR vg(b2) +y11(b2) 2626 
2812 IR 2 Xv7(a1) 2832 
2860 IR 2 X[1086 | +»10(d1) 2848 
2949 IR v6(41) +[1008] 2949 





[1086] and [1008] are the estimated positions of the 
unperturbed fundamentals v2(a,) and v3(a;). 

*May be due to CH;C=CD impurity which has its 
strongest infrared band at 498 cm.~ (4). 

tConsistent with corresponding assignments in allene 
and allene-d, (6). 


TABLE III 
PRODUCT RULE RATIOS 








Symmetry species used for ratio calculation 























Theor. Obs. CH »CCH. CH,CCD.2 CD.CCD, 
0.513 0.519 ay, by a 
0.512 0.527 a, a1, be 
0.866 _ by ae 
0. 816 seasons ae by 
0.382 0.426 e bi, be 
0.461 0.457 by, be e 
APPLICATION OF SUM RULES TO ALLENE-I,1-d2 
=v?, cm.~? =v, cm.~! 
Species Theor. Obs. Cale. Obs. 
a, 2.192 X 107 2.192 X 107 10,677 10,671 
bi, be 1.792 X10" 1.817 X10" 9,415 9,466 
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occurred and two bands in the 1000-1400 cm.~! region, one due to the sym- 
metrical CH» bending and the other due. to the symmetrical C—=C=C 
stretching frequency, have undergone major changes upon deuterium substi- 
tution. This change is explained by the coupling of the symmetrical CD. 
bending mode with the symmetrical C—C=C stretching mode, which give 
rise to frequencies in the 1000-1100 cm.~'! region, causing one frequency to 
increase to 1169 cm.~', and the other to drop to 924 cm.~'. This coupling is also 
present in allene-d,, but to a slightly greater extent, the two resultant fre- 
quencies being at 1228 cm.—! and 874 cm."!. 

The sample of allene-1,1-d, was prepared for us through the courtesy of 
Morse and Leitch (7). 
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A NOTE ON THE SYSTEM WOLFRAM-TELLURIUM 
By OsvAL>D KNoPp AND HAAKON HARALDSEN 


Goldschmidt (2) considered wolfram ditelluride to be isostructural with 
wolfram disulphide and wolfram diselenide, but this has never been confirmed 
by a structural investigation. In fact, in a preliminary study of the binary 
systems of wolfram with sulphur, selenium, and tellurium, Stangebye (7) 
showed that the compound WTe,, prepared by heating stoichiometric amounts 
of the elements im vacuo, did not give X-ray powder diffraction photographs 
similar to those of the disulphide and diselenide. 

This observation was fully confirmed in the present investigation. The rota- 
tion and Weissenberg photographs of a single crystal found in the mass of fine 
ditelluride powder showed that the compound had an orthorhombic structure 
with the unit-cell dimensions a = 3.490+0.006 A, b = 6.277+0.015 A, and 
c = 14.07+0.02 A. The volume of the unit cell based on these dimensions, 
308+2 A’, and the experimental density, qd” = 9.40, led to four WTez units 
per unit cell. Reflections of the type (h0/) were absent when h+/ = 2n+1. 
The space groups consistent with these systematic absences are D),—Pmnm, 
& — Pmn, and D3} — P2122. 

Unless WTez is polymorphic—and we obtained no evidence to that effect— 
it cannot be considered isostructural with the C7-type compounds, i.e., with 
MoS:, MoSez, WS2, and WSeo. There are indications, however, that it has a 
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layer-like structure closely related to the C7 type. The axial ratio of WTes, 
c/a = 4.03, does not differ greatly from the values 3.92-3.93 and 3.94 for the 
corresponding sulphur and selenium compounds (9, 10), and a value for 6 
about 3.7% lower than the above value would give an axial ratio b/a = ./3 
required for an orthohexagonal cell. This suggests a modified double unit cell 
of the C7 type. A Patterson projection on the (010) plane seemed to confirm 
this similarity but the intensity data on which it was based were rather 
inaccurate. The single crystal used was a very thin plate of irregular shape, 
with the approximate dimensions 0.15, 0.08, and 0.02 mm. along the a-, b-, 
and c-axes, respectively, and there must have been heavy irregular absorption 
of X-rays in the crystal. Powder photographs taken both with an ordinary 
and with a microfocus camera contained too many diffuse and overlapping 
lines to permit an intensity evaluation. Attempts are being made to secure 
more suitable single crystals of WTe, that will provide reliable data for a 
structure determination. 

No line shift outside the usual measurement spread was detected in the 
powder photographs of samples having the gross compositions WTe;.9. and 
WTe;. These samples were prepared in the same manner as WTee and were 
annealed in vacuo at 500°C. for several days. It seems likely that the homo- 
geneity range of wolfram ditelluride is very narrow (compare the slight vari- 
ation in the lattice constants reported for WS, (9)). This indicates that if the 
composition of the single crystal used deviated at all from the stoichiometric 
composition, it could have been only by a very small amount. 

In agreement with the results of Stangebye (7) and Morette (4, 5) no com- 
pound other than WTe, was detected in alloys obtained by direct synthesis in 
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Fic. 1. Phase relationships in the system wolfram-tellurium under orthobaric conditions. 
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vacuo. Combined evidence from X-ray powder photography, and from thermal 
analysis of samples sealed in small evacuated. quartz containers, suggests 
phase relationships as represented in Fig. 1, with the system under its own 
pressure. However, one cannot exclude the possibility that other compounds 
might be obtained by low-temperature annealing of intimate mixtures of WTe2 
and tellurium or wolfram for long periods of time. As is evident from Morette’s 
examination of polished sections of alloys rich in tellurium (5), the solubility 
of wolfram in tellurium at temperatures below 500°C. must be quite low, as 
a solidified alloy with 0.33 at. % W still contained crystals of the ditelluride 
embedded in a matrix of tellurium. 

A fuller account of the continued investigation, and of a similar investiga- 
tion in the system molybdenum-tellurium now in progress, will be published 
later. 


EXPERIMENTAL 


The wolfram used in this investigation was prepared by reducing calcined 
WO; (c.p. grade) with hydrogen at temperatures between 650° and 800°C. 
The lattice constant of the resulting lightly-sintered powder was found to be 
3.1647+0.0003 A at room temperature, in very good agreement with the 
value 3.16477 +0.00008 A reported by Cohen (1). Both reduced metal and a 
commercial wolfram powder reported to contain 99.8% W and 0.2% Mo were 
used in the preparation of samples with the composition WTe: oo. No differ- 
ences in the properties of the resulting powders were observed. 

A grant of tellurium of about 99.9% purity was kindly made by the Research 
and Development Department of the Canadian Copper Refiners, Limited, 
Montreal. According to the analysis supplied the main impurity was selenium. 
The lattice constants of the material were determined as 4.456+0.002 A 
and 5.925+0.002 A, in satisfactory agreement with the values of Straumanis 
(8), 4.45591+0.00012 A and 5.92687-+0.00002 A. 

Accurately weighed amounts of the elements totalling 15-20 gm. were intro- 
duced into degassed clear-quartz tubes which were then evacuated to 10-§ mm. 
Hg, sealed off, and put into a deep vertical tubular furnace at 700-800°C. 
After several hours the power was turned off and the samples allowed to cool 
in the furnace. The alloys were then crushed and ground to 100 mesh, thor- 
oughly mixed, and annealed at 500°C. for several days in sealed evacuated 
quartz tubes. Samples of the gross compositions WTe1.50, WTe1.90, WTe2.00, 
WTe;, WTes, WTes, and WTe,; were prepared. 

The density of the ditelluride was determined pycnometrically with water 
as displacement liquid, about 1 cm.* of powder being used in separate deter- 
minations. The average of several determinations gave d*? = 9.40. Morette 
(4, 5) reported d'§ = 9.44+0.02, without mentioning details of his method, 
while the value calculated from the above lattice constants is 9.46+0.06 
gm./cm.? 

Although there was practically no loss of material during the handling of 
the powders, all samples were analyzed. The determination of tellurium using 
sulphur dioxide and hydrazine hydrate (3) gave satisfactory results. Wolfram 
was determined as WO;. The analytical results were in close agreement with 
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the intended gross compositions. Indeed, any discrepancies between the ana- 
lytical results and the compositions of the samples as obtained by weighing 
are likely to have originated in the inaccuracies inherent in the analytical 
method. 

The samples for thermal analysis were sealed in evacuated containers of 
clear quartz having internal free volumes of about 3 cm.’ An average of about 
5 gm. of sample was used in individual runs. There was a pronounced tendency 
to supercooling, and the geometry of the heating curves depended considerably 
on the heating rates. 

Filtered or unfiltered CuK-radiation (Aa; = 1.53740 kX, Aa, = 1.54123 kX) 
from a hot-cathode tube was used throughout. The powder photographs were 
taken in a Straumanis-type camera with a nominal diameter of 114.6 mm., 
and the rotation and equi-inclination Weissenberg photographs of the single 
crystal about the a- and b-axes (three and four zones, respectively) were taken 
in a Unicam instrument, using the multiple-film technique. The intensities 
were estimated visually. Nelson and Riley’s method (6) was used in the treat- 
ment of the lattice constants. Difficulties were encountered in the measurement 
of the films, as WTe: did not give good powder photographs. 
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